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Page  7. 

Preface. 

New  physical  concepts  are  created  not  only  in  the  process  of 
generalizing  the  physical  theories,  but  also  by  the  return  route: 
they  can  arise  as  a  result  of  applying  the  approximation  methods  to  a 
more  precise  physical  theory.  On  the  fundamental  value  of 
approximation  methods  in  theoretical  physics  for  us  it  was  already 
necessary  to  write  in  connection  with  quantum  mechanics  [44].  Since 
this  book  is  dedicated  to  electromagnetic  theory,  we  will  examine 
here  examples  of  the  emergence  of  new  concepts  in  the  theory  of 
light/world  (both  in  the  electromagnetic  theory  and  in  the  simpler 
wave  theories  of  light/world  preceded  it). 

So,  the  concept  of  ray/beam,  and  equally  all  geometric 
optic/optics  can  be  derived  from  the  wave  theory  of  light/world  as 
the  idealizations,  suitable  in  the  extreme  case  of  very  small 
wavelength  (in  the  region  near  the  boundary'  of  light/world  and  shadow 
these  idealizations  are  already  unsuitable).  During  the  less  complete 
idealization  are  considered  the  divergences  from  the  geometric 
optic/optics,  in  other  words,  it  is  considered  the  diffraction,  such 
as  it  is  also  the  new  physical  concept  (diffractive  phenomena  are 
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most  clearly  seen  exactly  near  the  boundary  between  the  light/world 
and  the  shadow) . 

In  the  first  investigations  of  diffractive  phenomena  the 
specific  properties  of  the  material  of  the  diffracting  body  usually 
were  not  considered:  the  body  was  accepted  as  the  absolutely  black 
(absorbing).  Laws  of  reflection  which  consider  these  properties,  in 
the  analysis  of  diffractive  phenomena  were  not  utilized.  As  the  basis 
of  the  theory  of  diffraction  was  assumed/set  Huygens-Presnel ’ s 
principle;  on  this  basis  light  field  near  the  boundary  of  light/world 
and  shadow  was  described  by  means  of  Fresnel's  integrals,  and  field 
near  caustics  -  by  means  of  Airy's  integrals. 

With  the  emergence,  about  hundred  years  ago,  electromagnetic 
theory  of  light  (applied  also  to  the  radio  waves)  it  became  necessary 
to  anew  formulate  the  theory  of  diffraction.  As  reliable  theoretical 
basis  for  describing  the  diffraction  phenomena  serve  the  equations  of 
Maxwell  with  the  appropriate  limiting  conditions  and  the  conditions 
for  radiation/emission. 

The  contemporary  asymptotic  theory  of  diffraction  can  be  defined 
as  the  approximation/approach  to  Maxwell’s  theory,  suitable  in  the 
extreme  case  of  small  wavelengths  and  large  radii  of  curvature  of  the 
diffracting  bodies  (we  eliminate  from  the  examination  of  case  of 
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chisel  edges). 

Page  8. 

The  asymptotic  theory  of  diffraction  introduced  its  inherent 
principles.  Is  established/ installed  the  principle  of  local  field 
into  the  region  of  penumbra  on  the  surface  of  convex  body,  and  also 
the  generalization  of  this  principle,  which  makes  it  possible  to 
apply  it,  also,  to  the  region,  which  adjoins  the  body.  Further  we 
have  the  approximate  form  of  limiting  conditions  on  the  surface  of 
the  well  conductive  body  (Leontovich) ;  Leontovich’s  conditions  are 
impedance  type  conditions  and  escape/ensue  from  the  local  character 
of  field  in  the  surface  layer  of  conductive  body.  These  principles 
allowed  us  to  obtain  explicit  field  expressions  close  to  and  on  the 
very  surface  of  the  well  conducting  convex  body  of  arbitrary  form. 
(Our  formulas,  although  completely  general/common/total,  were 
obtained  by  applying  the  principle  of  local  field  to  the  case  of 
diffraction  on  the  paraboloid  of  revolution;  their  direct 
conclusion/output  from  the  2  integral  equations  established/installed 
by  us  in  chapter  it  was  subsequently  given  by  Kellon  [28].). 

Further  principle  is  the  introduction  of  different  scales  for 
the  horizontal  and  vertical  distances  above  the  conductive  body 
(above  the  surface  of  the  Earth);  this  is  allowed  (Leontovich) .  The 
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latter  takes  the  form  of  the  quantum-mechanical  equation  of 
Schroedinger  (or  the  equation  of  diffusion  with  the  alleged  diffusion 
coefficient),  in  whom  the  time  is  replaced  to  the  horizontal 
coordinate.  Parabolic  equation  makes  it  possible  to  introduce  the 
concept  of  transverse  diffusion  (Malyuzhinets) . 

In  the  first  part  of  this  book  is  constructed  the  asymptotic 
theory  of  diffraction.  The  construction  of  theory  is  realized  in  two 
ways:  first,  on  the  basis  of  strict  solutions  of  equations  of  Maxwell 
(these  solutions  are  converted  to  the  approximate  form,  which 
allows/assumes  physical  interpretation  and  suitable  for  the  numerical 
calculations)  and,  in  the  second  place,  by  the  direct  use/application 
of  principles  indicated  above,  in  particular  parabolic  equation.  In 
certain  cases  are  applied  both  methods;  this  makes  it  possible  to 
compare  them  between  themselves  and  gives  supplementary 
substantiation  to  "asymptotic"  method.  In  the  majority  of  the  cases 
the  field  is  examined  only  near  the  surface  of  the  diffracting  body, 
but  some  formulas  of  latter/last  Chapter  One  part  relate  also  to  the 
large  distances  from  the  surface.  The  necessary  tables  are  given  in 
addition  3. 

In  the  second  part  of  the  book  is  developed  the  theory  of 
radiowave  propagation.  In  first  Chapter  Two  part  is  examined  the 
propagation  in  homogeneous  atmosphere;  the  object/subject  of  these 
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chapters  pertains,  thus,  to  the  diffraction  in  the  true  sense,  and 
they  are  the  straight/direct  continuation  of  the  first  part.  In 
further  chapters  is  examined  the  heterogeneous  (laminar)  atmosphere 
with  the  refractive  index,  which  depends  only  on  the  height/altitude. 

Page  9. 

It  is  here  appropriate  to  note  that  the  case  of  stratified 
atmosphere  gives  another  series/row  of  examples  of  the  new  physical 
concepts,  which  appear  as  a  result  of  applying  the  approximation 
methods.  Such  concepts  are:  the  concept  of  the  given  refractive  index 
(combining  the  effect  of  refraction  with  the  effect  of  the  curvature 
of  the  earth/ground),  the  concept  of  radio  duct  and  finally  the 
concept  of  the  horizons/ levels,  determining  the  range  of  propagation 
in  the  case  of  superrefraction.  Let  us  note  that  with  the  concept  of 
radio  duct  are  connected  two  different  concepts,  namely:  the  seized 
waves  (oscillation)  and  the  waves  (rays/beams) ,  reflected  from  the 
boundaries  of  waveguide.  These  two  concepts  are  supplementary,  in  the 
sense  that  about  the  reflection  in  the  sense  of  geometric 
optic/optics  it  is  possible  to  speak  only  if  excited  not  one,  but 
there  are  many  oscillations. 

The  latter/last  chapter,  dedicated  to  the  phenomenon  of  coastal 
refraction,  is  at  the  same  time  an  example  of  use/appiication  to  the 
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tasks  of  the  diffraction  of  the  integral  equations  of  the  specific 
form  (namely:  equations  with  the  semi- infinite  limits  and  vith  the 
nucleus,  depending  on  the  absolute  value  of  a  difference  in  the 
arguments).  Since  integral  equations  of  this  form  are  encountered 
both  in  the  task  about  the  radiation/emission  of  waveguide  vith  the 
open  end/leafi  and  also  in  other  tasks  of  mathematical  physics,  we 
give  in  addition  1  complete  mathematical  theory  of  such  integral 
equations.  In  addition  2  is  given  the  survey/ coverage  of  properties 
and  uses/appli cat ions  of  Airy’s  functions,  which  play  very  large  role 
in  the  asymptotic  theory  of  diffraction,  and  are  also  given  the 
four-place  tables  of  these  functions. 

The  target  of  this  book  is,  in  the  first  place,  the  presentation 
of  the  general  theory,  developed  in  the  investigations  of  the  author, 
and  only  in  the  second  turn  -  the  analysis  of  numerical  results. 
Nevertheless,  when  the  demonstrative  quantitative  interpretation  of 
general  formulas  is  difficult,  we  give  also  numerical  results  in  the 
form  of  graphs/curves  and  tables. 

For  facilitating  the  ’use  of  the  book  we  premise  to  each  chapter 
and  each  addition  short  abstract;  this  will  allow  the  reader  with  the 
fugitive  survey  of  the  book  to  already  obtain  representation  about 


her  content 
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Since  this  hook  represents  the  collector/collection  of  the 
original  work  of  the  author  in  their  initial  form  (with  the 
insignificant  changes),  in  it  are  unavoidable  the  repetitions.  We 
think,  however,  that  this  book  with  a  sufficient  clarity  shows  and 
logical  the  development  of  the  ideas,  which  compose  the 
object /subject  of  the  works  of  the  author  on  the  theory  of 
diffraction  and  propagation  of  electromagnetic  waves. 

Some  works,  connected  with  this  collector/collection,  are 
written  by  me  with  the  co-authors  whose  names  are  indicated  in  the 
appropriate  chapters.  Me  it  would  like  to  express  here  to  all  my 
co-authors  her  deep  appreciation  for  the  collaboration. 

Page  10. 

The  first  edition  of  this  book  left  in  1965  in  the  English 
language  as  the  first  volume  of  the  series  of  monographs  on  the 
electromagnetic  waves,  the  international  editorial  board  published  on 
the  resolution.  (V.A.Fock,  Electromagnetic  Diffraction  and 
Propagation  Problems,  Pergamon  Press,  London,  1965,  Vol.  I  of  the 
International  Series  of  Monographs  in  Electromagnetic  Waves). 

This  Russian  edition  differs  from  English  in  terms  of  the 
inclusion/ connect ion  of  chapter  16,  dedicated  to  coastal  refraction. 
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and  addition  1,  dedicated  to  the  theory  of  integral  equations.  It 
represents,  thus,  the  expanded  (in  comparison  with  English)  meeting 
of  the  work  of  the  author  on  this  question.  Although  all,  or  almost 
all,  these  works  were  already  printed  in  Russian  in  the  form  of 
various  articles,  to  use  them  is  sometimes  difficult,  since  these 
articles  were  scattered  on  the  different  journals,  and  the  time  of 
their  publishing  covers  the  almost  twenty-year  period  (from  1944 
through  1963).  Meanwhile  how  we  can  judge,  our  results  did  not  cease 
to  be  urgent/actual.  Therefore  we  hope  that  this  edition  of  our  works 
on  the  theory  of  diffraction  and  propagation  of  electromagnetic  waves 
answers  the  actually/really  ripened  necessity. 
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Chapter  13. 


PROPAGATION  OF  DIRECT  WAVE  AROUND  THE  EARTH  WITH  CONSIDERATION  OF 
DIFFRACTION  AND  REFRACTION.  1 . 


FOOTNOTE  l.  Fok,  1948.  ENDFOOTNOTE. 

Is  derived/concluded  approximate  solution  of  the  equations  of 
Maxwell  for  the  vector  of  hertz  of  vertical  dipole,  which  considers 
both  the  diffraction,  and  refraction.  This  solution  befits  under  the 
very  general/common/total  assumptions  about  the  course  of  the 
refractive  index  of  air  in  the  dependence  on  the  height/altitude.  In 
some  practically  important  cases  the  solution  can  be  expressed 
through  the  functions,  introduced  for  uniform  atmosphere.  Here  the 
significant  role  plays  the  concept  of  an  equivalent  radius  of  the 
Earth.  This  concept,  initially  introduced  on  the  basis  of 
considerations  about  the  curvature  of  ray/beam,  can  be  substantiated 
and  generalized.  It  proves  to  be  applicable  also  in  the  shadow  zone 
and  penumbra  where  the  geometric  optic/optics  is  not  applicable. 


Introduction. 
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Under  the  assumption  of  the  uniformity  of  the  earth’s  surface 
the  radiovave  propagation  around  the  Earth  is  caused,  in  essence,  by 
the  following  three  facts:  diffraction  around  the  convexity  of  the 
Earth,  by  refraction  in  the  lower  layers  of  the  atmosphere  and  by 
reflection  from  the  ionosphere.  At  the  small  distances,  the  order  of 
hundred  or  several  hundred  kilometers,  the  reflection  from  the 
ionosphere  role  does  not  play.  However,  at  the  distances  of  the  order 
of  thousand  or  several  thousand  kilometers  the  reflection  from  the 
ionosphere  begins  to  play  the  significant  role,  since  to  the 
straight/direct  (diffractive)  wave,  which  did  not  undergo  reflection, 
begin  to  be  superimposed  the  waves  reflected,  which  can  have 
considerably  greater  intensity  than  direct  wave. 

However,  at  these  large  distances  it  is  possible,  under  the 
known  conditions,  to  isolate  direct  wave  and  to  observe  it 
separately.  Its  study  is  of  great  practical  interest  for  interference 
methods  of  the  determination  of  distances. 

Page  263. 

Therefore  the  construction  of  theory,  which  would  give  amplitude  and 
phase  of  direct  wave  up  to  the  largest  distances,  represents  task 
very  important  for  the  practice. 
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The  theory  of  direct  wave  must  consider  both  the  diffraction  and 
refraction.  Meanwhile  in  view  of  the  complexity  of  task  in  the 
majority  of  theoretical  studies  atmospheric  refraction  either  was  not 
considered  completely  or  it  was  considered  very  roughly,  by  the 
methods  of  geometric  optic/optics.  Very  important  in  this  question 
concept  about  equivalent  radius  of  the  Earth  did  not  obtain  a 
sufficient  theoretical  substantiation.  It  was  introduced  on  the  basis 
of  considerations  about  the  curvature  of  ray/beam,  while  in  the 
region  of  penumbra,  and  those  more  in  the  shadow  zone,  the  concept  of 
ray/beam  generally  becomes  meaningless.  In  connection  with  this  were 
not  explained  those  conditions,  with  executing  of  which  the 
replacement  of  a  radius  of  the  Earth  by  an  equivalent  radius  is 
legal. 

In  the  present  work  we  will  give  approximate  solution  of  the 
Maxwell  equations  for  the  vector  of  hertz,  that  considers  as 
diffraction,  so  refraction.  This  solution  befits  under  the  very 
general/common/total  assumptions  about  the  course  of  tna  refractive 
index  of  air  in  the  dependence  on  the  height/altitude. 

In  some  practically  important  cases  this  solution  can  be 
expressed  through  the  functions,  introduced  by  us  in  our  solution  of 
the  problem  about  the  radiowave  propagation  in  homogeneous 
atmosphere.  These  functions  partly  are  already  tabulated;  when  there 
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are  tables,  the  calculation  of  field  taking  into  account  refraction 
does  not  compose  large  labor/work.  Incidentally  we  will  give  the 
substantiation  of  the  concept  of  an  equivalent  radius  of  the  Earth, 
will  show  that  this  concept  is  applicable  in  the  shadow  zone  and 
penumbra  (where  the  geometric  optic/optics  is  not  applied),  and  let 
us  explain  those  conditions  when  the  use  of  it  is  lawful. 

1.  Differential  equations  and  limiting  conditions  of  task. 

Let  us  designate  through  r,  •&,  o  spherical  coordinates  with  the 

beginning  in  the  center  of  terrestial  globe  and  with  the  polar  axis, 

passing  through  the  radiating  dipole.  Dipole  we  will  assume 

arranged/located  on  the  earth's  surface  and  will  study  the  field ^of 
in  the  air.  We  designate  the  Earth's  radius  by  a.  The  dielectric 

constant  of  air  we  will  consider  function  from  height/altitude 

b=r-a  -above  the  earths  surface . 

t  =  t  ( h ),  h  —  r  —  a.  (l.nu 

Page  264. 

As  in  the  case  of  homogeneous  atmosphere,  field  component  in  the 
air  can  be  expressed  through  the  function  of  hertz  U.  We  have: 
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c  i  n  /  di  \ 

£^-7T-0-(frlF)’ 


W  _  _ ;b  c  gv 

nl  ‘E0E  1 


(1.02) 

(1.03) 

(1.04) 


whereas  remaining  field  component  are  equal  to  zero.  The  dependence 
of  field  component  on  the  time  we  assume  in  the  form  where. 


u»  _ ^  _  2.i 


(1.05) 


Here  X„  -  wavelength  in  the  void  (in  its  our  task  it  is  necessary  to 
differ  from  wavelenath  in  the  air). 


The  value  of  the  dielectric  constant  of  air  on  the  earth’s 
surface  we  will  designate  through  et.  -e(0)  and  let  us  assume. 


—  E(r 


(1.06) 


The  field,  expressed  by  formulas  (i.02)-(1.04) ,  will  satisfy  the 
equations  of  Maxwell,  if  function  U  satisfies  equation 
—  {  -1-  - - - -r  ( S-'.n  0  1  -  klarl'  =  0.  (1.07) 

Let  us  introduce  new  function 

L‘,  =  erV  sinft  V.  (1.08) 

This  function  must  satisfy  equation 


Jl  i 

f  1  8L't\  . 

1  \ 

or  ' 

k  £  or  ) 

•  t  r-  [  db*  ' 

\  4  ' 

4  sin*  i)J 

Field  must  satisfy  on  the  earth’s  surface  Leontovich’s 


£«  =  — 


conditions: 


(1.10) 
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where 


t)  =  es  -  z 


.  4JUI; 


(1.111 


there  is  composite  dielectric  constant  of  soil. 


Page  265. 


Leontovich's  condition  will  be  made,  if  function  U2  satisfies 
condition 


(npu  r  =  a). 


(1.12) 


Let  us  isolate  in  function  U*  the  rapidly  changing  factor,  after 
placing 


.isovir 

U  C;  =  e  Lj 


(1.13) 


where  k  has  value  (1.06)  and  s=a-6  is  the  arc  length  of  terrestial 
globe  from  the  point  where  is  located  dipole,  to  the  point,  above 
which  is  calculated  the  field. 


For  function  U2  is  obtained  the  equation 

i  —  A®  i  —  -  K)  V.  = 

-iLHlh. [J _ 1 _ \  1 

e  dr  r~  [  '  \  4  '4  sin  -i))  5j  * 


(1-14) 


Equation  (1.14)  is  written  so  that  in  its  left  side  stand  the 
dominant  terms,  and  in  the  right  -  correcting,  which,  as  we  will 
show,  can  be  replaced  with  zero. 
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During  the  estimation  of  the  order  of  magnitude  of  derivatives 
we  can  use  the  results,  obtained  for  the  case  of  homogeneous 
atmosphere.  If  we  introduce  the  "high  parameter" 


m  ~  V  2  ’ 


(1.15) 


then  it  will  be 

where  symbol  0(x)  indicates,  as  it  is  accepted,  the  "value  of  order 


x". 


On  the  other  hand,  if  we  exclude  from  the  examination  the 
ionosphere  (in  which  e  it  can  turn  into  zero),  then  the  gradient  of 
logarithm  there  will  be  the  order  of  the  curvature  of  the  earth’s 
surface,  so  that 


(1.17) 


Page  266. 

Hence  it  is  apparent  that  the  individual  members  of  left  side  (1.14) 
will  be  by  order  of  value  not  less  ■£■(/*  whereas  in  the  right  side 
the  members,  who  contain  derivatives,  will  be  on  the  order  of 
However,  concerning  the  term,  which  contains  sinJ-d  in  the 
denominator,  then  under  the  condition  Us. 

ks  >  m  (1.18j 

it  will  be  also  small.  Thus,  throwing/rejecting  the  values  of  order 
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1/m*  in  comparison  with  unity,  we  can  right  side  (1.14)  replace  with 
zero,  after  which  we  will  obtain 


dr- 


21 


(1.19) 


This  is  a  parabolic  equation  of  our  task,  the  reminiscent  of  in  form 
equation  of  Schroedinger  of  quantum  mechanics. 


In  the  obtained  equation  we  can  make  further  simplifications, 
after  using  the  approximate  equality 

1-4  =  2  4-  {l-20> 

After  introducing,  furthermore,  instead  of  the  angle  -d  arc  length 
s=a-»  and  by  considering  s  and  h  as  the  independent  variables,  we  will 
obtain . 

<U1> 

Limiting  condition  for  U2  on  the  earth's  surface  there  will  be  the 
same  as  for  Ulf  namely 

2jj±  =  —iky±Ut  (npu  h  =  0).  (1.22) 

Condition  at  infinity  (h-*«)  can  be  obtained  from  the 
examination  of  the  phase  of  the  function  of  hertz.  Let  us  assume 

U  =  \V\4*  h  Ut  e*  <*-*•’.  (1.23) 

Since  we  examine  the  wave,  which  goes  om  source,  then  phase  Q>  must 
grow/rise  with  an  increase  in  altitude  h.  Hence  we  obtain  the 
condition 

£>«.  (1-24, 

which  must  be  made  at  least  with  sufficiently  large  h. 
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Furthermore,  the  function  of  hertz  U,  and  also  function  U,  must 
remain  final  and  continuous  in  entire  space,  with  exception  of  the 
region,  which  adjoins  the  source. 

Page  267. 

For  the  single-valued  determination  of  the  solution  of  equation 
(1.21)  to  us  it  remains  to  formulate  the  condition,  equivalent  to  the 
volume  which  must  satisfy  function  U2  in  the  region,  close  to  the 
source.  First  of  all,  it  is  obvious  that  in  directly  the  nearness  to 
the  source  is  disrupted  inequality  (1.18)  and  equation  itself  (1.21) 
becomes  wrong.  Therefore  condition  in  the  region,  close  to  the 
source,  must  be  replaced  by  equivalent  condition  in  the  "wave  zone". 
We  can  take,  for  example,  the  region  where  is  applicable  "reflecting 
formula",  and  obtain  the  desired  condition,  after  requiring  so  that 
the  unknown  solution  in  this  region  would  convert/transfer  into  it. 

Reflecting  formula  takes  the  form 

U  =  U-25> 

where  f  -  Fresnel’s  coefficient. 

Since  we  use  the  limiting  conditions  of  Leontovich  (1.10),  then 
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we  thereby  assume  V  >  1-  If,  furthermore,  ve  will  consider  that  h«s, 
i.e.  to  examine  small  inclinations  of  ray/beam  to  the  earth's 

surface,  then  ve  can  asstime. 

h‘‘  n.26» 

K  s  Hi’  1  k\'l ;-s‘ 

Substituting  these  expressions  in  (1.25),  we  come  to  the  conclusion 
that  in  that  region  where  is  applicable  "reflecting  formula",  the 
function 


U.  =  e~'**er  1  sin  i<  L’ 
must  be  reduced  to  the  form 


(i-27) 


V.  = 


1. 1 


I  i  I:  l  V 


L1:  -» 


I  I.-2S: 


Mathematically  this  condition  is  equivalent  to  requirement  so  that 
with  s->0  and  h>0  the  function  U,  would  have  a  special 
feature/peculiarity,  characterized  by  the  formula 


.  h* 

iim  !  —  ~  ” e  j  =  0. 

!-*•  ■  l'  s 


•>lli 


The  more  detailed  substantiation  of  condition  (1.29)  is  given  in 
Chapter  11. 
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Let  us  note  that  instead  of  conditions  (1.29)  or  (1.28)  we  could 
supply  even  the  more  rigorous  condition,  after  requiring  so  that  in 
that  region  where  the  effect  of  the  curvature  of  the  earth's  surface 
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and  heterogeneity  of  the  atmosphere  no  longer  is  manifested  and  vhere 
is  applicable  Weyl-van  der  Pol’s  formula  x,  our  solution  would 
convert/transfer  into  the  decision  of  Weyl-van  der  Pol. 

FOOTNOTE  1 .  The  region  of  the  applicability  of  Weyl-van  der  Pol's 
formula  is  in  detail  investigated  in  work  [223  and  in  Chapters  10  and 
11.  ENDFOOTKOTE. 

2.  Transition  to  are  dimensionless  values  - 


The  differential  equation  derived  by  us  for  function  ls  takes 


the  form. 


§>  -  2 (l~^ Ut  =  0.  (2.C 

ok -  ■  Oi  \  c,  •  a  J 


Let  us  examine  coefficient  with  U2  in  this  equation.  After 
designating  through  to  the  value  of  the  gradient  of  dielectric 
constant  on  the  earth's  surface,  we  can  isolate  in  the  expression  for 
linear  term  and  write  coefficient  with  U2  in  the  form 

**  C-tt  -  *] '  <2  02> 


Let  us  assume  now 


a  •  ae,' 


(2.03) 


Value  (2.03)  is  a  difference  between  the  curvature  of  the  earth’s 


surface  and  the  curvature  of  ray/beam,  and  value  a*  is  conventionally 
designated  as  an  equivalent  radius  of  the  Earth.  After  accepting 
designation  (2.03),  we  can  formula  (2.02)  write  in  the  form 
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where 


f-f.. 


«  (i  -g). 


(2.04; 


(2.05; 


As  can  be  seen  from  (2.05),  value  g  is  the  difference  expressed 
in  the  arbitrary  (deprived  of  dimensionality)  units  between  the 
average/mean  in  the  height/altitude  value  of  the  gradient  of  the 
dielectric  constant  of  air  (undertaken  in  the  section  from  the 
earth’s  surface  to  the  given  height/altitude)  and  its  value  on  the 
earth’s  surface.  In  the  standard  atmosphere  value  g  is  positive, 
whereas  in  the  case  of  temperature  inversion  it  can  become  negative 
and,  only  beginning  from  certain  height/altitude,  it  again  becomes 
positive.  In  terms  of  the  absolute  value  value  g  usually  is  nor  more 
than  two-three  tenths. 

Page  26S. 

With  h-*®  the  theoretical  value  g  exists  -l~a.  while  with  h=0  there 

O 

will  be  g=0.  With  the  standard  atmosphere  value  g  is  changed  very 
slowly,  and  in  the  case  of  inversion  its  change  occurs  considerably 
more  rapid. 

Substituting  expression  (2.04)  in  differential  equation  (2.01), 
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ve  will  obtain 


-g)L\  =  Q. 


(2.06j 


For  the  investigation  of  equation  (2.06)  it  is  convenient  to 
switch  over  from  h  and  s  to  dimensionless  quantities.  For  this 
purpose  let  us  introduce  the  vertical  and  horizontal  scales 


k 


t 


('2.07 ) 


and  let  us  assume 


(2.0?! 


In  order  to  simplify  condition  (1.29),  we  will  switch  over  also  to 

the  new  dimensionless  function  W, ,  after  assuming 

V.  =  ft",.  i2  u?i 

i  a 


Furthermore , 


let  us  assume 


(2-!0j 


In  the  new  designations  differential  equation,  limiting 
condition  and  condition,  which  is  determining  special 
feature/peculiarity,  will  be  written: 
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_t/(i _g)r,  =  0,  (2.1!) 

oir  ox  v  v  ° 

^Ll  _  q\\\  =  0  (npis  it  =  0).  (2  12) 

oy  n  1  " 

limflTj - l=e  ^\  =  0  (y>0).  (2.13) 

x-»l'  \  1  A 

Key:  (1).  with, 

Furthermore,  remains  valid  condition  for  phase  <1>  =  *>  —  arc  U',,  namely: 

f  >0  (npM  y »  1).  (2.1-D 

Key:  (1).  with. 

Page  270. 

The  entering  equation  (2.11)  value  g  was  defined  above  [formula 
(2.05)]  as  function  from  height/altitude  h.  Let  us  designate  through 
h,  certain  height/altitude,  which  characterizes  the  speed  of  a  change 
in  the  gradient  of  the  dielectric  constant  of  air,  for  example  that 
height/altitude  for  elongation/extent  of  which  the  gradient  is 
changed  in  e=2.718...  once.  (For  standard  atmosphere  ho=7400  m;  in 
other  cases  it  is  possible  to  indicate  only  the  order  of  magnitude 
/i0.,  which  to  us  only  is  needed.)  The  value  of  g  we  can  consider  as 
the  function  of  ratio  h/h0: 

£  =  £(-£),  £(0)  =  0,  (2.15) 

considering  that  the  derivative  of  this  function  of  its  argument 
there  will  be  the  order  of  one.  Upon  transfer  to  dimensionless 
quantities  (2,08)  we  had  to  consider  g  as  function  from  y.  So  h^l^y, 
we  will  have 


e  =  g  (M. 


(2.16; 
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where 


(2.17) 


Subsequently  we  will  consider  the  parameter  £  low  value.  In 
order  to  rate/estimate  its  order,  let  us  assume  ho=7400  m  (standard 
atmosphere)  and  let  us  replace  the  equivalent  radius  a*  with  the 
geometric  radius  a.  Then  for  X=1  m,  X=10  m,  X=100  m,  X=1000  m  is 
obtained  respectively  £=0.006,  £=0.027,  /3-0 - 13 r  /3=0.58.  In  the  case 
of  inversion  value  he  will  be  considerably  less,  and  the  parameter  £ 
will  prove  to  be  small  only  for  the  respectively  smaller  wavelengths. 


3.  Solution  of  equations. 


If  we  in  the  equation 

- «' ^  -  y  1 1  “  *<&)]  =  0  <3-0!> 

will  count  £=0,  then,  since  g  (0)=0,  function  g  identically  will 
become  zero,  and  equation  will  be  led  to  that  such  as  was  examined 
and  solved  [together  with  the  limiting  conditions  (2.12)  and  (2.13)3 
in  the  preceding/previous  chapters  of  this  book,  dedicated  to  the 
investigation  of  the  case  of  homogeneous  atmosphere.  It  is  important, 
however,  to  note  that  the  condition  £=0  corresponds  not  to  assumption 
about  the  uniformity  of  the  atmosphere,  but  to  more  general 
assumption  about  the  constancy  of  the  gradient  of  dielectric 
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constant.  In  this  more  general  case  of  formula  are  obtained  the  same 
as  in  the  case  of  homogeneous  atmosphere,  only  in  the  expressions  for 
x,  y  and  q  enters  instead  of  a  radius  of  Earth  a  the  equivalent 
radius  a*. 

Page  271. 


Thus,  the  smallness  of  value  jS  characterizes  the  degree  of  accuracy 
with  which  it  is  possible  to  use  (for  the  finite  values  of  y)  the 
concept  of  an  equivalent  radius. 


Solution  for  /3=0  have  obtained  we  in  the  form  of  the  integral, 
which  contains  composite  Airy's  function.  The  latter  is  that  solution 

of  the  differential  equation 

wr  (/)  =  tw  (/),  (3.02) 

which  has  with  large  negative  t  the  asymptotic  expression 


c(/)  =  e‘ 


_ 1_ 

4  (—  (I  Ae‘ 


(3.03) 


Solution  for  g(/3y)=0  takes  the  form 


W 


w  (I  —  y) 

if'  U)  —  qzi'(l) 


dt. 


(3.04) 


where  duct/contour  C  goes  from  t=i®  to  t=0  and  from  t=0  to 
t  =  oo-ela(0  <  a  <—  J,  covering  all  roots  of  the  denominator  of 
integrand.  (This  duct/contour  it  is  possible  it  goes  without  saying 
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to  replace  with  some  other  equivalent  contour.)  this  solution 
coincides  with  the  fact  which  was  obtained  earlier  in  Chapter  10  and 
in  book  [22] . 

Let  us  attempt  analogously  to  find  the  solution  of  our  equations 
for  the  general  case  /3*=0.  In  this  case  we  will  not  first  assume  )3 
small  and  only  then  for  the  purpose  of  simplification  in  the  finished 
solution  we  will  use  assumption  about  the  smallness  /3. 

Equation  (3.01)  allows/assumes  separation  of  variables  the 
particular  solution  of  equation  (3.01),  which  has  the  form  of  the 
product  of  function  from  x  for  the  function  from  y  and  containing  the 
arbitrary  parameter  t,  will  be  written 

r  j  =  e11'/  (//,  /).  (3.05) 

where  f(y,  t)  satisfies  the  equation 

■  -  [y-t  -ygimi-v.  (3.06) 

From  the  theory  of  differential  equations  it  is  known  that  if  the 
initial  values  (i.e.  the  value  with  y=0)  of  function  f  and  its 
derivative  on  y  are  the  integral  functions  of  parameter  t,  then 
integral  of  equation  (3.06)  will  be  whole  transcendental  function 
from  t. 
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We  will  understand  by  f(y,  t)  that  integral  of  equation  (3.06),  which 
is  whole  transcendental  function  from  t  and  allows/assumes  at  the 
high  values  of  difference  y-t  (or  its  real  part)  the  asymptotic 
representation 


/O'-  0  =  — 


Ce 


>  y  —  <  t  y&  (Pi/) 


—  exp 


¥ 

i  f  Yu  —  t  -i-ug{g>u)du 


.  (3.07) 


Lower  limit  z  in  the  integral,  which  stands  in  the  index,  can  be 
undertaken  arbitrarily.  Factor  C  can  be  function  from  parameter  t. 
Phase  factor  </  <  is  added  so  that  with  g=0  and  r^t  expression 
(3.07)  would  convert/transfer  into  the  asymptotic  expression  for  the 
function 

}{y,  t)  =Cw(l  —  y).  (3.08J 


Expression  (3.07)  is  undertaken  in  accordance  with  requirement 
aa> 

~g£->0,  by  that  superimposed  for  the  phase. 


Understanding  by  f(y,  t)  the  recently  definite  integral  of 
equation  (3.06),  let  us  examine  the  expression 


(3.09) 


where  the  duct/contour  r  .  takes  the  form,  analogous  to  duct/contour 
C  in  integral  (3.04). 
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Let  us  first  of  all  note  that  the  integrand  in  it  is  determined 
by  the  conditions  superimposed  above  unambiguously,  since  remaining 
indefinite  in  (3.07)  factor  C  in  it  it  is  shortened. 

Further,  integrand  in  (3.09)  presents  meromorphic  complex 
variable  function  t:  unique  singular  points  in  it  are  the  roots  of 
denominator. 

The  investigation  of  the  roots  of  denominator  in  (3.09)  is 
difficult  to  conduct  witn  complete  strictness.  For  this  investigation 
it  is  necessary  to  know  the  behavior  of  function  g(j3y)  at  the 
composite  values  of  y  near  the  half -line  arc  y~t/Z.  However,  on  the 
basis  of  some  not  completely  strict  considerations  which  we  do  not 
here  give,  it  is  possible  to  expect  that  if  function  g(/3y)  in  complex 
domain  indicated  will  remain  small  (for  example,  Igl<-g-)»  then 
roots  will  be  arranged/located  in  the  manner  that  in  the  case  of  g=0, 
i.e.,  in  the  first  quadrant  of  plane  t  near  the  half-line  ars  t=ir/3. 
In  any  case  this  will  be  so  at  the  low  values  of  the  parameter  /3. 

Page  273. 


We  should  also  know  the  behavior  of  function  f(y,  t)  at  the 
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positive  values  of  t-y  (and  also  in  certain  sector  of  plane  t,  which 
switches  on  positive  real  axis).  The  unknown  asymptotic  eqpression 
will  be  obtained  by  the  analytical  continuation  of  expression  (3.07) 
through  the  third  and  fourth  quadrants  of  plane  t,  since  in  the  first 
are  arranged/located  roots  f(y,  t).  It  will  take  the  form 

"  T 

Hu.  i)—  1  -  - exp  [|  /  —  u  —  ug  (f5u)  du  .  (3. 1 0) 

,  /  —  u  —  ygtfu)  Li 

If  we  assume  here  g=0  and  to  take  r=t,  then  this  expression  will 
be  led  as  (3.07),  to  the  asymptotic  expression  for  function  (3.08). 

Knowing  the  location  of  roots  and  the  behavior  of  integrand  on 
both  sides  of  the  region  where  these  roots  are  arranged/located, 
already  it  is  possible  to  conduct  in  integral  (3.09)  duct/contour  f. 
so  that  it  would  cover  all  roots  of  denominator  and  would  become  by 
both  branches  infinite.  On  the  initial  branch  of  duct/contour  (going 
from  infinity)  there  will  be  correctly  asymptotic  expression  (3.07), 
while  on  the  final  branch  (exiting  to  infinity)  -  expression  (3.10). 
In  this  case  the  integral,  undertaken  on  this  duct/contour,  will  be 
that  converging. 

The  preceding/previous  reasonings  had  as  a  goal  to  show  that 
expression  (3.09)  for  function  makes  thinned  mathematical  sense. 

Let  us  show  now  that  it  satisfies  all  conditions  presented. 


i 

i 

i 
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It  is  first  of  all  clear  that  it  satisfies  differential  equation 

(3.01)  because  it  is  satisfied  by  the  integrand.  Further  it  satisfies 
limiting  condition  (2.12): 

/JIT  *' 

+  qWi  =  0  (ripii  y  =  0).  (3. 1 1)  • 

Key:  (1).  with. 

Actually/really,  producing  in  (3.09)  the  differentiation  under 
the  integral  sign  and  assuming/setting  then  y=0,  we  will  see,  that 
the  numerator  of  fraction  will  be  shortened  with  the  denominator  and 
integrand  will  be  holomorphic,  in  consequence  of  which  the  integral 
will  become  zero.  Then,  integral  will  be  convergent  and,  therefore, 
final  at  all  positive  values  of  x  and  y.  It  is  not  difficult  to  also 
check  that  it  will  satisfy  condition  for  phase  (^>°)‘ 

To  us  it  remains  to  check,  does  have  expression  (3.09)  required 
by  condition  (2.13)  special  feature/peculiarity  near  x=0,  or  which  is 
equivalent  to  this,  to  be  convinced  of  the  fact  that  at  the  small 
distances  from  the  source  it  gives  the  formulas  of  Weyl-van  der  Pol 
or  reflecting  formula. 
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With  the  aid  of  asymptotic  expressions  (3.07)  and  (3.10)  for 

f(y,  t)  it  is  possible  to  show  that  if  x  and  y  are  small,  and  ratio 

y/x  is  great,  then  the  main  section  of  integration  will  lie/rest  at 

the  region  of  the  high  negative  values  t.  (Initial  contour  r  can  be 

deformed  so  that  it  would  pass  through  this  region.).  Using 
expression  (3.07),  we  will  obtain  for  large  negative  t: 


Hence 


and 


/«».  t)  i*  —t 

i<Q.t)-\  u-t~yg(py)  eXp 


/  j  J  u  —  t  -  ugdu 


.  (3.12) 


-7-  i  ‘ 1  y—t  -ryg($y) 


(3.13) 

(3.14) 


But,  when  y  is  small,  terms  yg(/3y)  are  small  in  comparison  with  y, 
and  we  can  instead  of  (3.12)  write 


fjt  H  _  \Y 


/  (0.  I ) 


V 


y  —  t 


exp 


(3.15) 


i  f  ]'u  —  tdu 

L  o 

Let  us  note  now  that  the  same  asymptotic  expressions  will  be  obtained 
for  this  region,  if  we  instead  of  f(y,  t)  substitute 

f  {y.  f)  =w{t  —  y).  (3.16) 


But  after  this  substitution  integral  (3. OS)  will  become  (3.04),  and 
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this  latter  gives  for  small  x,  y  the  formula  of  Weyl  -  van  der  Polya, 
reflecting  formula  and  limiting  condition  (2.13). 

We  can  check  this  and  it  is  direct.  By  introducing  the 
variable/alternating  of  integration  p  =  V—  t  and  disregarding 

values  y  and  ys  in  comparison  with  p,  we  will  obtain 

Hy.  -r*)  _e>„  (3.17) 

i  (0.  -p1) 

and 

(rf +«)«--*’+*•  <318> 

The  substitution  of  these  values  into  integral  (3.05)  gives 

*'»  =  e"  I e_fxp"'i'P  T=W'  (3‘I9} 

r  * 

where  duct/contour  f.  •  it  intersects  positive  real  axis  in  plane  p 
from  below  upward  (near  point  p=y/2x). 

Page  275. 

If  we  will  calculate  integral  (3.19)  without  the  neglects,  then 
will  be  obtained  the  formula  of  Weyl  -  van  der  Polya.  If  we  will 
calculate  it  according  to  the  method  of  steady  state,  obtain 
reflecting  formula.  If  we  disregard/neglect  value^qjin  comparison 
with  y/2x  will  be  obtained  the  expression,  which  is  turned  into  zero 
left  sides  (2.13)  even  to  the  transition  to  the  limit. 

Thereby  it  is  proved  that  expression  (3. IS)  for  V1  presents  the 
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unknown  solution  of  our  problem. 


4.  Investigation  of  the  solution  for  the  zone  of  straight/direcl 
visibility. 


Instead  of  function  Wj  to  more  conveniently  examine  another 
function,  which  differs  from  W,  in  terms  of  factor  1  *•  Me  will  assume 


Vu,  y.  <j)  =  e‘  4  ]/^  l  e'*‘  7^~T  du  (4‘0l) 


Recollecting  connection/ communication  between  functions  U,  Ulr 
U5,  given  by  formulas  (1.08),  (1.13),  (2.09),  and  disregarding 
the  difference  between  r  and  a  ana  between  e  and  £<»  when  these 
values  enter  as  the  factor  with  U,  we  can  write 

u  = - V  (x.  if.  q).  (4-02) 

"I  ^  cs  sir. 

where  s,  as  before  is  the  horizontal  distance,  counted  on  the  arc 
along  the  earth’s  surface,  and  x,  y,  q  are  connected  with  s,  h , 
with  the  relationships/ratios 

*  ..-ft  /.im, 


s  ft  . far  ■«  •  e0 

y=  ?=i1  —  1  t* 


(4.03) 


moreover 


2a*1  .  3 

*»“ V  *>=1 


If  s  is  small  in  comparison  with  a  radius  of  terrestial  globe,  then 
instead  of  sin  s/a  it  is  possible  to  write  simply  s/a  (as  usually  and 
they  write). 
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However,  since  our  formulas  remain  valid  up  to  the  very  large 
distances  where  the  difference  between  sine  and  another  becomes 
noticeable,  we  leave  in  (4.02)  under  root  of  sin  s/a. 


Function  V(x,  y,  q)  can  be  named/called  attenuation  factor;  when 
is  possible  to  count  g=0  and  to  use  the  concept  of  an  equivalent 
radius,  formula  (4.01)  for  V  converts/transfers  into  the  formula 


V  (,v. 


Si'  It - L-l 


u’  c,  —  c- 


a;. 


(4.05) 


that  derived  by  us  earlier  for  the  case  of  homogeneous  atmosphere 
[formula  (6.02)  of  chapter  103. 


The  function  (4.05)  it  is  in  detail  investigated  in  chapter  10 
and  in  our  work  [22]  and  it  is  partly  tabulated  (for  q=0). 
Investigation  will  flow/occur/last  in  many  respects  in  parallel  to 
analogous  investigation  in  chapter  10. 


In  the  present  paragraph  we  will  examine  the  zone  of 
straight/direct  visibility. 


In  the  zone  of  straight/direct  visibility,  it  is  not  too  close 
to  the  horizon,  comes  into  force  geometric  optic/optics.  If  we  will 
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use  expression  (3.12)  and  let  us  introduce  the  variable/alternatinc' 

of  integration  p,  we  will  obtain  for  the  V  integral  of  the  form 

(4.06) 


4/— 


i  =  e 


=L_  ,  •"  v  _ 

'r.  »  xJe  I  y~r*  —  yg(teip  —  u< 


where  for  the  brevity  is  placed  1 


O,  =  —  Xfr  -  I  |  u  -  pr  —  ug($u)du.  (4.07) 

i 

FOOTNOTE  1 .  Phase  u  should  not  be  mixed  with  the  angular  frequency, 
designated  in  paragraph  1  by  the  same  letter.  ENDFOOTNOTE. 


Calculating  integral  according  to  the  method  of  steady  states, 


we  find  from  equation  ~  -  v  or 

f 


du 


2  .*  1  u  —  fr—ug  (?*! 
t- 


(4-08) 


the  extremum  of  phase  and  after  some  linings/calculations  come  to  the 
expression 

(4.09) 


r 


=  e**  — ~~  1/  2x  . 
p-4  r  oy 
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In  this  formula  under  p  is  understood  the  function  from  x  and  y, 
determined  from  equality  (4.08).  With  g=0,  and  also  with  small  x  and 


y 


(4.10) 


and  the  expression  -under  the  sign  of  square  root  in  (4.09)  will 
become  unity.  Formula  (4.09)  is  valid  when  value  p  is  great  and 
positive. 
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Our  formulas  allow/assume  simple  interpretation  from  the  point 
of  viev  of  geometric  optic/optics.  Actually/ really ,  the  complete 
phase 

<1»  =  ks  —  *t>  (4.1 1) 


of  function  U  presents  the  solution  of  equation  for  the  eikonal 

( i  ^  JLy  (i± i:  -  (*£.  f  =  A-  ( l  - ± Y-L. ,  (4  12.i 

which  after  neglect  of  lov  values  reduces  to  the  following  equation 
for  ui 

( -^r-  j~  —  2k  -  h\\  —  gi.  <4.23i 


Here  to  the  right  stands  value  (2.04).  After  transition  to  to  the 
variable/alterna~inc  x,  y  we  obtain  from  (4.13) 

<4.14) 

Relationship/ ratio  (4.06)  is  an  equation  of  the  trajectory  of 
the  ray /beam,  passing  through  the  origin  of  coordinates,  and  value  p 
-  the  parameter  of  this  trajectory.  The  geometric  value  of  parameter 
p  exists 

/>  =  1  —cosy.  (4-15) 


where  7  -  angle  between  the  ray/beam  and  the  vertical  line  near  the 
source . 


Complete  phase  a*  is  optical  path  length  of  ray/beam,  counted 


from  the  source  to  ooint  x,  v.  Value  2t>/p-ic  is  eoual  to 
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where  f  -  Fresnel's  coefficient. 

Thus,  when  is  applicable  geometric  optic/optics,  our  formulas 
convert/transfer  into  the  formulas  of  geometric  optic/optics. 

Page  276. 

Formula  (4.09)  is  applicable  for  final  x  and  y  when  parameter  p 
is  positive  and  great.  But  if  x  and  y  are  small,  is  necessary  another 
condition 

£-■£»'•  (4,7) 

With  the  nonfulfillment  of  condition  (4.17),  in  the  case  small  x 
and  y  and  large  p,  expression  (4.06)  remains  valid  but  integral  must 
be  calculated  otherwise,  namely:  it  is  necessary  to  replace  u  by 
-xp*+yp  and  the  root  of  the  fourth  degree  to  one,  after  which  the 
integral  will  be  reduced  to  form  (3.19)  (with  factor  ] 'x)  and  it  will 
give  the  formula  of  Weyl  -  van  der  Polya. 

Let  us  note  that  if  x  and  ]  'y  are  great,  and  parameter  p  is  low 
in  comparison  with  these  values,  then  the  equation  of  trajectory 
(4.08)  can  be  approximately  solved  relative  to  p.  We  will  have 
approximately 
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P 


-  ,*  \"u-ug{ t*, 


(4.18) 


Under  the  same  conditions  it  is  obtained 


u)  =  <o0Cv)  +-§-  P3-  (4‘19) 

where 

0)0(£/)=  J]  u~ug$u)  (4-20) 

C 

but  symbol  p  it  is  necessary  to  understand  as  the  abbreviation  for 
value  (4.18). 


Equality  p=0  gives  the  geometric  boundary  of  shadow.  If  right 
side  (4.16)  becomes  negative,  then  equation  (4.06)  does  not  have  real 
solution  relative  to  p,  whereas  function  (4.18)  [and  also  (4.10)] 
retains  sense,  also,  in  this  case.  This  apparent  disagreement  is 
explained  by  the  fret  that  right  side  (4.08)  is  not  analytic  function 
from  p  near  p=0. 

Expressions  (4.18)  and  (4.19)  will  be  encountered  to  us  in  the 
region  of  the  penumbra  where  the  geometric  optic/optics  is  already 
not  applied. 


5.  Investigation  of  solution  for  the  region  of  penumbra  (final  x  and 
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where  hj  -  height  scale,  and  t#  '■  *-ue  of  the  second  derivative  of 
on  the  height/altitude  on  the  earth's  surface. 


With  small  /3C  and  finite  y  and  t  as  approximate  solution  of 
equation  (3.06)  it  is  possible  to  take  the  function 

f  (it,  t)  =  (/  —  y)  — 

-  -IT-  [(%  -r  2/)  z £»(/  —  y)~  (3i f  -f-  4 yt  -f  8/*)  w'  (/  —  y)] .  (5.06i 
Substituting  this  expression  in  (5.03),  we  will  obtain  for  the 
unknown  root  the  approximation 


3  -  4/*  <7 


f. 


(5.07) 


where  i°„  -  root  of  the  equation 

(5-08) 

in  detail  investigated  in  chapter  10  and  in  work  [22] . 
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For  function  D(t)  is  obtained  the  expression 

D  (/)  =  (/  —  q')  ( 1  -  -5-  Po» )  +  -f  M-  (5-09t 

The  entering  formula  (5.01)  high-altitude  factors 

<5101 

can  be  obtained  by  the  numerical  integration  of  the  differential 
equation 

~^~[y—!n~ys(Py))ln  =  o 


(5.11) 
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under  the  initial  conditions 

fr. (0)  =  1  H  f„(0)=-q.  (5.12) 

Thus  far  y  is  certain  (whereas  x  can  be  great),  obtained  in  this 
method  values  V(x,  y,  q)  they  will,  with  small  ones  0,  differ  little 
from  the  values  for  /3=0.  More  or  less  essential  difference  can  be 
revealed  only  in  exponential  factors  e,x'«.  of  those  giving  fading  and 
additional  phase.  Therefore  correction  it  suffices  to  introduce  into 
these  exponential  factors. 

But  if  for  special  accuracy  it  is  not  required,  it  is  possible 
to  disregard  this  correction  and  it  is  simple  to  consider  rhat  in  the 
case  in  question  the  expression  for  V(x,  y,  q)  coincides  with  that 
derived  for  the  case  of  homogeneous  atmosphere  (with  condition  of  the 
replacement  of  a  geometric  radius  of  the  Earth  by  an  equivalent 
radius).  Then  it  is  possible  to  use  all  formulas  and  tables,  obtained 
for  this  case. 

6.  Investigation  of  solution  for  the  region  of  penumbra  (large  x  and 

y). 

There  is  greatest  practical  interest  in  the  case,  whereas  when  x 

and  y  is  very  great,  value 

n  =  -L  l' 

2  J  V u  -r  ug  (pm 


X 


(6.01) 
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With  small  and  final  y  and  t 

X  —  I  —  u  —  t="  (3jr  -  4/1/  -  8/*j.  (6.05) 

JO 

Then  the  function 

(6-06) 

will  present  the  solution  of  equation  (3.06)  with  error  in  the  order 

/3:  with  final  y  and  t  and  order  £Vl  with  large  y  and  final  t.  With 

the  aid  of  expression  (6.05)  it  is  not  difficult  to  check  that  in 

expansion  (6.06)  according  to  degrees  />B  the  members  of  order  to  /2„ 

inclusively  coincide  with  (5.06).  However,  expression  (6.06)  is 

correct  even  when  (with  large  y)  expansion  (5.06)  is  not  applicable. 

If  value  X  is  very  great  and  negative  (which  will  be  with  large  y) , 

then  expression  (6.06)  is  led  to  the  following: 

t  f _ ‘ 

U J/-?>  =4 - - - exp  i  i  |  u  — ;  —  ug (pu) du  .  iG.OT, 

\ry—  <-ryg(fc>)  t 

The  latter  coincides  with  (3.07),  if  we  there  assume  C=1  and  to 
take  as  r  the  root  of  equation  (6.04).  Thus,  by  means  of  formula 
(6.06)  we  ascertained  that  one  and  the  same  solution  of  equation 
(3.06)  has  with  final  y  expression  (5.06).  And  with  large  y  - 
expression  (6.07). 
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We  can  now  use  during  the  calculation  of  the  integral 

:.n  r~ 

*  -r  *»  / 


v  1  -T  1  /  x  \  ixt  I  (w. 

Ux.y.q)  =  e  )  —  e  ™ - - 

y  \dt i~q!) 


d! 


(6.08) 


both  expressions  (5.06)  and  (6.07)  simultaneously,  namely  substitute 
in  the  numerator  expression  (6.07)  and  in  the  denominator  - 
expression  (5.06).  In  this  case  it  is  possible  both  expressions  to 
somewhat  simplify.  Throving/re jecting  small  corrections,  we  will 
write  instead  of  (5.06)  simply 

f  <«-  t)  =  c  (/  —  y).  (6  09) 

However,  in  the  formula  (6.07)  in  the  factor  before  the  exponential 
function  let  us  disregard/neglect  value  t  in  comparison  with  y,  and 
irylex  is  replaced  by  the  approximation 

V 

\  1  u  —  t  —  ug  (pu)  du  = 


-  \  I  u  -  ug (P«| du  — x_ /  l  ■— — ul_. 

.  *  |  U  ~ 


(6.10i 


Using  designations  (4.18)  and  (4.20),  we  can  write 

/  («/.  0=]  2  ^  e‘  T  (I+pl .  (6.1U 

In  the  final  analysis  we  replace  function  f(y,  t)  in  the 
denominator  by  Airy’s  function,  while  in  the  numerator  -  by 
exponential  function. 


Substituting  (6.09)  and  (6.11)  into  integral  (6.08),  we  will 


obtain 
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»•  . . ,612) 

r 

The  remaining  integral  can  be  considered  known  function.  In 
chapter  10  [formula  (6.20)  it  is  designated  through 

i  f  —4  pi  _ dt 


...  .  1  |  —ipi  Ul 


(6.13) 


and  it  is  in  detail  investigated.  For  cases  of  p=0  and  <7  =  1  i 
integral  V,  is  tabulated  * . 


FOOTNOTE  1 .  Tables  for  q=0  are  given  in  addition  3.  ENDFOOTNOTE. 


formula  (6.12)  gives  attenuation  factor  for  the  region,  close  to 
the  horizon. 

Page  283. 

It  is  interesting  to  compare  this  formula  with  formula  (4.08),  valid 
in  the  region  of  the  applicability  of  geometric  optic/optics.  Using 
(4.19),  let  us  write  expression  (4.09)  in  the  form 

(6.14) 

But  in  chapter  10  it  is  shown  that  function  (6.13)  has  with  large 
positive  p  the  asymptotic  expression 

VA-p-  <?)  =  -^r7-,eTr*  (6-151 

[formula  (6.24)  of  chapter  10 3.  Thus,  our  formula  (6.12) 
converts/transf ers  in  the  zone  of  straight/direct  visibility  into  the 
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formula  of  geometric  optic/optics. 


With  negative  p  the  expression  for  V j ( -p ,  q)  can  be  written  in 
the  form 


'1  (-P-  =  12 1  n  S  „nI^)?i77 


n.s.1 


(6.16) 


Wnen/pjis  great  (p<0),  this  series/rov  is  reduced  to  the  first 
term  which  gives  the  attenuation  of  wave  in  the  shadow  zone  according 
to  the  exponential  law. 


Function  \\i-p,  q)  was  for  the  first  time  introduced  in  our 
works,  dedicated  to  diffraction  from  the  body  of  the  arbitrary  form 
(see  chapters  1,  2  and  5).  In  these  works  was  established/ installed 
the  principle  of  local  field  in  the  region  of  penumbra  anc  it  was 
shown  that  in  this  region  the  field  is  expressed  as  function  Vj(-p, 
q),  which  has  universal  character. 

The  comparison  of  formulas  (6.12)  and  (6.14)  makes  it  possible 
to  in  a  certain  sense  say  that  the  wave  reaches  the  horizon/level 
with  the  amplitude  and  the  phase,  which  corresponds  to  the  laws  of 
geometric  optic/optics,  and  at  the  horizon/level  it  undergoes 
diffraction  according  to  the  lav  of  local  field  in  the  region  of 
penumbra . 
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This  picture  is  in  complete  agreement  with  L.  I.  Mandelstam's 
ideas  about  the  fact  that  during  the  propagation  of  electromagnetic 
waves  along  the  earth’s  surface  the  properties  of  soil  are  essential 
not  along  the  entire  trajectory  of  ray/beam,  but  only  in  that  region 
where  is  arranged/located  the  transmitter  located  on  the  earth/ground 
or  receiver  ("takeoff"  or  "landing"  area/site). 

If  we  accept  this  picture,  then  solution  obtained  in  this 
paragraph  can  be  applied  also  to  that'  case  when  the  properties  of  the 
earth’s  surface  in  the  different  sections  are  not  identical,  when  in 
function  V.(-p,  q)  the  composite  parameter  q  corresponds  to  the 
properties  of  soil  in  that  region  where  the  ray/beam  concerns  the 
earth/ground. 
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Chapter  14. 

Theory  of  radiowave  propagation  in  an  unhomrgeneous  atmosphere  for 
the  elevated  source  1 .  ENDFOOTNOTE. 

FOOTNOTE  1 .  Foch,  1950.  ENDFOOTNOTE. 

The  theory  developed  in  the  preceding/previous  chapters  of  radiovave 
propagation  is  generalized  to  the  case  of  the  elevated  dipole  in  the 
heterogeneous  (laminar)  atmosphere,  moreover  is  considered  as  normal 
refraction,  so  also  super ref raction. 

In  the  beginning  of  chapter  are  extracted  fundamental  equations 
and  limiting  conditions  of  task.  Then  is  examined  the  approximate 
form  of  equations  (parabolic  equation  of  Leontovich)  with  the 
appropriate  limiting  conditions  and  the  conditions,  which  are 
determining  special  feature/peculiarity.  Further  is  drawn  analogy 
between  the  formulated  task  and  the  unsteady  task  of  quantum 
mechanics.  After  transition  to  the  dimensionless  quantities  are 
studied  the  properties  of  the  particular  solutions  of  differential 
equation  from  which  is  constructed  the  general  solution  in  the  form 
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of  contour  integral  and  in  the  font  of  series/ row.  General  theory  is 
applied  to  the  case  of  superrefraction  and  is  examined  a  schematic 
example.  In  conclusion  ere  derived  the  approximation  formulas  for 
determining  the  attenuation  factors  and  nigh-altitude  factors.  These 
formulas  are  analogous  to  the  semi -classical  formulas  of  quantum 
mechanics. 

Introduction. 

The  theory  of  radiovave  propagation  in  the  atmosphere  vith  the 
dielectric  constant,  which  depends  on  height/altitude,  have  developed 
we  in  chapter  13  for  the  case,  when  source  is  the  vertical  electric 
dipole,  arranged/iocatec  on  the  earth’s  surface.  On  the  other  hand, 
the  case  of  the  elevated  source  (horizontal  and  vertical,  electrical 
and  magnetic  dipoles)  was  examined  by  us  in  chapter  12  under  the 
assumption  of  homogeneous  atmosphere.  At  present  to  chapter  is 
studied  the  "combined"  case:  the  elevated  cipoie  and  unhomogeneous 
atmosphere . 

The  formulas,  derived  in  chapter  13  for  the  general  case  of  the 
arbitrary  course  of  refractive  index,  were  there  in  more  detail 
developed  under  the  assumption  of  the  normal  refraction  when 
radiowa*  e  propagation  has  qualitatively  the  same  character  as  in 
homogeneous  atmosphere. 
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The  case  of  the  superrefraction  when  the  lower  layer  of  the 
atmosphere  acquires  the  character  of  waveguide,  it  is  of  independent 
interest  and  deserves  special  examination.  In  present  chapter  we 
examine  this  case  in  more  detail.  For  its  tualitative  characteristic 
of  the  highly  useful  proves  to  be  the  analogy  with  the  unsteady  task 
of  quantum  mechanics  about  spreading  wave  packet  in  the  assigned 
force  field:  this  analogy,  apparently,  remained,  until  nov,  of  that 
not  noticed. 

1.  Fundamental  equations  and  limiting  conditions. 

Let  us  designate  through  r,  f>,  ©  spherical  coordinates  with  the 
beginning  in  the  center  of  terrestial  globe  and  with  the  polar  axis, 
passing  through  the  radiating  dipole.  A  radius  of  the  Earth  let  us 
designate  through  a.  Dipole  we  will  assume  finding  on  height/altitude 
h'=b-a  above  the  earth's  surface,  so  that  its  coordinates  will  be 
r=b,  -6=0.  The  dielectric  constant  of  air  t  we  will  consider 
function  from  height/altitude  h«r-a  above  the  earth’s  surface. 

Field  in  the  air  can  be  expressed  according  to  the  known 
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ormulas  through  Debye's  potentials  u,  v. 


We  have: 


_  dv_  _  | 


Er  =  - 

£.  =  _ 

1  <■>'  (era) 

10)  dv 

er  dr  dd 

c  sin  d  dq  ’ 

£_  =r  - 

1  d 

r-  (tru)  iio  dv 

er  sin  0 

dr  o<{  ~  ub 

Hr  = 

_ L  A*r. 

r  J  ," 

//*  = 

Uu  c  d:, 

__  1  d-  (r;>) 

1 1  v  — 

c  sm  d  0<i 

'  r  dr  d{)  ’ 

U; 

-n 

II 

lui  ou 

C  Oil 

1  d-  (/;.) 
r  sin  o  or  d(j 

(1.01) 


(1.02) 


The  same  expressions  are  applicable  for  the  field  lower  than  earth's 
surface,  if  we  by  s  understand  the  composite  dielectric  constant  of 
soil.  Dependence  on  the  time  is  assumed  here  in  the  form  of  factor 


—lot 
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Symbol  A*  indicates  Laplace’s  operator  on  the  sphere: 

A*“  =  lmTdF  +  l^Td^'  (l*03) 

The  equations  of  Maxwell  will  be  satisfied,  if  functions  u  and  v 

satisfy  the  equations 


and 


1  6  !  I  d(Enn\ 
r  dr  \  e  dr  ) 


A*u 

r2 


i  (,|J  n 

+  —  ««  =  0 


(1.04) 


r 


&■  Irv)  |  A*t>  , 
dr 2  '  r*  ' 


c* 


eu  =  0. 


(1.05) 
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The  continuity  of  the  tangential  components  of  field  will  be 


provided,  if  there  will  be  also  continuous  values 


1  o(tru)  dirt') 

ErU *  1  6r~  *  ~dr' 


(1.06) 


Hence,  via  known  reasonings,  is  obtained  the  approximate  form  of 
limiting  conditions  (Leontovich ' s  condition).  If  we  will  assume 
k=u/c,  let  us  designate  the  composite  dielectric  constant  of  soil 
through  t).  and  designation  e  let  us  retain  for  the  dielectric 
constant  of  air,  then  we  will  have 


—  —  ik~:(eru)  (npH  r  =  a) 


Key:  (1).  with. 


=  —  ik  1  ‘ti  (rv)  (npn  r  =  a). 


Key:  (1).  with. 


In  further  field  for  which  uj^O,  v=0,  we  will  call  vertically 
polarized,  and  the  field  for  which  u=0,  v^O,  horizontally  polarized. 
The  field  of  vertical  electric  dipole  remains,  in  this  sense, 
vertically  polarized  in  entire  space.  The  field  of  vertical  magnetic 
dipole  (horizontal  framework)  possesses  everywhere  horizontal 
polarization.  However,  horizontal  electric  dipole  excites  the  fields 
both  form:  as  it  is  horizontal,  so  also  vertically  polarized.  In  the 
case  of  homogeneous  atmosphere  the  vertically  polarized  field 
decreases  with  the  increase  of  distance  more  slowly  than  horizontally 
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polarized  (see  Chapter  12).  Therefore  field  from  the  horizontal 
electric  dipole  at  the  small  distances  from  the  source  will  possess 
predominantly  horizontal  polarization,  and  at  large  distances  (in  the 
region  far  beyond  the  horizon/level)  polarization  will  be 
predominantly  vertical. 

Page  287. 

The  vertically  polarized  field  can  be  expressed  through  function 
U  (function  of  hertz  of  vertical  electric  dipole),  which  possesses 
the  following  properties:  function  U  satisfies  tne  same  differential 
equation  (1.04)  and  the  same  limiting  conditions  (1.07)that  also  u, 
and  has  near  the  source  a  special  feature/peculiarity  of  the  form 

JkR 

V  =  •— — r  U*,  (\xm 

where  U*  remains  final.  Values  R  and  k  have  values 

R  =  ]  rs  —  b*  —  2rfccos#,  /:  =  ■—.  (1.10) 

Analogously,  horizontally  polarized  field  can  be  expressed 
through  function  W  (function  of  hertz  of  vertical  magnetic  dipole), 
which  satisfies  the  same  differential  equation  (1.05)  and  the  same 
limiting  conditions  (1.08)that  also  v,  and  has  near  the  source  a 
special  feature/peculiarity  of  the  form 

«"  =  -V  r*.  (1.11) 

where  W*  remains  final. 
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Fields  from  the  vertical  and  horizontal  electrical  and  magnetic 
dipoles  with  moment/torque  M  are  expressed  as  the  specific  above 
functions  U  and  W. 


For  the  vertical  electric  dipole  we  must  assume 

(1.12) 


Af  A 

u  =  —  l'.  v  —  0. 


For  the  vertical  magnetic  dipole  (horizontal  framework)  we  have 

u  =  0.  (1.131 

For  the  horizontal  electric  dipole,  directed  along  the  axis  x, 
entering  formulas  (1.01)  and  (1.02)  functions  u  and  v  are  determined 
from  the  equations 

A*u  =  -A1lr(4--irjCOS(r' 


id 


A*u  =  —  ikM  sin  <f. 


(1.14) 


where  A*  -  operator  of  Laplace  on  sphere  (1.03). 
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Finally,  for  the  horizontal  magnetic  dipole,  directed  along  the 
axis  x,  we  have 


and 


&*u  =  —  ikM  sin  <p 


.*  > ,  c<  /aw  ,  u~\ 

if  =  %U'T)C0S,f" 


( 1 - 1 3» 


Thus,  into  all  four  cases  the  study  of  field  is  reduced  to  the  study 
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of  functions  U  and  W. 

2.  Approximate  form  of  equations. 

Converting/transferring  to  the  approximate  form  of  equations, 
let  us  designate  through  et  the  value  of  the  dielectric  constant  of 
air  near  the  source  (practically  it  is  possible  to  assume  £j =  1)  and 
let  us  assume 

s  =  ad,  (2.01) 

so  that  s  -  horizontal  distance  between  source  and  observation  point, 
counted  on  the  arc. 

Let  us  introduce  instead  of  U  and  W  slowly  changing  functions  U, 
and  W,,  after  assuming 


U  =  — 1 U: 

f.r  V  sin  \i 

(2.02) 

Jts 

r  =  ,e  -  R7S. 

r  J  sin  it 

(2.03) 

As  shown  in  chapter  13,  after  neglect  of  low  values  the  equation 
for  Uj  takes  the  form 

=  (2.04) 

on-  os  \  oj 

As  the  independent  variables  are  here  accepted  instead  of  r  and  0 
values  h  (height)  and  s  (horizontal  distance).  Equation  for  W3  in  our 
approximation/approach  will  take  the  same  form,  namely 
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Equations  (2.04)  and  (2.05)  we  will  call  the  parabolic  equations  of 
Leontov  i  c'n . 

During  the  composition  of  limiting  conditions  on  the  earth's 
surface  (h=0)  we  can  disregard/neglect  difference  between  the 
dielectric  constant  of  air  and  unity. 

On  the  other  hand,  we  can  somewhat  refine  these  conditions, 
using  our  results,  obtained  according  to  the  method  of  the  addition 
of  the  series/rows  (see  chapters  12  and  6). 

Page  28S. 

This  refinement  is  reduced  to  the  fact  that  in  (1.07)  we  replace 
and  <|t  1,  and  in  (1.08)  we  replace  n  on  tj  —  1.  As  a  result  we  will 
obtain 


d'J. 

dh 

r  lk  U. 

V 11-1  * 

LO 

(npii  h  =  0) 

(2.06) 

and 

II 

SI* 

—  ik  J  ii—l  R75 

0 

(npii  h  =  0). 

(2.07) 

Key: 

(1).  with. 

Further  we  must  formulate  requirement  so  that  in  the  region  near  the 
source  where  it  is  possible  to  disregard  the  curvature  of  the  earth's 
surface  and  the  curvature  of  ray/beam,  would  occur  reflecting  formula 
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for  plane  earth.  If  the  height/alt ituae  of  the  source  above  the 
earth/ground  exists  h’=b-a,  then  this  requirement  means  that  in  the 
region  indicated  must  be 


and 


c 


k  (k  ->!-)• 
2$ 


(2.08) 


.  t  (fl— ■«!*)» 

e 


-  i| 

/(  —  /l  -  S  l  l|  —  I  I  ‘ 


(2.09i 


Factors  with  second  exponential  function  represent  approximate  values 
of  Fresnel's  coefficients  for  the  vertical  and  for  the  horizontal 
polarization.  Latter/last  two  formulas  are  the  generalization  of 
formula  (1.28)  of  chapter  13. 


Let  us  note  that  expressions  (2.08)  and  (2.09)  themselves 
approximately  satisfy  limiting  conditions  (2.06)  and  (2.07). 


In  the  case  of  the  field  above  absolutely  conducting  surface 
hi  *=®)  limiting  conditions  (2.06)  and  (2.07)  take  the  form 


and 


dU.  „  .O'* 

-  0  (npn  h  =  0) 

r,  =  o  (i$h  h  =  o). 


(2.10) 

(2.11) 


Key:  (1).  with. 


and  reflecting  formulas  will  be  written  in  the  form 
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r - ,  ik  th-h')' 

v,=  1  -f-U  “ 

ik  i 

—  e  2s  , 

(2.12) 

r -  i  iklh-k'V 

* 

ik  (h+h-)‘  i 

—  e  *s  )  • 

(2.13) 

Page  290. 

3.  Analogy  with  the  unsteady  task  of  quantum  mechanics. 

The  problem  formulated  in  the  preceding/previous  paragraph  about 
the  propagation  of  wave  in  the  spherical  layer  with  the 
variable/alternating  refractive  index  presents  analogy  with  the 
quantum-mechanical  problem  about  the  motion  of  wave  packet  in  the 
assigned  force  field. 

Let  us  write  the  equation  of  Schroedinger  for  the  particle 
motion  of  mass  m0  in  the  force  field  with  potential  energy  o. 
Designating  through  x  the  coordinate  of  particle,  after  t  -  the  time, 
through  •fi  -  divided  on  2r  Planck's  constant,  we  will  have 

<30,) 

Comparing  the  equation  of  Schroedinger  (3.01)  with  the  equation  of 
Leontovich  (2.04)  or  (2.05)  for  Ua  and  for  ,  we  see  that  these 
equations  take  one  and  the  same  form,  moreover  coordinate  x  is 
proportional  to  height/altitude  h,  time  t  is  proportional  to  the 
horizontal  distance  of  s  and  potential  energy  d»  is  proportional  to 
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value  e  —  1  —  undertaken  with  the  opposite  sign  which  is 
characterized  by  from  the  so-called  given  (or  modified)  refractive 
index 

M  =  106  (1^1 4-  4)  =  (n  -  1  -  4)  <3-02) 

only  constant  factor. 

Thus,  the  parabolic  equation  of  Leontovich  for  the  amplitude  of 
stationary  process  coincides  in  form  with  the  unsteady  equation  of 
Schroedinger. 

The  resemblance  between  both  problems  is  not  limited  to  the 
coincidence  of  differential  equations,  but  it  is  spread  also  to  the 
maximum  and  "initial"  conditions. 

To  the  case  of  self-adjoint  differential  equations  and  limiting 
conditions  in  question  in  quantum  mechanics  corresponds  in  the 
electromagnetic  problem  the  case  of  the  absence  of  absorption  in  the 
air  and  in  the  earth/ground,  i.e.,  that  case  wher  the  refractive 
index  of  air  is  real,  and  the  earth/ground  -  absolute  conductor.  This 
case  is  of  greatest  interest,  also,  for  the  problem  about  the 
superrefraction.  However,  the  quantum-mechanical  methods  can  be 
generalized,  also,  in  the  case  of  the  presence  of  absorption. 

Page  291. 
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If  the  earth/ground  -  absolute  conductor,  then  limiting 
conditions  for  U2  and  for  Ws  take  form  (2.10)  and  (2.11),  and  to  them 
they  correspond  in  the  quantum-mechanical  problem  of  the  condition 

~  =  0  (npn  x  =  0)  (3.03) 

Key:  (1).  with. 


or 


Key:  (1).  with. 


( i ) 

!{•  =  0  (npii  x  =  0). 


(3.0-J) 


However,  as  far  as  initial  condit-^ns  are  concerned,  their  general 
view  consists  of  the  assignment  of  the  initial  value  of  the  wave 
function 

*  =  %  (*)  (npii  /  =  0.  0;<  x  <>=).  (3.05 1 

Key:  (1).  with. 


Function  which  satisfies  differential  equation,  initial  and 
limiting  conditions,  can  be  searched  for  '.n  the  form 

«c 

ty(x,  ()  =  f  F (x.  x'.  t)$v(x')dx‘.  (3.06) 

6 

Function  F  must  with  everyone  x’  satisfy  the  differential  equation 

-i.  2/-^  JL  _  q,f  .  0  (3.07i 

ox-  t  o:  l-  ' 

and  the  limiting  condition  of  form  (3.03)  or  £3.04)  (to  the  same  as  t}-i 
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So  that  expression  (3.06)  with  t=0  would  be  led  to  to  (*).  function 
F  must  with  t-»0  have  a  special  feature/peculiarity  whose  character 
was  connected  with  the  limiting  condition.  In  the  case  of  the 
condition 

—  =  0  (npli  x  =  0)  (3.08) 

Key:  (1).  with. 


special  feature/peculiarity  F  must  be  the  form 

_  /  t  Hi.  (»— 1*)»  .  n,  \ 

F (x.  x',  ()  —  £  4  V  -sin-l.e  -e  |* 


(3.09) 


However,  in  the  case  of  the  condition 


i‘\ 


F  —  0  (npii  x  =  0) 


(3.10) 


Key:  (1).  with. 


special  feature/peculiarity  F  must  be  the  form 

f.(«.x',0  =  e  ’lwl'  j’ 

(3.11) 

Comparing  these  formulas  with  (2.12)  and  (2.13),  we  see  that  for  the 
corresponding  limiting  conditions  special  feature/peculiarity  F  in 
the  accuracy  coincides  with  special  features/peculiarities  of  Us  and 
W, .  Actually/really,  equalizing  height/altitude  h  to  coordinate  x,  we 


must  assume 
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Expressing  F  through  the  variable/alternating  h,  h'  s,  we  will  have 
for  limiting  condition  (3.08) 

F  =  F.  {h,  h',  s),  (3.13) 

where  F2  satisfies  the  same  equation,  as  U2,  to  the  limiting 


condition 


Key:  (1).  with. 


~  =  0  (npii  h  =  0) 


and  has  the  special  f eature/pecuiiarity 

F.  (h.  k'.  s)  =  e  4  j  —  (e  **  —  e  )•  (3.15) 

For  limiting  condition  (3.10)  we  assume/set 


F  =  G.  ik.  k’  si. 


(3-16) 


where  Ga  satisfies  the  same  differential  equation,  as  war  to  the 


limiting  condition 


G.  =  0  (npi:  k  =  0) 


key:  (1).  with. 


and  has  the  special  feature/peculiarity 

.  n  _  .  *  ih-n -,»  .  k  Ih-f)1  t 

G;  ( h ,  A\  s)  =  e~  4  ]  !  ^e*  i*  -  e‘  =*  I .  (3. 18) 


We  see  that  Fa  differs  only  in  terms  of  constant  factor  from  U 
Gs  -  from  Kj ,  namely  we  have: 


and 
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L't  =  eT  ]  /-~  (3.19) 

and 

1>e'TV/5^  (3-20) 

If  we  designate  through  f(h,  s)  the  function,  which  satisfies 
the  same  equation  and  the  same  limiting  conditions  as  Ua,  and  taking 
with  s=0  the  value 

/  (h.  si  =  ft  ih)  (ripM  s  =  0),  (3.21) 

Key:  (1).  with. 

that  on  the  basis  (3.19)  we  can  write 

*  -  « 

f  ih.  S)  =  e  ’  ]  |  L'.  (A,  h\  s)  }b  (k  )  dh’.  (3.22) 

It  is  analogous,  if  f(h,  s)  satisfies  the  same  limiting  conditions 
that  also  Wa,  will  be 

.  3T  _  __  m 

f(h,  s)  =  e  "  4  j  lP.(h.  /i',  s)fb(h')dh'.  (3.23) 

Page  293 . 

Latter/last  two  formulas  are  valid  not  only  for  the  limiting 
conditions,  which  correspond  absolutely  to  the  conducting 
earth/ground  (when  there  is  an  analogy  with  quantum  mechanics),  but 
also  for  more  general/more  common/more  total  maximum  conditions 
(2.06)  and  (2.07),  moreover  then  special  features/peculiarities  U. 
and  Wa  are  given  by  formulas  (2.08)  and  (2.09). 

If  function  fe(h)  is  different  from  zero  only  in  the  region 


point  h*h',  moreover  integral  of  f, 


on  this  region  is  a  value  final, 
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then  f(  h,  s)  will  be,  with  not  too  small  s,  is  proportional  Ua  or 
respectively  Thus,  functions  Ua  and  Wa  correspond  to  point 
source,  arranged/located  on  height/altitude  h'  as  this  and  must  be. 

In  the  language  of  quantum  mechanics  it  is  possible  to  say  that 
function  v.  proportional  Us  or  K-,  represents  the  solution  of  the 
problem  about  the  deliquescence  of  wave  packet,  initially 
concentrated  in  the  vicinities  of  one  point. 

From  quantum  mechanics  it  is  known  that  the  speed  of 
deliquescence  significantly  depends  on  the  form  of  potential  energy. 
Let  us  visualize  that  the  particle  motion  is  limited  on  the  one  hand 
by  the  impenetrable  wall.  If  potential  energy  is  such,  that  the  force 
always  is  directed  from  the  wall,  then  deliquescence  occurs  rapidly. 
But  if  force  holds  particles  in  certain  region  where  the  potential 
energy  has  a  minimum,  or  near  the  wall,  against  the  deliquescence  it 
occurs  slowly  or  does  not  occur  completely.  In  this  case  the  equation 
of  Schroedingsr  admits  the  solutions,  which  correspond  to  steady 
states . 

The  wave  function  of  almost  steady  state  at  zero  time  is 
different  from  zero  only  in  the  region  of  the  minimum  of  potential 
energy.  In  the  course  of  time  the  amplitude  of  wave  function  in  this 
region  decreases  and  occurs  the  decomposition  of  initial  almost 
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steady  state.  The  decrease  of  amplitude  occurs  according  to  the 
exponential  lav,  and  the  speed  of  decomposition  is  characterized  by 
coefficient  in  the  index  which  is  called  of  disintegration  constant. 

If  initial  wave  function  itself  is  not  the  wave  function  of 
almost  steady  state,  then  in  the  resolution  of  it  it  is  possible  to 
isolate  member,  who  corresponds  to  almost  steady  state,  and  for  the 
high  values  of  time  this  member  will  be  the  main  thing. 

In  our  electromagnetic  problem  the  horizontal  distance  s  plays 
the  role  of  time  t  of  the  quantum-mechanical  problem.  To  the 
deliquescence  of  wave  packet  corresponds  the  decrease  of  the 
amplitude  of  field  with  an  increase  in  the  horizontal  distance.  The 
role  of  wall  plays  earth's  surface  h=0. 

Page  294. 

Wall  will  be  impenetrable,  if  the  earth/ground  is  absolute  conductor; 
with  the  final  conductivity  the  wall  will  be  absorbing  and  the 
decrease  of  amplitude  will  occur  not  only  as  a  result  of  wave  leak 
into  the  upper  air,  but  also  as  a  result  of  absorption  by  its 
earth/ground.  The  role  of  potential  energy  plays,  as  we  saw,  the 
reduced  refractive  index  M  undertaken  with  the  opposite  sign.  The 
course  of  the  given  refractive  index  depending  on  height/altitude  is 
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shown  schematically  in  Fig.  1. 

Unbroken  curve  showed  course  M  in  the  presence  of 
superrefraction.  The  broken  continuation  of  rectilinear  segment  of  a 
curve  corresponds  to  that  case  when  there  is  no  superrefraction  and 
it  is  possible  to  introduce  an  "equivalent  radius"  of  the  Earth, 
straight  line  proportional  to  angular  coefficient  with  respect  to 
axis  M. 

If  we  consider  the  curve  of  Fig.  1  as  the  curve  of  potential 
energy,  then  i*  will  be  clear  that  the  characteristic  for  the 
superrefraction  presence  of  maximum  for  (-M)  (minimum  for  M)  is  the 
necessary  condition  for  existence  almost  steady  state. 
Actually/really,  if  we  designate  through  A*  the  height/altitude, 
which  corresponds  to  the  maximum  potential  energy,  rhen  in  region 
h  <hm  force  will  seemingly  force  wave  packet  against  wall  and  not 
give  to  it  to  depart  to  region  A  >  hm. 

But  in  our  electromagnetic  problem  the  presence  of  almost  steady 
state  indicates  such  propagation  of  wave,  in  which  its  amplitude 
decreases  with  an  increase  in  the  distance  anomalously  slowly,  so 
that  attenuation  factor  of  its  (corresponding  to  disintegration 
constant)  is  anomalously  low.  Hence  it  follows  that  the  condition  for 
existence  Of  almost  steady  states  -  this  is  the  condition  of  the 
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possibility  of  hyperdistant  radiowave  propagation. 

The  analogy  given  here  with  quantum  mechanics  makes  it  possible 
to  compose  the  qualitative  picture  of  the  phenomenon  of  hyperdistant 
radiowave  propagation.  This  analogy  is  useful  for  those  that  on  its 
basis  some,  used  in  quantum  mechanics,  the  mathematical  methods  can 
be  transferred  into  the  region  of  radiophysics. 


Fig.  1.  Given  refractive  index  during  the  superrefraction. 

Page  295. 

On  the  other  hand,  developed  by  us  the  methods  of  the  solution  of  the 
problem  about  the  radiowave  propagation  can  find  use,  also,  in 
quantum  mechanics  1 . 

FOOTNOTE  1 .  Thus,  the  contemporary  theory  of  the  composite  moment  of 
momentum  is  directly  connected  with  our  methods  of  the  addition  of 
series/rows,  based  on  the  integration  in  the  plane  of  complex 
variable  *  (chapter  10).  (Note,  introduced  into  this  edition). 
ENDFOOTNOTE. 

However,  these  questions  emerge  from  the  scope  of  this  book. 

4.  Transition  to  the  dimensionless  quantities. 
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Let  us  return  to  the  solution  of  the  problem,  formulated  in 
paragraph  2.  Me  should  determine  functions  U,  and  Wa,  which  satisfy 
differential  equations  (2.04)  and  (2.05),  to  limiting  conditions 
(2.06)  and  (2.07)  and  to  conditions. (2.08)  and  (2.09),  that 
characterizes  special  feature/peculiarity.  This  problem  have  solved 
we  earlier  for  two  cases:  a)  unhomogeneous  atmosphere,  source  on  the 
earth/ground  and  b)  homogeneous  atmosphere,  source  was  elevated.  We 
will  now  show  that  this  problem  can* be  solved  and  for  t*e  general 
case  of  unhomogeneous  atmosphere  and  elevated  source. 

Let  us  switch  over  in  our  equations  to  the  dimensionless 
quantities,  used  in  the  preceding/previous  chapters.  For  this  let  us 
examine  coefficient  with  Uj  in  equation  (2.04).  This  coefficient  is 
proportional  tc-  value 

JLZ±  -L.  JL  =  10'6.\!  (h).  (-5.011 

where  M(h)  -  the  "modified"  refractive  index. 

We  let  us  assume  that,  beginning  from  certain  height/altitude 
h*=H,  this  value  it  is  possible  to  approximate  by  linear  function  frc 
h  and  to  place 


where  a*  -  the  so-called  equivalent  radius  of  the  Earth  and  a  - 
certain  low  constant  (for  example,  a<0.0005). 
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In  the  simplest  case  it  is  possible  to  consider  that  with  h>H 
(where  H  -  certain  height/altitude(will  be  t  =  l;  then  it  is  necessary 
to  assume  o=0  and  a*=a. 

In  that  region  where  is  valid  formula  (4.02),  equation  for  Us 
takes  the  form 


(4.03i 
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In  order  to  get  rid  of  value  a  in  the  latter/last  term,  let  us 
produce  in  it  the  substitution 

U.  =  Ceio*sMr.  (4.04) 

where  C  -  constant  by  which  we  will  be  ordered  subsequently . 


Then  equation  (4.03)  is  reduced  to  the  form 


<j/l-  OS  u 

Let  us  introduce  the  abbreviation 

feu*  \  1/3 


_  2ik  —  -  k2  ~  V  =  0. 

OS 

Dbrevial 


(4.05) 


m 


,(4.06) 


and  let  us  assume 


ks  =  2m2x,  kh  =  my ,  kh'  =  my'. 
Then  equation  (4.05)  is  written  in  the  form 


(4.07) 


(4.08) 


Oir1  Ox 

The  same  substitutions  reduce  more  precise  equation  (2.04)  to  the 
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form 

±3L-/j5L-ly-r(*)iy  =  0,  (4.09) 

oy-  ox 

where 

r(y)  =  m1!(e—  2a  —  -fr).  (4-10) 

Value  r(y)  characterizes  the  anomalous  course  of  refractive  index 
near  the  earth’s  surface;  beginning  from  certain  value  of  y,  value 
r(y)  it  can  be  placed  equal  to  zero.  If  we  consider  that  a=0  and 
a*=a ,  then 


r  (j/)  =  mz  (e  —  1).  (4.11) 

To  us  it  remains  to  express  under  new  the  variable/alternating 
limiting  conditions  and  the  conditions,  which  characterize  special 
f eature/pecul iar i ty .  Assuming/setting 


9  ~ 


im 


(4.12) 


we  will  have 


Key:  (1).  with. 


oy 


-  q'V  =  0 


(A 

(itp'H  y  =  0). 


(4.13i 


The  constant  C  in  equation  (4.04)  we  will  select  so  that  the 

equation,  analogous  (3.22),  could  be  written  in  the  form 

•» 

/(x,  y)  =  j'Hx,  y,  y')f0W)dij.  (4.14) 

o 
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Then  the  equation,  which  is  determining  special  feature/peculiarity 
T,  takes  the  form 


f* 
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—i  —r  i  .  (//—  </')’  .  (</-T 


u  —  y  ~  2n/x 


2  t'  ox  -  “  y‘  —  2/</.v 

From  comparison  (2.06)  with  (4.15)  we  obtain 


j.  <4. 15.) 


C  =  e 


‘  T  I'S-ifeT 


(4-16) 


Function  W2  differs  from  U3  only  in  terms  of  the  fact  that  under  the 


limiting  conditions  and  in  the  equation,  which  is  determining  special 


feature/peculiarity,  value  p==  is  replaced  on  1  *i— 1.  This 
corresponds  to  the  replacement  of  value  a  on 


qx  =  im  1  ij—  1.  <4.i/) 

Virtually  it  is  possible  in  all  cases  to  assume/set  qx=*. 


Together  with  function  y  we  will  examine  the  function 
V(x.y.y‘.q)  —  21  nxe  *  T.  (4.18) 

tfhich  we  will  call  attenuation  factor.  Values  U-  and  U  are  expressed 
as  V  as  follows: 

U~  =  -j-  V  (4.;‘M 

and 

i  <1— ai  Ss 

U  =  . . .  It-v.  y.  y.  q).  (4.2uj 

j/  sc  sin 

Function  W  is  obtained. from  (4.20)  by  replacement  q  on  q,. 


5.  Properties  of  the  particular  solutions  of  differential  equation. 


For  the  construction  of  function  y.  of  that  satisfyinc  the 
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conditions  presented,  it  is  necessary  to  investigate  the  properties 
of  the  particular  solutions  of  equation  (4.09),  obtained  by 


separation  of  variables.  Assuming/setting 

Y  =  e‘*'/(i/'  /)• 


(5.0h 


we  obtain  for  f(y,  t)  the  equation 


3-  +  Ur  +  »,Uf)— 


(5.02) 


DOC  *  82001203 


PAGE 


Page  298. 

Let  us  designate  through  f#(y,  t)  and  f*  (y,  t)  the  solution  of 
equation  (5.02),  which  satisfy  the  initial  conditions 

HO.  rt-l.  (5.03) 

and 

/*(0,  t)  =  0,  (5.04) 

The  general  solution  of  equation  (5.02)  will  take  the  form 

/  ( u ,  t)  =  A°fu  ( y .  /)  -  A*} *  y,  J).  (5.05) 

On  the  other  hand,  if  function  r(y)  to  increases  y  decreases 
sufficiently  rapidly,  then  with  real  t  equation  (5.02)  will  have  one 
integral  (definition  with  an  accuracy  to  the  factor,  which  does  not 
depend  on  y),  which  behaves  with  large  y  as  Wj(t-y),  and  another 
integral  which  behaves  as  Wj(t-y),  where  w,  and  w,  -  composite  Airy’ 
functions,  which  represent  the  solutions  of  the  equation 

-f  (y  —  t)  w  =  0.  (5.06) 

obtained  from  (5.02)  by  replacement  r(  y)  by  zero.  Functions  w,  and 
w,  have  the  asymptotic  expressions 
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and 


4  _4  ,  i 

n\(t—y)  =  t  (y  —  t)  e  (5.07) 

-:(f—  l/)=e  (y~i)  4  e  8  .  (5.08) 


Therefore  the  behavior  of  the  general  solution  of  equation  (5.02) 
with  y->co  and  real  t  can  be  characterized  by  the  constants  C2  and  C2 
in  the  expression 

f  (y,  0  =  Citt'j  (t  —y)  +  C,vs,  (t  —  y).  (5.09) 

Let  us  establish  connect ion/communi cat ion  between  the  constants  A0, 
A*  ,  C2  and  C2  (which  can  be  functions  from  parameter  t). 


In  view  of  equations  (5.02)  and  (5.06)  wa  have: 

=  (510) 
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In  this  equality  we  can  assume  consecutively/serially  w=wlf  then  w=w2 
and  integrate  it  from  0  to  ».  As  a  result  of  the  relationship/ratio 

(5-H) 

oy  *  oy 

we  will  have 

ta(»v§-/'f-)  =  ac,  (5.12) 

and 

=  (5.13) 

equality  (5.10)  after  integration  gives 

m 

2iC\  =  A’vhU)  -r  A’wiU)-  \r(y)f(y ,  /)  ^;(:  —  y)dy  (5.14) 


and 


*  I 

—  2iC-  =  .4°  (?)  —  I  r  (?)  f  (?•  t) -  .  (?  —  h)  d?j  - 

i  i  1 

|  r  j 

.4*  |k'j  (?)  —  J  r  (i/)  r  (?.  0  ^  — y)  ty  j  ■  (5-1;> 


Let  us  pay  now  attention  to  the  fact  that  the  coefficients  with 
A*  and  of  A*  in  these  equations  are  whole  transcendental  functions 
from  parameter  t.  Actually/really,  function  f i* ,  w, ,  w3,  they  are 
goal  functions  from  t,  whereas  integration  is  conducted  actually  in 
the  final  limits,  since  r{y)  it  is  possible  to  assume  equal  -fa  ze£-0 
beginning  with  certain  y.  (The  same  conclusion  it  will  be  correct  and 
without  this  limitation  for  r(y),  if  only  r(y)  sufficiently  rapidly 
decreases  at  infinity). 

Hence  it  follows  that  if  the  constants  A*  and  A*  will  be  whole 


.  w 
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transcendental  functions  from  t,  then  the  same  character  will  have 
the  constants  C»  and  C2.  This  makes  it  possible  for  us  to  apply 
equations  (5.16)  and  (5.17),  derived  for  the  case  real  t,  also  in  the 
case  of  the  arbitrary  composite  values  t. 
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and 


If  we  will  assume 

•• 

y?°  =  (/)  =  ic'i  (/)  —  j  r  ( y )  f  ( y •  t)  tt*i  (/  —  y)  dy  t5-18) 

6 

•t 

x  =  A\  ij)  =  —  C'i  (? »  -  j  r  (y)  f  t)  S’;  0  —  .V>  (riA9) 


the  expression 

/,  [y.  t)~  A?  (!)  f  It/,  i)  -  A]  (/)  f  (y.  t)  (5.20.1 

will  be  such  solution  of  equation  (5.02),  which  behaves  with  y->®  as 
WjCt-y),  and  it  represents  at  the  same  time  whole  transcendental 
function  from  t.  It  is  analogous,  if  we  will  assume 

m 

A~  =  .4'J  1  =  «■*;  (/ )  —  \r(y)  f  (y.  t)  a*  (;  —  y)  dy  (5.21 ) 


and 


A  =  An(f)  =  —  (/)  -r  I  r (y)  f  (y.  1)^(1  —y)dy,  (5.22) 


the  expression 


f.  ( y -  ■'.»  -  A'i  U) r (y.  i)  -  a\ (I)  f  (y.  J)  (5.23., 


will  behave  with  as  w2(t-y),*and  it  will  be  integral  function 
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from  t. 


Integral  ft(y,  t)  will  have  the  asymptotic  expression 

.24) 


m  is-  n  =  — 


c  t 


.  si 

i  — 

4 


■exp 


i  1 1  u  —  t  —  r  (u)du 

T 


,  (5.2 


»  V  —  I  -f  r  (V) 

while  integral  f2(y,  t)  will  have  the  asymptotic  expression 


7s  W-  <>  =  - 


r  «■ 


•exp 


l  y-t-  r[tj) 

where  c‘,  c",  r  -  constants. 


—  i  J  la-/  -r  r  (u)  du\,  (5.25) 


j,  (5-5 


If  we  place  r(y)-0,  r=z ,  c'=c"  =1,  then  formulas  (5.24)  and 
(5.25)  will  pass  into  asymptotic  expressions  (5.07)  and  (5. OB)  for  w 
and  w,. 


Integral  f2(y,  z)  we  already  used  in  Chapter  13,  where,  however 
it  was  accepted  without  the  proof  that  there  is  such  integral,  which 
allows/assumes  asymptotic  representation  (5.24)  and  is  at  the  same 
time  whole  transcendental  function  from  t. 

i 
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Formulas  ( 5 . 18 > — ( 5 . 23 ) ,  derived  in  the  course  of  proof  of  this 
confirmation,  can  be  used  also  for  the  actual  (numerical) 
construction  of  this  integral. 
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Function  f  x ( y ,  t)  with  composite  t  will  grow/rise  with  increase 
y,  and  integral  of  square  fx(y,  t),  undertaken  on  y  from  0  to  »,  will 
be  divergent.  However,  under  some  assumptions  about  behavior  r(y)  in 
'-.he  composite  plane,  function  fj(y,  t)  will  behave  with  composite  y 
as  w^t-y),  and  will  vanish  on  ray/beam  y  =  rela  (where  a  =  so  that 
the  integral 

/  =  (  f\(y.t)dy  (5.26) 


will  already  converge.  Let  us  ce?culate  the  value  of  this  integral. 
Differentiating  equation  (5.02)  for  t,  we  will  obtain 

|r  ('£)-!!'  —  '<*>-'>&-/•  l5-27» 

Hence  and  from  (5.02)  we  obtain  the  relationship/ratio 
(l  =  \pdu.  (5. 28) 

Assuming/setting  here  f=fi(y,  t)  and  counting  upper  integration  limit 
by  equal  to  coe-,  we  will  have 


>  ...  (,  &'■! i  dj i  ft,  \ 

j  / 1  (y>  0 dy  —  \h  fad;  dt  oy  l 


6.  Construction  of  solution  in  the  form  counter  integral  and  in  the 
form  of  series/row. 
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In  the  preceding/previous  paragraph  we  established  the  existence 
of  two  integrals  of  the  ordinary  differential  equation 

-0--  l.M  —  r  (y)  —  l\  /  ==  0,  (6.01) 

which  are  whole  transcendental  functions  from  parameter  t  and  have 
asymptotic  expressions  (5.24)  and  (5.25).  Integrals  these,  designated 
by  us  through  f.(y,  t)  and  f*(y,  t)  are  determined  by  formulas  (5.20) 
and  (5.23). 
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We  will  show  now  that  with  the  aid  of  functions  f,  and  f ,  it  is 
possible  to  construct  for  the  V  and  for  ¥  contour  integral,  which  is 
the  solution  of  our  problem.  (Reason  ingsy5ur)  it  will  be  analogous  to 
the  reasonings,  presented  in  the  paragraph  3  o £  Chapter  13,  and 
final  formulas  will  present  analogy  to  formulas  (2.24)  and  (3.10) 


Chapter  12. 


Let  us  designate  through  D^(t)  Wronskian: 


d„<o /«■£■ 


and  it  is  placed 


f  (/.  y<  y'<  <?)  = 


=  rrW.fe'.o!f,<!/.  D-  ?■(»•<)]■  <«•« 


h  (0.  0  +  ?/j  (0,  t) 
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with 

where  primes  a  fi  and  f ,  designate  derivatives  on  y. 

Me  will  consider  that  y'>y,  and  let  us  compose  the  integral 
T  =  [  eix'F(t,  1 1,  y',  q)di,  (6.04) 

undertaken  on  the  duct/contour,  which  covers  in  the  positive 
direction  all  of  the  pole  of  integrand. 


From  the  determination  of  function  F  it  follows  that  it  is 
meromorphic  function  by  the  variable/alternating  t  (i.e.  it  has  with 
final  t  of  no  special  features/peculiarities,  except  poles).  Function 
F  is  completely  determined,  even  if  the  entering  it  functions  f,  and 
fj  are  determined  only  with  an  accuracy  to  the  factors,  which  do  not 
depend  on  y.  Since  integrals  fa  and  f,  all  values  t  are 

independent  (this  is  evident  from  their  asymptotic  expressions),  then 
Wronskian  DXI(t)  roots  does  not  have,  and  the  unique  poles  F  are  the 
roots  of  the  equation 

f'i  (0-  t)  -r  <?/,  (0.  /)  =  0.  (6.05) 

If  in  differential  equation  (6.01)  function  r(y)  is  equal  to  zero, 
then  it  is  possible  to  assume 

U  to.  t)  =  ©J  (t  —y),  f.  to,  t)  =  10,  (t  —  y).  (6.06) 


Then 


Dlt(t)  =  -2  i 


(6.07) 


and  expression  (6.03)  for  F  is  led  to  that  examined  in  Chapter  12 
[formula  (2.21)3. 
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Let  us  show  that  function  Y,  determined  by  contour  integral 
(6.04),  satisfies  all  conditions  presented. 

Page  303. 

It  is  first  of  all  obvious  that  function  Y  satisfies 
differential  equation  (4.09),  since  it  satisfies  expression  under  the 
integral.  Further,  function  *}f  satisfies  limiting  condition  (4.13), 
since  with  all  t  and  y'  we  have 

—  qF  =  0  (npii  y  =  0).  (6.08.1 

Key:  (1).  with. 

Us  it  remains  to  show  that  function  T  has  necessary  special 
f  eature/pecul i ar i ty . 

With  the  aid  of  asymptotic  expressions  (5.24)  and  (5.25)  it  is 
possible  to  show  that  if  x,  y  are  small,  and  ratio  y/x  is  great,  then 
the  main  section  of  integration  in  (6.04)  will  lie/rest  at  the  high 
negative  values  of  t.  But  if  t  is  great  and  is  negative,  then  in 
differential  equation  (6.01)  main  role  in  the  coefficient  with  f  will 
play  term  --t.  Therefore  with  large  negative  t  we  will  approximately 
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have 

itO-rte**'-*  <6’09) 

and 

fsto.0~fs(°-  (6‘10) 

After  substituting  these  expressions  into  formula  (6.03)  for  F,  we 
will  obtain 

r _ £ _ lei  u-'-v)  *  -<  _  q—'KSe*  «»*+«  ,  rT !  •  (6-1 1) 

F--tyT7f  q-riV-i  ) 

The  substitution  of  this  value  F  into  integral  (6.04)  gives  for  Y 
the  formula  of  Weyl  -  van  der  Polya,  who  after  neglect  of  the  low 
values  (in  the  second  term)  is  led  to  expression  (4.15),  which 
characterizes  special  feature/peculiarity  Y. 

Thus,  the  validity  of  expression  (6.04)  for  V  we  have 
established/ installed. 

From  contour  integral  (6.04)  is  not  difficult  to  switch  over  to 
the  series/row,  arranged/located  on  the  deductions,  relating  to  the 
roots  of  equation  (6.05).  Let  us  write  this  equation  in  somewhat  more 

detail. 

Utilizing  expression  (5.20)  for  fx(y,  t)  and  initial  values 
(5.03)  and  (5.04)  functions  f°  and  f*  we  obtain 

M0,  =  /i(0,  i)  =  Ai(t)  (6-12) 
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and  equation  (6.05)  takes  the  form 

A\(i)~qAUt)  =  0.  (6.13) 

Substituting  here  values  (5.18)  and  (5.19)  functions  A®  x  and  A*x,  we 


will  have 


-  J  r  to)  \f° to ,t)~qf(y,t )1  »,  (/ -y)dy  =  0.  (6.14) 


This  equation  we  will  call  characteristic. 
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For  us  it  is  significant  that  the  left  side  of  the 
characteristic  equation  is  whole  transcendental  function  from  t  and 
that  it  contains  only  functions  f#(y,  t)  and  f*  (y,t)  which  can  be 
obtained  for  all  values  of  t  by  the  numerical  integration  of 
differential  equation  (5.02)  with  initial  conditions  (  03)  and 

(5.04).  When  function  r(y),  beginning  from  certain  y=yx,  is  equal  to 
zero,  integration  in  (6.14)  can  be  carried  out,  and  characteristic 
equation  is  reduced  to  the  form 

tt’i  U  —  y)  l/°  ( y ,  0  —  qf  to<  0]  + 

-  a-,  (?  y)  if0  to'  0 — ql*  to'  01  =  o 

(,  > 

(npii  y  —  «,)•  (6.15) 

Key:  (1).  with. 


Characteristic  equation  for  the  case  of  homogeneous  atmosphere  takes 


UrtlUiftMBMV  *,,4.1  ,,.1,4, ,  u  Mty  lui  jl  \jLk  (j,  j 
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the  form 


(:i  —  <7-1  (t)  — 


This  equation  is  obtained  from  the  preceding/previous  formulas,  if  we 
in  (6.14)  place  r(y)=0  or  if  we  in  (6.15)  assume  y^y^O. 

The  roots  of  characteristic  equation  we  will  designate  through 


h  (<7).  U  (4).  •  • 


(6.17) 


These  roots  will  be  functions  from  parameter  q. 


Let  us  turn  to  the  calculation  of  the  deductions  of  integral 
(6.04).  From  equations  (6.02)  and  (6.05)  it  follows  that  with  y-0  and 


t^t^  will  be 


?;<o.  /»-«/;  io.  n  _  i 
“ — D^TTi  f>  <0'  »>  ‘ 


(6.18) 


Further,  derivative  on  t  of  the  denominator  in  (6.03)  can  be 
represented  in  the  form 

-zk  -  =  •  <6I9> 

Therefore  the  deduction  of  function  F  at  point  ft,  will  be  equal  to 


dti  /,  (y.  /»)  Ij  <Mi  hi 
dq  /,  <0.  *,»  /,  (0.  <»)  ‘ 


(6.20) 


Taking  the  sum  of  expressions  (6.20),  multiplied  by  e  ,  we  vAll 
obtain  the  unknown  expansion  of  function  4'  in  the  senes/rov 
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>r-y 


»=i 


jxts  dts  fi  (y  .it)  ft  ((/■  *«) 

6  dq  /,( 0.  I*)  /,  (0,  /,) 
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Values 


4Kf  «/•(».« 


(6.21) 


(6.22) 


can  be  named/called  hi  cjh-altitude  factors.  Let  us  note  that  the 
high-altitude  factors  are  expressed  according  to  (6.22)  through 
functions  f®  and  f*  which  are  directly  calculated  by  the  numerical 
integration  of  differential  equation  (5.02). 


When  value  q  is  very  great  or  equal  to  infinity  (horizontal 
polarization,  absolute  conductor),  formula  (6.21)  must  be  converted 
via  term-by-term  multiplication  by  value 

*•"  i  (->.  (,) 


Result  can  be  written  in  the  form 


(6.23) 


a* 

T-  V  e'"‘  (V 44 )  4^  —  •  (6  24 1 

Zi  aq  ;i(u.  g 


Value 


**  "ST  ”  d  /,(0.  M  , 
di  \  to.  ft  • 


(6.25) 


will  have  with  q->®  finite  value.  Let  us  note  that  from  formulas 
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(6.19)  and  (5.29)  escape/ensues  the  relationship/ratio 

K  C‘Q 

7=  f  f\&t)dy  (o  =  -2-).  (6.26) 

6 

Therefore  series/row  (6.21)  can  be  written  in  the  form 

(6.27) 


tjr  _  f\  (y  •  /i  (y.  ^ 

/  i 

|  (lu.U)dy 


In  this  form  it  resembles  eigenfunction  expansion.  In  series/row 
(6.27)  the  "eigenvalues"  are,  however,  composite  and  the 
standardizing  integral  standing  in  denominator  converges  only  with 
the  composite  way  of  integration. 


In  order  to  switch  over  from  function  Y  to  attenuation  factor 

V,  it  suffices  to  recollect  relationship/ratio  (4.18) 

.  « 

V(x,y.y\q)  =  2ViTxt  4  V. 


Page  306. 

In  the  case  of  homogeneous  atmosphere  when  it  is  possible  to 
place  r(y)=0  and  fj(y,  t)=Wj(t-y),  expressions  escape/ensuing  from 
our  formulas  for  V  are  led  to  those  such  as  were  derived  in  (Chapter 
12  by  another  method. 

7.  Use/application  of  general  theory  to  the  case  of  superrefraction. 
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(Schematic  example). 

• 

Examined  in  the  paragraph  3  analogies  with  the  unsteady  task  of 
quantum  mechanics  make  it  possible  to  compose  for  itself  the 
qualitative  picture  of  the  phenomenon  of  superrefraction  and  those 
conditions  under  which  this  phenomenon  can  occur.  On  the  other  hand, 
the  obtained  in  paragraph  6  general/common/total  expression  for 
attenuation  factor  is  useful  for  the  quantitative  calculations  which, 
true,  require  sufficiently  complicated  calculations. 

Lex  us  write  expression  for  function  M'  connected  with  the 
attenuation  factor.  Assuming/setting  for  the  brevity 

f{y.  0  —  (701) 

we  will  haMKon  the  basis  (6.21) 

«c 

V  =  £  e*'z,»  -g-  /  (§r\  t,)  f  (y,  tj,  (7.02) 

where  of  value  they  are  the  roots  of  transcendental  equation 
(6.14). 

If  r(y)=0  with  y>y,,  this  equation  can  be  according  to  (6.15) 
written  in  the  form 

®« 0  —  y%)  f  tei-  0  +  (t — y  i)  Y  (y,.  o  =  o.  (7.03) 

where  w*  >  indicates  derivative  on  argument  (t-y),  and  derivative  on 
y.  Parameter  q  enters  into  this  equation  through  value  (7.01). 
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The  determination  of  the  conditions  under  which  the  hyperdistant 
possible  propagation,  is  reduced  to  the  investigation  of  the  roots  of 
characteristic  equation (6. 14)  or  (7.03).  In  the  absence  of 
superref ract ion  the  alleged  part  of  the  roots  of  this  equation,  which 
according  to  (7.02)  gives  weakening  wave  with  an  increase  in  the 
distance,  there  will  be  of  the  same  order  as  real  part.  However,  in 
the  presence  of  superrefraction  there  are  one  or  several  roots  with 
the  anomalously  small  alleged  part. 

In  order  to  compose  to  itself  representation  about  that,  under 
what  conditions  can  occur  hyperdistant  propagation,  let  us  examine 
the  followin^schematic  example. 

Page  307. 

Let  function  r(y)  take  the  following  form: 

r(y)  =  (l-i-|»*)  (*»-*)  (h|h  0<y <*,).  1 
r{y)  =  0  (npH  H,<y).  I 

Key:  (1).  with. 

This  corresponds  to  the  assumption  that  the  graph/curve  of  the  given 
refractive  index  represents  the  broken  line,  depicted  in  Fig.  2. 
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If  we  consider  that  dielectric  constant  f  is  changed  according 
to  the  law 


ft) 


e  =  1 

Key:  (1).  with. 


e=  1  —gih  —  hj)  (npw  h  <  hx).  \ 
(npH  1 


(7.05) 


that  the  parameters  u3  and  y!  will  be  equal  to 
u*  =  -f-l,  »,  =  *,{  If.  (7.06) 


Thus,  the  parameter  u  on  wavelength  does  not  depend,  bur  parameter  yx 
(the  reduce  height  of  crank  point)  will  be  proportional  a-*/®. 

Equation  for  f  is  written  in  the  form 

J^L-Ui-r  ns) wi -  fy - M f  =  0  <nPI!  y < I 

**  * .  O  (7.07, 

—  —  0  (npu  I 

du- 


Key:  (1).  with. 


Let  us  introduce  instead  of  t  the  new  parameter 

*  _  /  — (i  -  m’IJ/i  _ 

so - p -  (*-08, 

and  instead  of  y  -  new  variable/alternating 


“  «•  +  vs- 


(7.09) 
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Fig.  2.  Schematic  course  of  the  given  refractive  index. 

Page  308. 

To  value  of  y-y.  will  correspond  value  £=£,,  where 

us|,  =  /  —  (7.10) 

Equation  (7.05)  will  take  the  form 

~<i§r  —  *f  (io  <C  5  <T  !i)-  (7. 1 1 ) 

Its  independent  solutions  it  will  be  Airy’s  function  u(£)  and  v{£). 
Functions  f*  and  f*  will  beequal  to 

/•'  < u)  =  i s  (lot  :■<;)  —  v‘  do)  u  d)  j 

and  '  (7.12} 

r  tt,)  =  -y  iv  do»  u  d)  -  u  do)  t-  an-  | 


In  view  of  the  relationship/ratio 

u' d)*  d>  -  i"  d) «  d>  =  l  (7.13) 

of  function  f*  and  f*  they  will  satisfy  initial  conditions  (5.03)  and 
(5.04).  Introducing  according  to  (7.01)  the  function 


i 
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i  u-i  =  — fif'  Soil  “  (I)  - 


-  —  iau  (Ic'i  -  uu‘  (1*1!  v  (i) , 


(7-14) 


ve  vill  obtain  characteristic  equation,  if  ve  substitute  values  of 
f(y)  and  f’(y)  with  y=yj  into  formula  (7.03).  This  characteristic 
equation  can  be  written  in  the  form 


u:-‘  ij*i  — 

m.'  i|*.  —  qu  i|4i 


Let  us  assume  that  value  y,  and  parameter  u  is  sufficiently  great. 
This  means  that  the  "potential  pit"  in  Fig.  2  is  sufficient!/ vide  and 
deep.  In  that  case  of  value  £,  and  v1%x  [arguments  of  the  functions, 
entering  the  right  side  (7.1S)3  they  vill  be  great.  In  view  of  the 
asymptotic  expressions 


_i_  j  _i.  _JLt3;2 

u(a  =  E  •‘e3*  ,  t'(|»  =  -vl  *  e  3’  (7.16) 

the  right  side  of  equation  (7.15)  will  be  very  small,  and 
characteristic  equation  approximately  is  led  to  the  following: 


u  (1#)  *:•  qv  =  0. 


(7.17) 


Page  309. 

This  the  case  will  occur  when  in  the  sufficiently  large  region  of 
heights/altitudes  the  gradient  of  tha  dielectric  constant  of  air  will 
be  negative  and  greater  than  2/a,  where  a  -  radius  of  the  Earth;  then 
the  curvature  of  ray/beam  will  be  more  than  the  curvature  of  the 
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Earth  and  formally  computed  "equivalent  radius"  is  negative. 

The  attenuation  of  wave  with  an  increase  in  the  horizontal 
distance  is  connected  with  the  alleged  part  of  value  t,  and  also, 
therefore,  with  the  alleged  part  of  value  £0;  if  £0  they  was  real, 
fading  would  be  absent.  Fading  can  originate  from  two  reasons:  from 
the  absorption  in  the  earth/ground  and  from  wave  leak  into  the  uppe 
air.  Absorption  in  the  earth/ground  is  characterized  by  the  composite 
parameter  q.  Equation  (7.17)  corresponds  to  that  case  when  fading 
occurs  only  due  to  the  absorption  in  the  earth/ground.  If  we  consider 
the  earth/ground  as  absolute  conductor,  it  is  necessary  to  assume  q=0 
for  the  horizontal  polarization  and  q=®  for  the  vertical 
polarization.  With  q=0  equation  (7.17)  is  reduced  to  the  form 

v  (l„)  -  0  (ii fpii  q  -  Oi. 

Key:  (1).  with. 

Its  roots  will  be  the  real  negative  numbers 

S0  =  — 1,019;  —3.248.  -4.820;  .  .  .  (~-19) 

With  q=®  equation(7.17)  takes  the  form 

.  '<■’>  .  (7.20) 

r(S*)»0  inpH9=*«)  1 


Key:  (1).  with. 


DOC  *  82001203 


PAGE 

and  it  has  by  the  roots  of  the  number 

=  -2.338;  -4,088;  -5.521;  .  .  .  (7.21) 

Since  in  these  cases  of  value  they  are  real,  the  fading  is  absent. 

Equation  (7.17)  will  present  a  good  approximation/approach  to 
(7.15),  if  value  £*  (or  its  real  part)  is  positive  and  sufficiently 
great.  Since 

fi  =  too  “ 

the  this  condition  will  cease  to  be  made  beginning  from  certain  root 
£0.  Therefore  a  number  of  roots  with  a  small  alleged  part  will  be 
final. 

It  can  derive  the  approx imatio^formula  for  the  correction  to 
value  {o,  obtained  via  the  account  of  right  side  (7.15).  Let  us 
designate  by  £oe  the  root  of  equation  (7.17),  which  we  will  consider 
as  the  defective  value  £„,  and  through  -  correction.  This 
correction  will  be  obtained,  if  we  substitute  into  right  side  (7.15) 
approximate  value  £  x ,  equal  to 

*  *  i 

vi  *"  vim  v  MI/i- 
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Approximate  value  of  correction  will  be  obtained  then  from  the 
equation 
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(7-22, 

where 

S  =  i-tf+l)ir.  (7.23) 

Let  us  note  that  the  alleged  part  of  the  correct ioA'is  positive. 
This  corresponds  to  that  fact  that  the  leakage  into  the  upper  layers 
increases  fading. 


-  I 


it 
3  *1 


J/'J 


|  10  !' 


Applicability  condition  for  these  approximation  formulas  is 
sufficiently  high  value  tty!.  Let  us  recall  that  according  to  (7.06) 
we  have 


(7.24) 


where  g  -  undertaken  with  the  opposite  sign  gradient  of  dielectric 
constant;  a  -  radius  of  the  Earth  and  hx  -  spot  height  of  fracture  in 
Fig.  2. 


The  greater  the  value  nylr  the  greater  the  number  of  almost 
steady  states  with  a  small  fading. Tentatively  it  is  possible  to  say 
that  a  number  of  such  states  is  equal  to  a  number  of  roots  which 

do  not  exceed  (in  the  absolute  value)  the  parameter  vyl. 


The  concept  about  the  ray/beam,  which  is  reflected  from  the 
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upper  bound  of  layer  and  from  the  earth’s  surface,  begins  to  become 
applicable  if  and  only  if  a  number  of  almost  steady  states  [number  of 
terms  of  series/ row  (7.02)  with  a  small  fading]  it  becomes  large. 
Generally,  the  necessary  condition  for  the  applicabiliteof  the 
concepts  of  geometric  optic/optics  is  slow  convergence  of  series 
(7.02),  when  in  it  plays  a  r»‘&  a  large  number  of  members.  But  if  in  it 
are  essential  one-two  members  (which  can  correspond  both  almost  to 
stationary  and  damped  states),  then  the  concept  of  ray /beam  is 
applicable  completely. 

^Approximation  formulas  for  the  terms  with  a  small  fading. 

Using  method,  by  the  analogous  to  the  volume,  such  as  it  is 
applied  in  quantum  mechanics,  it  is  possible  to  derive  approximation 
for  the  high-altitude  factors,  which  correspond  to  terms  with  a  small 
fading,  and  to  also  give  the  estimation  of  that  part  of  the 
attenuation  factor,  which  corresponds  to  leakage  into  the  upper 
layers. 

Let  us  assume  in  equation  (6.01) 

y  +  r  (y)  ~  p  ly)  (8-01) 

and  let  us  write  this  equation  in  the  form 


(8.02) 
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In  the  case  of  superrefraction  function  p(y),  proportional  to  the 
given  refractive  index,  it  will  have  a  minimum,  and  on  both  sides 
from  it  it  will  grow/rise;  to  the  left  of  the  minimu/jthe  greatest 
value  p(y)  will  be  p(0),  and  to  the  right  p(y)  will  increase  as  does 
value  y.  If  parameter  t  lies/rests  between  the  least  by  value  p(y) 
and  by  value  p(0),  then  coefficient  with  f  in.to  equation  (8.02)  will 
become  zero  at  two  values  of  y  which  we  will  designate  through  yx  and 
y,.  In  gap/interval  yi<y<ya  value  p(y)-t  will  be  negative,  and  out  of 
this  gap/ interval  -  positive. 

In  gap/ interval  yi<y<y1  the  solution  of  equation  (8.02)  can  be 
approximately  expressed  t'nroughthe  Airy's  functions.  Let  us  assume 

j  V7=m  dy-Jt  (8-()3) 

Vi 

and 

j  V't—piy) dy  =  —  Is"  (8.04) 

V 

and  designate  through  S  sum  of  these  values  which  does  not  depend  on 

y:  _ 

bt  „ 

s  =  j  l  !  —  p  (y)  dy- 

bi 

Value  S  we  will  consider  it  large.  With  such  designations  we  will 
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approximately  have 


f=\  7 


4/ 


■p<y) 


and  also 


moreover 


i/  — h—~  [Am  (y  -  Btv  d-)]. 
>  I  /—/My) 


(8.06) 

(8.07) 


±t 4-  i-tf  =  S  (8-OS) 

and  the  constants  A,,  Ba,  Aa,  Ba  are  connected  with  the 
relationships/ratios 

=  ~  Bfi-s,  B.  =  2A1es,  (8-09| 


which  escape/ensue  from  the  comparison  of  asymptotic  expressions  for 
(8.06)  and  (8.07)  at  the  high  values  Si  and 


Page  312. 

With  y>yj  we  can  determine  valu  Ss  by  means  of  the  equality 
I  I  p  (ii)  —1  dy  =  -j-  (—  li?  -  (8  - 1 0) 

i* 

and  use  previous  expression  (8.07)  for  f. 

It  is  analogous,  with  y<yi  we  can  instead  of  (8.03)  assume 
j  \7w=idy  =  ±i-illF  (8.11) 

V 

and  apply  for  f  expression  (8.06). 
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Let  us  select  the  constants  A,  B  so  that  the  function  f  would  be 
proportional  f*(y,  t).  We  must  assume 

A.  =  C,.  B-  =  /Cj  (8.12) 

anct,  therefore, 

A^  +  Cie-s.  Bl  =  2C1es.  (8.13) 

Then  formulas  (8.06)  and  (8.07)  take  the  form 


and 


U  :)  =  2C,es  -| 


(8.14) 


'!<'/■  -1=  C,  1 


(8.15) 


Analogously  are  obtained  the  following  approximations  for  f*(y,  t): 
/.,(«.  ? )  =  2C.es  ;t.(li)--(-e-2S«(I1)]  (8.16) 

and 

fiilJ'  0  —  }'  ^*)*  (8-17) 

In  this  approximation/approach the  Wronskian  Dj,(t)  proves  to  be  equal 
to 

£>,.  (0  =  ~2iC1Ct.  (8.18) 

With  y<y*  functions  v(  {x)  and  u(^)  will  be  one  order. 

Page  313. 
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As  a  result  of  the  smallness  of  factor  e-25  the  second  terms  in 
(8.14)  and  (8.16)  vill  present  small  corrections  [generally  speaking 
it  is  less  than  an  error  in  entire  expression  (8.14)  or  (8.16)]. 
Therefore  in  region  y<yi  functions  fA  and  fa  vill  be  almost 
proportional  to  each  other. 

Throwing/rejecting  small  corrections,  it  is  possible  to  write, 
equation  for  determining  of  t  in  the  form 


We  designated  here  through  £,  and  (~t)(  the  values  and  iiL  with 
y=0. 


This  equation  is  analogous  with  equation  (7.17).  It  gives  only 
that  part  of  the  attenuation  factor  of  wave,  which  originates  from 
absorption  in  the  earth/ground.  Since  the  composite  parameter  q, 
which  characterizes  the  properties  of  soil,  is  known  only  very 
roughly,  then  coefficient  (-^)c  it  suffices  to  take  in  the  rough 
approximation  and  to  assume  accordingly  (7.06)  and  (7.09) 


(§)„  =  "  =  f?->- 


(6.20) 


where  a  -  radius  of  the  Earth  and  g  -  undertaken  with  the  opposite 
sign  gradient  of  dielectric  constant. 


Then  equation  (8.19)  is  reduced  to  form  (7.17),  investigated  in 


the  preceding/previous  paragraph.  The  roots  of  equation  (8.19) 
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will  be  connected  with  the  appropriate  values  of  parameter  t  with  the 
relationship/ ratio 

x _ 

*f]  -t-l  -  -~-dh  (8.21 1 

»• 

where  h*  x  -  smaller  of  those  two  values  of  h*x  and  h*x 
height/altitude  h,  at  which  a  subradical  quantity  becomes  zero. 

If  if  is  real,  then  h*x  and  t  are  real;  if  £,  it  is  composite, 
then  the  calculation  of  integral  (8.21)  requires  the  analytical 
continuation  of  the  interpolation  of  formulas  for  *  into  complex 
domain. 


The  necessary  condition  for  the  applicability  of  the 
preceding/previous  formulas  is  the  smallness  of  value  e-^.  where  S 
has  value  (8.05).  In  usual  unity  the  integral,  which  expresses  S, 
will  be  written  in  the  form 


1 1 . 
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After  determining  t  from  (8.21),  it  is  necessary  to check  that 
integral  S  for  this  t  is  sufficiently  great. 
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In  the  case  of  absolute  conductor  (q=0  and  q=«)  approximate 
values  and  t  are  from  (8.19)  and  (8.21)  real.  In  this  case  it  is 
possible  to  indicate  approximate  value  of  the  alleged  part  of  the 
correction  to  £,. 


Assuming/setting 


;u  =  f/.  —  <;<•, 


(8.23) 


we  will  have 


|  St"  =  4-  e--\ 


On  the  conclusion/output  this  formula  we  stop  will  not  be. 


Since  ll  -  low  value,  then  to  increase  A{0=i£"  ,  will 
correspond  increase  A/  =  //*  =■  JL  But  value  (8.24)  multiplied  by  i, 

c»o 

presents  an  increase  in  integral  (8.21).  Therefore  we  can  determine 
t"  (alleged  part  t)  from  the  equation 

1  ”"(*  i  ^  ~  ~  6  ~  1  '  -T"  dh)  =  — -f  e_iS-  (8-25) 

Since  derivative  of  integral  is  negative,  then  for  t"  is  obtained 
positive  value,  which  corresponds  to  fading. 


The  obtained  formulas  make  it  oossible  to  derive  also 
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approximation  for  value  According  to  (6.19)  we  have 


CC  _  S'-  ig 

<1:  o?  ot 


(8.26) 


Substituting  here  value  +x  from  (8.14)  and  disregarding  low 
values,  we  will  obtain 


_  (4±\  dj*  (  f'-  foi  f'  g.\  \ 
lit  \  du  Jo  dt  \  r  (j,)  /  ■ 


(8.27) 
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Here  we  can  assume  in  auto  the  rough  approximation 

,  si 


{it),.**!1'  ~ir  =  -p-  (8.28) 

[see  formulas  (7.08)  and  (7.09)].  Using  differential  equation  and 
limiting  conditions  for  v,  we  will  obtain  hence 


sq 

ji 


c- 

M* 


(fe.l'iM 


The  first  terms  of  series/row  (6.21)  for  S',  possessing  a  small 
fading,  will  be,  in  our  approximation/approach,  they  are  equal  to 


2j  Z, 


t“<  ,  .  (8.30) 


***  ii*.'  —  s.^  (i* 


where  it  relates  to  reduced  height  y’. 


If  it  is  so  great  in  the  absolute  value  that  it  is  possible 
to  use  for  *?(£•>  and  *  v  (l»V  asymptotic  expressions,  then  denominator 
in  this  formula  will  be  approximately  equal  to 
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»**<£.» (8'31) 

In  conclusion  it  is  necessary  to  emphasize  that  formulas  derived 
in  this  paragraph  are  based  on  the  sufficiently  rough  approximations 
and  are  intended  for  the  tentative  calculations.  More  precise 
calculations  must  be  based  on  a  strict  theory,  presented  in  the 
preceding/previous  paragraphs. 
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Chapter  15. 

Approximation  formula  for  the  horizon  distance  in  the  presence  of 
superrefraction  1 . 

FOOTNOTE  1 .  Foch,  1956.  ENDFOOTNOTE. 

On  the  basis  of  the  general  formulas,  obtained  in  Chapter  14, 
are  investigated  the  cases  of  normal  refraction  and  superrefraction. 
For  the  functions,  which  stand  under  the  contour  integral  in  the 
expression  for  the  attenuation  factor,  are  derived  the  approximations 
by  the  asymptotic  integration  of  the  corresponding  differential 
equations.  In  the  case  of  normal  refraction  these  functions  are 
approximated  through  the  Airy's  functions,  while  in  the  case  of 
superrefraction  -  through  the  functions  of  parabolic  cylinder.  The 
qualitative  investigation  of  contour  integral  leads  to  the 
approximation  formulas  for  that  distance,  beginning  from  which  field 
rapidly  decreases,  i.e,,  for  the  horizon  distance.  During  the 
superrefraction  when  there  are  many  waves  reflected,  under  the 
"horizon  distance"  is  understood  the  distance  for  the  first  wave 
reflected.  General  formulas  are'  applied  to  the  case  when  the  given 
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refractive  index  depends  on  height/altitude  according  to  the 
hyperboli  law. 

1.  Introduction. 

In  Chapter  14  is  derived  for  the  attenuation  factor  general 
formula  in  the  form  of  contour  integral.  The  obtained  expression  is 
applicable  for  the  very  general  case  of  the  arbitrary  course  of 
refractive  index  depending  on  height/altitude.  Basic  difficulty  of 
applying  our  general  formula  consists  of  the  solution  of  differential 
equation  for  the  high-altitude  factor.  This  of  difficulty  can  be 
avoided  by  using  the  asymptotic  solution  of  the  equation  (this 
reception/procedure  is  based  on  the  presence  in  equation  of  the  high 
parameter) . A  f ter  obtaining  approximations  for  the  high-altitude 
factor,  it  is  possible  to  write  in  the  explicit  form  integrand  in  the 
contour  integral  and  then  it  to  investigate. 

Page  317. 

The  qualitative  investigation  of  integrand  makes  it  possible  to  give 
the  estimation  of  those  distances  beginning  from  which  attenuation 
factor  begins  rapidly  to  decrease,  in  other  words,  the  estimation  of 
horizon  distance. 
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2.  Initial  formulas. 

As  shown  in  Chapter  14,  field  from  the  vertical  and  horizontal 

electrical  and  magnetic  dipole  is  expressed  in  the  general  case 
through  two  functions  of  /^ertz:  U  and  Wj which  satisfy  identical 
differential  equations;  limiting  conditions  for  U  and  W  also  of  the 
identical  type,  but  with  the  different  values  of  coefficients.  Each 
of  the  functions  of/^ertz  can  be  expressed  through  the  attenuation 


In  the  tasks,  connected  with  the  superrefraction,  equivalent 
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tiX 


radius  of the  Earth  does  not  play  that  role  such  as  it  plays  in  the 
case  of  normal  refraction;  therefore  we  it  do  not  here  introduce. 
Besides  the  enumerated  values,  the  factor  of  attenuation  V  depends  on 
parameter  q,  entering  the  limiting  conditions.  For  functioning  the 
l4e rtz  U  (vertical  polarization)  parameter  q  is  equal  to 


a  = 


(2.05; 


l^PTi  ’ 

where  t|  -  composite  dielectric  constant  of  medium. 
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For  functioning  the/-/ertz  W  (horizontal  polarization)  parameter 
q  is  equal  to 

q  =  im]  t)  —  1.  (2-06) 

In  the  latter  case  it  is  possible  to  virtually  assume/set  q=®,  since 
parameters  m  and  i)  are  great. 


Thus,  attenuation  factor  V  is  a  function  from  dimensionless 
quantities  x,  y,  y',  q: 

V  =  V(x  ]y,y'\q).  (2.07) 

Together  with  the  attenuation  factor  V  convenient  to  examine  function 
'F,  connected  with  it  through  which  V  is  expressed  according  to  the 

formula  n 

V  =  2  1  nx  e  TY.  (2-08) 


Function  T  satisfies  the  differential  equation 
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ffj 


where 


liL-ry^r^lT-O. 

Ol~  OX  U  ' J 


r  (y)  =  m-  (e  —  1). 


moreover  e  =  e  (A)  is  dielectric  constant  of  air  as  function  from  the 
height/altitude.  Equation  (2.09)  isobtained  by  transition  to  the 
dimensionless  quantities  from  the  equation 

_£I_  _  2ik  -  Jr  (2-  -  e  —  1  W  =4.  (2-11) 

Ok-  Pi  '  a  / 

in  which  coefficient  with  T  proportional  to  the  given  refractive 


index 


Coefficient  when  in  equation  (2.09)  it  is  convenient  to  designate 
by  one  letter  and  to  assume/set 


We  have 


p(y)  =  r  (y). 


p  (£/>  =  /fi*  (e  —  14  — 


so  that  p(y)  there  is  actually  the  same  given  refractive  index,  that 
only  expressed  through  the  dimensionless  height/altitude  y. 


During  the  use  of  designation  (2.13)  equation  (2.09)  will  be 


written 


ir  +  '-£-+p(y)v= o. 
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Besides  differential  equation  (2.15),  function  V  satisfies  the 
limiting  condition 


dy 


O 


-  qV  =  0  (npH  y  =  0). 


(2-16) 


Key:  (1).  with. 

With  x=0  it  has  a  special  feature/peculiarity  of  form. 


T  =  -z-yss-  ;  e 
2  f.-u  l 


• - -*»e  -  —  i.  (2.17, 


y  -  y’  —  2i<?.v  J ' 


In  Chapter  14  there  was  given  general/common/total  expression  foi 
function  Mr  in  the  form  of  contour  integral.  Integrand  in  it  is 
expressed  through  the  solutions  of  the  equation 


<r-i 


(2.18) 


where  t  -  composite  parameter.  (These  solutions  and  they  were  named 
above  high-altitude  factors). 


For  the  compositiovof  integrand  it  is  necessary  to  know  both 
solutions  of  equation  (2.18);  let  us  designate  them  through  fx(y,  t) 
and  f a (y,  t) . 


With  large  y  these  functions  have  the  asymptotic  expressions 
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A  to.  -> =  - 


c  t 


exp 


,  p  (l/.  - : 


i  \  ]  p  (//)  —  t  du 


A  to-  0  =  — 


exp 


—  i  I  I  p(u)  —  idu 


(2.19) 


(2.20) 


»  P  if)  —  ’■ 

Here  c',  c",  r  -  constants  whose  values  are  unessential,  since  they 
drop  out  from  the  expression  for  Y.  In  the  case  of  homogeneous 
atmosphere  when  p(y)-y,  functions  fj(y,  t)  and  f2(y ,  t)  are  led  to 
the  composite  Airy's  functions  wa{t-y)  and  ws(t-y),  moreover  then  it 
is  possible  t»  assume  c ' =c"  =1  and  r=t. 


Let  us  assume 


D«  «  =  h%- 


I: 


oh 

oy 


(2.21) 


On  the  strength  of  equation  (2.18)to  which  satisfy  ix  and  f2, 
this  value  will  not  depend  on  y. 
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Let  us  designate  values  and  with  y=0  through  f’a(0,  t) 

and  f'j(0,  t)  and  will  compose  the  function 

F(i.  y,  y',  q)  = 

i  !  lUO.  1) -=r  q/.iO.  t)  ) 

=  Iff  /•  n  I n  ~  f;,o.n-,/,,o.n  ’•  <*  "  i  ■  l-2-' 

Function  is  determined  with  y’>y  by  the  contour  integral 


y  =  4r  \*t:F(i'  y>  y '  du 


(2.23) 
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Undertaken  on  the  duct/contour,  which  covers  in  the  positive 
direction  all  of  the  pole  of  integrand.  As  shown  in  Chapter  14, 
function  T  satisfies  all  conditions  set  above  and  gives  the  solution 
of  our  problem. 

3.  Case  of  normal  refraction. 


The  case  of  normal  refraction  is  characterized  by  the  fact  that 
the  given  refractive  index  M(h)  is  the  monotonically  increasing 
function  from  height/altitude  h,  and  therefore  coefficient  p(y)  is 
that  monotonically  increasing  function  from  y.  In  this  case  it  is 
possible  to  approximately  express  f^y,  t)  and  f2(y,  t)  through  the 
composite  Airy's  functions  from  the  argumerif'£,  determined  by  the 
equalities 


I  I  pun —  <*<*«  =  -5- 

t 


\  \  t  —  P(ti)  dll 


(3.01) 


(3.02) 


where  b  is  a  root  of  the  equation 


P  (b)  =  i- 


(3.03) 


Near  y=b  resolution  of  value  $  according  to  degrees  of  y  -  b  it 


will  begin  from  the  linear  terms,  namely 


=  82001203 
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jf 


I  =  }71F)(b-y)  + 
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We  can  approximately  assume 


/.-v 


-  ■£■*© 


also,  with  the  same  approximation/approach 


l«r_  I  a?  ■ u  <*.•  I  d.v 


(I).  (3-06) 


whence 


D,.*  -2:. 


Replacing  here  y  on  y*  and  |  on  £’,  we  will  obtain  expressions 
for  f  i  ( y ’ ,  t)  and  fa(y’  t).  The  value  £,  which  corresponds  y=0,  we 
will  designate  through  £».  With  these  designations  vewill  obtain  fo: 
function  F,  determined  by  formula  (2.22),  the  following 
approximation: 

r~  i  tf  d7~  1/  dy  . 


r  •  1/  dy  1 '  ■ 

f=-|/  -t!' 


W.<5D)-*  1 

(■*•) 

e-j  (;pi  -  <7 1 

f  dy  \ 

U  J 

Being  substituted  into  formula  (2.23),  this  expression  can  be 
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used  for  calculating  the  field  both  in  the  shadow  zone  and  in  the 
illuminated  region.  In  the  shadow  zone  attenuation  factor  (and  also 
function  Y)  is  calculated  from  a  number  of  the  deductions,  which 
correspond  to  the  roots  of  the  denominator 

“klo)  ~  Q  (-§-)v  ti'i  (Ip)  =  0 .  (3.09) 


In  the  illuminated  region  function  Mf  is  calculated  directly 
with  th  aid  of  the  contour  integral,  moreover  the  main  section  of 
integration  will  lie/rest  near  the  real  negative  values  of  t.  But 
with  negative  t  of  value  £e,  £  and  £*  they  will  be  also  negative. 
Assuming  that  these  values  are  sufficiently  great,  it  is  possible  to 
replace  for  function  w,  and  w,  with  their  asymptotic  expressions 


_ ]_  i  •  k- 

4  <-l)  <e  3 


(3.10) 

(3.11) 


This  replacement  is  reduced  to  the  fact  that  for  functions  fj(y, 
t)  and  fj(y,  t)  are  utilized  asymptotic  expressions  (2.19)  and 
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As  a  result  for  function  F  is  obtained  [according  to  the  formula 
(2.22)]  the  following  expression: 


F  =  - 

*  n 


i  pu/i  —  <  J  p(y')  —  t 


exp 


1  j  l  p{u)  —  tdu 


L  y 


//  —  i  |  p  iOi  — 


Q  —  /  I  /■  <0j  —  1 


exp 


'  |l  P  (M)  —  /  du  —  :  j  )  p  (U)  ~  t  du 


(3.12» 


This  formula  presents  the  generalization  of  formula  (6.11)  of  Chapter 
14.  The  latter  is  obtained  from  (3.12)  after  replacement  of  p(y)  by 


zero. 


The  substitution  of  value  (3.12)  into  the  contour  integral  gives 
for  the  attenuation  factor  the  expression,  which  consists  of  two 
members  of  whom  the  first  corresponds  to  the  incident  wave,  and  the 
second  -  to  wave,  once  reflected  from  the  earth’s  surface  with 
Fresnel’s  coefficient.  The  incident  wave  presents  the  imposition  of 
waves  with  the  phase 
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but  the  wave  reflected  presents  the  imposition  of  waves  with  the 
phase 

i 

<5  U!  =  x:  -  j  )  p{u)-tdu-r  \  Vp(u)  —  t  da.  (3.14) 


These  expressions  correspond  to  geometric  optic/optics. 
Integrals  can  be  calculated  according  to  the  method  of  steady  state, 
moreover  the  phase  of  the  incident  wave  will  be  equal  to  outer  limit 
u(t) ,  and  the  phase  of  the  wave  reflected  -  to  outer  limit  e(t). 
Function  uit)  reaches  its  outer  limit  with  t,  determined  from  the 
equation 

0/  (U  =  x~4r\  -=$L=  =  o.  (3.15) 

t 


while  function  o(t)  -  with  t,  determined  from  the  equation 
<T'  (/)  =  x - 1  \  -A—.  -  4r  f  d“  =  =  0-  (3-16) 

-J  1  p(ui~!  2  }  1  pf.u)-i 

0  b 

From  the  point  of  view  of  geometric  optic/optics  the  horizon 
distance  is  determined  from  the  condition  so  that  this  point  could 

u;'»+k. 

reach  the  wave  reflected  a  Extreme  value  t,  at  which  this  still 

occurs,  exists  t=p(0).  This  value  must  be  at  the  same  time  the  root 
of  equation  (3.16). 
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Consequently,  between  >:,  y  and  y*  must  exist  the 
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relationship/ ratio 

r  =  _L\  du  -  --  ;--LV - in - (3.1?) 

-  |  ptut  —  plUi  ‘  -  |  r'“!  —  P<0> 

t  O 

which  gives  formula  for  the  horizon  distance  during  the  normal 
refraction. 


More  exact  expression  (3.08)  for  F  shows  that  at  value  of  t=p(C) 
to  apply  formula  (3.12)  already  is  impossible.  In  fact,  at  this  value 
of  t  valuer,  turns  into  zero,  and  to  use  formulas  (3.10)  and  (3.11) it 
goes  without  saying  inacmissibly .  Nevertheless,  it  is  possible  to 
consider  that  value  x,  determined  from  (3.17),  approximately  gives 
the  boundary,  which  is  determining  the  illuminated  region  where  is 
applicable  reflecting  formula,  from  the  shadow  region  where  let  us 
use  the  series/row  of  deductions.  As  the  in  other  words,  it  is 
possible  to  consider  that  the  amplitude  of  field  begins  rapidly  to 
decrease,  when  x,  growing/rising,  i  passed  through  value  (3.17).  In 
this  sense  it  is  possible  to  apply  term  "horizon  distance"  and  in  the 
diffractive  theory. 


4.  Asymptotic  integration  of  differential  equation  with  the 
coefficient,  which  has  the  minimum. 


In  the  presence  of  superrefraction  the  given  refractive  index 


K(h)  there  will  not  be  by  monotonic  function  from  the 
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height/altitude,  but  will  have  one  or  several  minimum;,  which 
correspond  to  separate  waveguide  channels.  We  will  examine  the  case 
of  one  minimum;  th  appropriate  height /altitude  we  will  call  the 
height/altitude  of  inversion  and  designate  through  A.-- 

Coefficient  p(y)  of  the  differential  equation 
■jpr  -  P  U'>  /  =  if 

is  proportional  to  M(h)  and  therefore  it  will  also  have  one  minimum 
at  value  y  =  Ff.  appropriate  h  —  /«; 

We  will  count  p(y)  analytic  function  from  y.  Equation  p{y)=t 
will  have  in  the  region  interesting  us  two  roots:  y*b,  and  y=bs. 

With  t  real  and  lying  between  p(0)  and  both  roots  will  be 

real;  at  other  values  o  t  the  roots  can  be  composite. 
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We  should  have  this  asymptotic  expression  for  functions  fi(y,t) 
and  f j(y,t)which  would  be  correctly  evenly  at  all  values  of  y  and  t 
in  question,  switching  on  value  >  —  P  iiu),  with  which  roots  b.  and  b- 
coinciae. 
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The  approximation  for  fj  and  f,  used  in  paragraph  2  through  the  j 
Airy's  functions  here  narrower  is  not  applicable.  Its  validity  was 
based  on  what  with  the  aid  of  substitution  (3.01)— (3.02) ,  that 
defines  £  as  holomorphic  functio  from  y,  equation  (4.01) 
approximately  was  reduced  to  the  equation 

-~fr  __ |a.  =  0,  (4-02) 

in  which  coefficient  with  the  unknown  function  it  had  the  same 
monotonic  character  as  in  the  initial  equation,  l^ov  we  should  take  as 
the  standard  equation  no  longer  equation  (4.02)  for  the  Airy's 
functions,  but  the  equation 

-2- (-T  ( 4-03, 

for  the  functions  of  parabolic  cylinder,  since  this  is  the  simple 
equation  in  which  the  coefficient  with  the  unknown  function  has  the 
same  character  (with  one  minimum)and  coefficient  p(y).  Th 
substitution,  which  connects  (,  and  y,  must  be  selected  so  that  the 
value  p(y)-t  would  become  zero  simultaneously  with  value  -j-  —  v  and  so 

that  at  the  high  values  of  these  values  would  be  obtained  correct 
asymptotic  expressions.  These  conditions  satisfies  the  substitution 

j  1  p  </j)  —  i  dy  —  —  I  ]  Is  —  4v  di,  (4.0ii 

when  parameter  v  is  selected  so  that  it  would  be 
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_  -■]' _ 

|  I  piy\  —  idu=  ~  I  t  Is  -4\  d'T.  (4.05) 

1 1  *  -2;  I  % 

Integral  in  right  side  (4.05)  is  equal  to 

‘.7  \‘\ 

4-  f  !  5*  -  4v  dZ  =  uiv.  (4.06)  ^ 

i  v 

Therefore  equation  (4.05)  can  be  written  in  the  form 

i; iv  =  [  \  p(y)  —  t  dy.  (4-07) 

ii 
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It  gives  v  as  function  cr  t.  This  function  will  be  holomorpnic 
near  t  --  p  u/.),  namely  we  will  have 


Y 

_  /’(hi  —  * 

(4.08) 

Assuming/sett ing 

l'2rr  (i/.i 

S  =  |  1  p  (u)  —  t  du. 

(* 

(4  uy> 

f, 

So  ~  —  f 

tl 

1  p  f  dij  — T  |  |  p  (in  —  t  dy. 

<4.10* 

0  f* 


we  can  write  substitution  (4.04)  in  the  form 

* 

s-s0  =  4-j  IF-Tva;.  H-ii) 

o 

the  right  side  of  this  expression  is  equal  to 


—  II  dz  -  -r  Z  \  Z2  -  4v  - 


v  ig(;  —  i 


4v)  -  4r  1?  4v. 


(4.12; 
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Hence  we  can  conclude  that  with  $>0  value  S  —  So  —  -r lg vwill  be 
holomorphic  function  from  v  near  v  =  0.  ana  with  $<0  holomorphic 
will  be  values  S  — S0 — 1-lgv  and  S.  However,  since  with  y=0  (on  the 
earth's  surface)  knowingly  $<0,  then  we  deal  concerning  £he  second 
case,  and  then  holomorphic  functionfrom  v  will  be  sum  SD~4-lgv. 
This  observation  will  be  necessary  by  us  subsequently. 


In  the  asymptotic  approximation  of  solution  of  equations  (4.01) 
in  question  and  (4.03)  they  are  connected  with  the  relationship/ratio 


(4.13) 


The  solutions  of  equation  (4.03)  are  the  functions  which  are 
expressed  as  the  functions  of  parabolic  cylinder  Dn  (z).  satisfying 
the  equation 

-r  ( n  +  -f  -  4  2’)  °n  W  =  0. 


(4.14) 
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Functions  D„ ( z )  are  well  investigated.  We  will  not  enumerate  their 
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properties,  but  let  us  refer  to  the  book  of  Whittaker  and  Watson  the 
"bourse  of  contemporary  analysis”,  Vol.  2,  GITTL,  1934,  where  are 
given  the  principal  formulas.  As  determination  Dn  ^z)  it  is  possible 
to  take  the  series/rov 


D.  W  -  f  S' - r4 1  {)■■■ 2":  <4-*5> 

Bind 

Equation  (4.03)  is  obtained  from  (4.14)  by  replacement  z  on 


-i  ■ 


£e  *  and  n+i/2  on  1  !V-  The  solutions  of  equation  (4.03)  will  be  the 
functions 


g.Jl)-  D  i  (e  '  4  :)  ■ 


-i\  - 


.(.*+;)• 


(4.16i 

(4.17) 


With  real  ones  v  and  $  values  g,($)  and  g*($)  will  be  composite 
conjugated/combined. 


From  properties  D,,  (z)  it  escape/ensue 

Si  (-;.)  =  e~vn  ^  s,  Q  +  —4-——  e“  ^  ^  g,  (;).  (4 ■  18) 


(t-/v) 


and  also 


£.,(-;»  =  e' 


V«  J}_ 

8*0 - /l"'  "v  e  «  '  «  g1  (;).  (4.19) 


(t  +  IV) 


For  us  are  essential  asymptotic  expressions  for  g A ( $ )  and  g:($). 


In  the  region,  which  adjoins  the  positive  real  axis,  we  have 
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Using  formula  (4 
Si 
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Latter/last  expression  correctly  also  with  large  ones  v. 
Asymptotic  expression  for  g;($)  is  obtained  hence  by  replacement  of  i 
on  -i. 

In  order  to  obtain  the  formula,  valid  near  the  negative  real 
axis,  we  they  must  use  relationship/ratio  (4.18). 


1  - 


iv®  —  2v  —  i"  — 

4 


2? 


(4.20) 


.12),  we  can  also  write 


Jiv  n  v  .  .  v  . 

Q  =  e— ' 


x 


i-  --  4v 


•exp 


■f  [IT 


(4.21) 


We  will  have 
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p4 


3v« 

a  iv  , 

£,C)  =  e  4 

« 

h  e  ;  - 

» 

- — - exp 

,  -  4v 

-  4-  J  i  ;s  -  4v  d±  ■ 

i 

»'E  r 

va  .  .  n  v  i\ 

-‘Tp  2  — r ig x 

X  . . exp  ~  I  1  ;s  -  4 V  <f£  .  (4.22) 

i  ;2-4v  l  ‘  i  .1 

Nov  we  can  construct  the  solution  of  equation  (4.01),  which 
satisfies  all  stated  requirements. 

Let  us  assume 

ci  (v)  =  e‘  8  rlR'-s‘l.  (4.23, 

In  view  of  the  property  of  value  S„  noted  above,  index  in  (4.23) 
is  holomorphic  function  from  v  •  also  near  v  -  o. 

The  proper  solution  of  equation  for  the  high-altitude  factor 
will  be  the  function 


Ate.  /)  =  f1(v)]/2-^ft(;). 


Higher  than  the  layer  of  inversion  (with  S-SB>>1)  this  function  has 
the  asymptotic  expression 


.  n 

/»0r.  0  =  ~.e  =e' 

i  p  (.w)  —  ‘ 


iS-VS, 


k-k.v  5.. 


ij  ^n'iiiffiir  ■  iij.'  i  r,  ~iiini  y  Hr  J-  c ..  i  ^  >7i » ■ ,  ,“i ■  ri  .  „  '■(  ii  'ii;  i  lrfi  i'ni  -j  ii  iirSt 
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ensuing/escaping/f lowin  out  from  (4.21). 
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Lower  than  the  layer  of  inversion  (with  S0-S>>1)  asymptotic 
expression  for  ix{y,t)  will  have  the  form 


h  (y-  t)  =  Xi  (v)  tt: 


i  p(y)  —  t 


re 


i  p\y)  —  t 


:e~ls  ;(4.2C) 


where  we  placed 


I  2.i 

7.i  hj  =  — -7 - -e 


It-) 


V.T 

-  (J  If  V) 


(4.2,  , 


With  the  aid  of  the  known  asymptotic  expression  for  function 
t(\t  —  iv)  it  is  easy  to  show  that  at  high  positive  values  Y  function 
y.i  approaches  unity.  Since  when  v.  l  the  second  term  in  (4.26) 
becomes  small  in  comparison  with  the  first,  both  expressions  for 
fj(y,t)  will  then  in  form  coincide  It  is  substantial,  however,  then 
that  our  expressions  for  fi(y,t)  are  valid  not  only  with  the  large 
ones,  but  also  at  low  values  v'  up  to  v  «  0.  and  that  they 
represent  holomorphic  functions  from  v  near  v  ~  u- 


The  corresponding  expressions  for  f.(y,t)  will  be  obtained  from 


preceding/previous  by  replacement  of  i  on  -i.  In  order  to  extract 
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them  clearly,  let  us  assume 


a-.  _  _ \_  j  5  , 

c„  (vi  =  t  '  e 


1  2- 

7.i  ,v»  - - : - t 


/  <  —  .'I 


Then  it  vill  be 


(4.28> 

(4.29i 


}:(y.  :i  =•  c.cviV  2-^- £„(;). 

f  *■. 


(4.30) 


and  asymptotic  expressions  for  f2(y,t)  will  take  the  following  form 


with  S-SC>>1 


/..  <  a.  ?  i  =  -- 


(4.31) 


i  r 


with  S,-S>>1 


/.(«.  fl  = 


— :..s0 


;> « ;/ 1  —  ; 


—  e 


-v.i  t 


-  e‘s . 


(4.32) 


l  V  (</»  -  : 


Thus,  the  problem  of  the  asymptotic  integration  of  equation  fo 
the  high-altitude  factors  we  have  solved. 
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5.  Investigation  of  attenuation  factor 
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To  us  is  to  now  substitute  those  found  for  f2(y,t)  and  f2(y,t) 
of  expression  in  formula  (2.22)  for  F  and  to  investigate  attenuation 
factor  V  or  connected  with  it  function  For  simplicity  we  will  be 

restricted  to  case  <?  =  pc.  which  corresponds  to  horizontal  polarization 
In  this  case  function  F  is  reduced  to  the  form 

F(t.  y,  y',  oc)  =  -5-~  f,(|T.  t)  \h(y,  0~t^oT /.(!/•  0|  •  (5-01) 

For  functions  (4.24)  and  (4.30)  the  Wronskian  D22  is  equal  to  the 
constant  value 


D1;  =  -2  i. 


which  more  easily  in  all  is  derived/concluded  from  asymptotic 
expressions  (4.25)  and  (4.21).  We  will  asstime  that  so  that 

S(y’)-S,»l,  and  let  us  examine  two  cases:  when  the  second 
height/altitude  is  also  higher  and  when  it  is  lower  than  the  layer  of 
inversion.  In  the  first  case  we  will  consider  S(y)-S„>>l  that  makes 
it  possible  to  us  expressions  (4.25)  and  (4.31)  for  f2  and  f,.  In  the 
second  case  we  will  count  S„-S(y)»l  and  use  expressions  (4.26)  and 
(4.32) . 


In  the  first  case  we  will  have 


JS 


‘  ‘r  >  — -  ?  i  —  ; 


,  ..  _ V-- 


— nv  r—Z:S,. 

f  -  :>:,t 
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The  individual  members  of  this  expression  allow/assume 
interpretation  on  the  basis  of  geometric  optic/optics.  It  is  obvious 
that  the  wave,  which  goes  downward,  must  have  phase  factor  e-,s,  and 
the  wave,  which  goes  from  bottom  to  top  -  phase  factor  e;-s\ 

Expression  (5.03)  shows  that  there  is  only  one  wave,  whic  goes 
downward,  namely  incident  wave  with  the  complete  phase 

i-\  (?)  =  x :  —  S  if/*)  —  S  <«)  i5  04 1 

[  they  added  here  member  xt  from  the  exponential  factor  in  the 
integral  (2.23)].  This  phase  coincides  with  phase  (3.13)  of  the  case 
of  normal  refraction,  which  is  natural,  since  this  wave  yet  did  not 
reach  the  layer  of  inversion. 

As  far  as  the  waves,  which  go  from  bottom  to  top,  are  concerned, 
them  there  will  be  countless  set;  these  waves  will  be  obtained  by  the 
resolution  of  second  term  (5.03)  in  the  series/ro  according  to 
degrees  e_nv-  They  will  correspond  to  the  waves,  repeatedly  reflected 
from  the  earth’s  surface  and  from  the  layer  of  inversion. 

Page  330. 

The  phase  of  the  wave,  once  reflected  from  the  earth's  surface,  will 
be  equal  to 

<p  (0  =  xt  -r  S  iy')  -  S  («)  -  arc . 


(5.05) 


. . . . .  . . . I . in 
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This  expression  differs  from  12.14)  in  terms  of  the  latter/last 
term  which  cannot  be  obtained  from  the  geometric  optic/optics.  This 
term  is  equal  to 

r(4r-< .) 

arc  —  arc - 7 - —  2v  lev  —  2v.  (d.06) 

With  large  positive  ones  v  it  becomes  zero,  but  with  small  ones 
v  it  play  important  role,  since  because  of  it  entire/all  phase  ®(t) 
remains  hoiomorphic  function  from  v  near  v  =  0,  in  other  words, 
near 

Let  us  examine  now  the  case  when  point  y  is  located  below  the 
layer  of  inversion,  moreover  S*-S>>1. 

Using  expressions  (4.26)  and  (4.32)  and  utilizing  the  equality 


7.1  <v)  Xi  (v>  —  e  'nv  =  1 


we  obtain  after  some  linings/calculations 


js  <*••)  -■as. 


sin  S  (y) 


»  p  <y')  —  t  i  ~p  (y)  —  1  *»e  ’jS*  —  ,e  ”v 


In  this  case  there  is  not  by  one,  but  a  countless  quantity  of 
waves,  which  go  downward,  since  to  the  incident  wave  ar  connected  the 


1 
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waves,  reflected  from  the  layer  of  inversion  as  from  upper  boundary. 
Furthermore,  is  an  infinite  quantity  of  waves,  reflected  from  the 
earth/ground  and  which  go  from  bottom  to  top.  All  these  waves  are 
formally  obtained  by  expansion  (5.08)  into  the  geometric  progression 
according  to  degrees  e-3V- 

The  complete  phase  of  the  wave,  which  did  not  undergo  terrain 
echo,  is  equal  to 

co  (/)  X!  -  S  (»/')  —  S  (ij)  —  arc  (5.09) 

or 

to  (/) «  xt  -  S  (;/)  -  S  (v)  -4-arc-J-,  <5- l0) 

while  the  complete  phase  of  the  once  reflected  wave  is  equal  to 

<i  (;)  =  JC l  -  S  (in  -  S(1>)  -r  4-  arc-—-  (51  >) 


Pag  331. 

Expression  for  «(t)  does  not  coincide  with  (3.13)  or  (5.04),  which  is 
natural,  since  the  incident  wave  traversed  the  layer  of  inversion. 
Expression  (5.11)  differs  from  (5.05)  in  terms  of  the  fact  that  the 
additional  term  enters  with  factor  1/2. 
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Me  spoke,  until  now,  about  the  phases  of  different  members  of 
integrand.  To  each  such  term  corresponds  in  the  attenuation  factor 
integral  on  t.  If  we  calculate  these  integrals  according  to  the 
method  of  steady  state,  then  each  of  the  gives  in  the  attenuation 
factor  the  term,  which  presents  wave  with  the  phase,  equal  to  the 
outer  limit  of  the  phase  of  integrand. 

Of  course  this  method  of  calculating  the  attenuation  factor  is 
applicable  only  in  the  illuminated  region,  whereas  in  the  region  of 
shadow  it  is  necessary  to  use  the  series/rov  of  deductions. 

6.  Formula  for  the  horizon  distance. 

During  the  normal  refraction  (paragraph  2)  we  estimated  distance 
of  the  norizon/Ievei  as  this  value  of  the  horizonta  distance  of  >: 
which  gives  the  boundary  between  the  region  of  the  applicability  of 
reflecting  formula  and  the  region  of  the  applicability  of  the 
series/row  of  deductions.  At  this  value  of  x  the  extremum  of  the 
phase  of  the  wave  reflected  falls  at  extreme  value  of  t  at  which 
phase  itself  is  still  real. 

In  the  presence  of  superrefraction  are  many  waves  reflected.  But 
we  can  expect  that  main  role  plays  waves,  once  reflected  from  the 
earth’s  surface.  Since  the  "horizon  distance  is  not  the  strictly 
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defined  concept,  we  are  right  to  refine  it,  understanding  under  it 
horizon  distance  for  the  once  reflected  wave. 


The  phases  of  the  once  reflected  wave  have  found  we  in  paragraph 
5.  According  to  (5.05)  and  (5.11)  we  have  when  ?/'  > y  ;/. 

t'  i 

q  (?)  =  xi  —  [  |  ptui  —  i  du  —  ]  V*p<a‘‘  -tdu-  arc—  (6.05 1 

(•  £ 

and  when  y">y-.  y<y.. 

_ ? _ . 

q  (? I  =  x;  —  j  i  piu)  —  :du—  j  |  p  w  —  du  — rr  zr c  .  <6 

£  i 


These  formulas  car.  be  combined,  after  assuming 


S*  iu.  i)  =  i  1  p  (tit  —  i  du  —  rr  arc 


X; 

>.i 


.')  —  i  I  P  tu)  — ;  du 


(J/  >  Uft. 

«6.0>» 

iy  <!',>■ 

(0.04  • 

Page  332. 

Then  it  will  be  both  when  y  >  yit  and  when  tj  <  y- 

q  (/)  =  xi  ~  S*  (u',  t)  —  S*  (u,  /).  (6.05) 

Let  us  note  that  S*  is  holomorphic  function  from  t  near  t  —  p  (w,). 

By  discussin  as  in  paragraph  2,  we  will  obtain  for  the  horizon 
distance  the  expression 
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■  _  _  r  <*'  ■  0  dS*  (“■  1 1  1  _  (6  r> 5| 

A  —  i.  or  1  at  j  (=p  (»f  I  ‘ 

Let  us  write  this  expression  in  the  more  explicit  form.  According  to 
(4.08)  near  t  =  p  (y,)  we  have 

v  =  jjM±z±  .  (6.07) 

V  2;*  <;■-> 

On  the  other  hand,  near  v  =  o 


».rc 


r  (t  - iv) 


=  (C  -  21g2jv  — 


(.6.08) 


and,  therefore, 


—  arc  ~  =  v  (C  —  1  —  lg  4v) 


(6.09) 


where  C=0.577...  -  Euler's  constant.  Therefore  when  y  >  Hi 


a  S' 
at 


-t\ 


du 


2  .  V  P(u)  ~ 


(C-lg4v).  (5.10) 


This  expression  has  a  limit  when /  —  p  ({/,),  v  —  0.  when  y  <  y,  latter/last 
term  is  absent,  and  into  the  integral  it  is  possible  to  directly 
substitute  value  Therefore  during  y<yt  Budde's  lifetime. 


as*  _  j_  :• 
—  a  * 


du 


at 


V p  (u)  —  p  [yi) 


(6. i 1 ) 


The  presence  of  the  second  term  in  formula  ($.02)  causes  the 
dependence  of  horizon  distance  on  the  wavelength.  In  order  to  come  to 
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1 ight/detect/expose  this  dependence,  let  us  return  from  reduced 
coordinates  x,  y  to  the 'usual  coordinates  s,  h,  where  s  -  horizontal 
distance  and  h  -  height/altitude.  Designating  through  y(  h)  the  given 
refractive  index  without  factor  10‘,  we  will  have 

p  (y)  =  2msn  (h),  (6.12) 

where  m  is  value  (2.04). 

Page  333. 

Instead  of  t  let  us  introduce  the  paramet err,  connected  with  t 
with  the  relationship/ratio 

i  —  2m-\.  (6.13) 

Then 

y  h 

f  1  dj  -  k\\  2u  (h)  —  2t  dh.  (6.1-h 

6  £ 

xt  =  ksx. 

Value  v  will  be  now  approximately  equal  to 

v  —  -_L,^lu(fe,)-T).  (6.16) 

1' 

Formula  for  the  horizon  distance  will  be  obtained  from  the  condition 
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Composing  a  difference  in  expressions  .  0),  we  obtain  when  h  >  hi 

s„  —  s,  =  — t=L=^- lg  41  =  — -le^2-  (6.21i 

IV  (hj  8  *>  )V(- ,  8  **  1  ' 

and  when  h  <  h( 


Sa  ~ ~  Si  — 


1 


—  - In  it 

lV  (*d  s  *> 


)V  (hi) 


lgit.  (6.22) 


This  difference  depends,  besides  the  relation  of  wavelengths, 
only  on  the  course  of  the  given  refractive  index  near  its  minimum 


Let  us  apply  our  general  formulas  to  the  case  when  the  given 
refractive  index  u( h)  depends  on  height/altitude  according  to  the 
hyperbolic  law 

H(A)  =  fi(A,-)^4-  iTCTL*  ..  (6-23» 

where  a  -  radius  of  terrestial  globe?  parameter. 


In  this  case 


ft*  (A/» 


2 

u  (h,  -  1 1 


(6.24) 


The  integrals,  entering  in  ©(t),  will  be  elliptical,  but  when 
t  =  P  (/>.)  they  are  calculated  elementarily,  and  we  obtain  for  F(h)  the 
following  expressions.  When  h<rh,’ 
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F  (h)~  —  1  2a  (h  —l)  -  }  2al  + 

_  i  nm  HE  £5ZLzLigZ2  _  ig  (6-25) 

I  2  I  ®  )Oi.  —  /  —  lA/>  -r  /  &Vhi~l—  y 1  j 


and  when  *  > 


where 


F(A)  =  I  2a  (h  -/)  -r]  2al  — 


_]  '  £  <*.  ~  '•>  ilp  .  lr  K'ft,--/  -  VT| 

\  2  ,  -rXTT  K  \rhITT-  iTi 


*>-1 


■  As, 
(6.26) 


fa(h-  -  I l  („  •  1  2k"-(h,-l)3  i 

- 2 - ICi - 2“  lg - a  |* 


(6.27) 


Here 


C1  =  7  1g2  —  4-f-C  =  1,429. 


(6.28) 


For  the  comparison  let  us  not  that  the  horizon  distance  in  the 
absence  of  refraction  is  equal,  as  is  known, 


s'  =  ) r2ahr  -  I  *2 ah. 


(6.29) 


Thus,  an  increase  in  the  horizon  distance  as  a  result  of  the 
refraction  is  equal 

s _  s'  =  \F  (ti)  —  I  "2aFj  t  \F(h)  —  }  r2ah).  (6.30) 
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In  all  preceding/previous  reasonings  we  assumed  that 
heights/altitudes  h  and  h'  were  small  in  comparison  with  a  radius  of 
Earth  a.  But  the  preceding/previous  formulas  are  applicable  also  to 
the  case  of  the  wave,  which  goes  from  infinity  (for  example,  from  the 
Sun).  A  difference  F(h’)  —  )'2ah'  has  when  k' —  o c  the  final  limit,  equal 
to 


lim  \F  (A*)  -  |  2ah' }  = 

h  '-»»o 


=  12 


a  ,g  > 


VT 


1  ft,  -  /  -  V  ! 


—  As. 


(6.31) 


Replacing  in  (6.30)  the  first  two  members  by  their  limiting 
value,  we  will  obtain  for  an  increase  in-  the  horizon  distance  the 
following  expressions: 


when  h  <  h; 


s  —  s'  =  2  1  2al  —  1  '2a  (h  -  /)  —  I  2ah  -  As 


1  fa  (ft  -  h  |,n  1  'hi- l  -  \r>>~  ‘  _  olr  LliSJ-ZJLL  I  ifi.3'21 

-]  —2“  lg ',gr— -l'i* 


and  when  h  >  /i, 


s-s'  =  21  2ai  -  1  2a (/!—/)  —  1  2aA  -  2As 


O  (ft;  *r  /(  |,„  I'ft-l 
2  i  8  i/-rr 


j  8  jOTTT  -  fft^T 


2ig  Ui-  _LL  I 
g  -7  -  r'/  ) 


(6.33) 


To  this  increase  in  the  distance  corresponds  "lead  angle" 

t>=.LrjL.  ^6.34) 

C 
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Since  present  theory  does  not  consider  refraction  in  the  upper  levels 
of  atmosphere,  then  for  the  comparison  with  the  observed  lead  angle 
it  is  necessary  to  value  (6.34)  to  add  the  value  o  normal  refraction 
at  the  horizon/level. 
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Chapter  16. 

The  radiowave  propagation  near  the  horizon/level  during  the 
superrefraction  1 . 

FOOTNOTE  1 .  Foch,  Weinstein  Belkin,  1956.  ENDFOOTNOTE. 

In  the  beginning  of  chapter  is  more  precisely  formulated  the 
concept  of  the  horizon/level  in  the  presence  of  surface  tropospheric 
waveguide.  Further  are  investigated  approximations  (of  type  of 
reflecting  formula)  for  the  attenuation  factor.  Under  the  assumption 
of  the  hyperboli  law  of  a  change  of  the  refractive  index  with  the 
height/altitude  are  derived  the  formulas  for  the  horizon  distance  of 
the  straight/direct  and  reflected  waves.  In  second  half  of  chapter 
are  given  numerical  results  for  several  typical  examples,  moreover  it 
is  assumed  that  the  transmitting  antenna  is  arranged/located  highly 
above  the  layer  of  inversion,  and  receiving  -  within  the  layer  of 
inversion  on  the  small  height/altitude.  The  obtained  results  give  the 
estimation  of  the  possible  values  of  attenuation  factor  at  the 
horizons/levels  and  is  shown  its  dependence  on  the  distance  and  on 
the  wavelength.  The  conducted  investigation  shows  the  advisability  of 
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the  introduction  of  the  concept  of  the  horizons/levels  during  the 
analysis  of  hyperdistant  propagation. 

1.  Introduction. 

The  theory  of  radiowave  propagation  above  the  spherical  earth's 
surface  in  the  presence  of  the  unhomogeneous  atmosphere  the 
refractive  index  of  which  depends  only  on  height/altitude  was 
developed  in  Chapters  14  and  15.  In  Chapter  15  is  given  investigation 
of  attenuation  factor  in  the  unhomogeneous  atmosphere  near  the 
horizon,  moreover  the  concept  of  the  horizon/level  is  determined  for 
the  heterogeneous  (laminar)  atmosphere  of  any  type.  The  determination 
of  the  horizon/level  introduced  there  coincides,  in  the  case  of 
unhomogeneous  atmosphere  without  the  Aversion,  with  the 
determination  of  the  boundary  of  shadow,  which  escape/ensues  from  the 
laws  of  geometric.- optic/optics.  But  if  is  an  inversion  of  the  given 
refractive  index,  then  the  horizon/level  is  necessary  to  find  of  the 
thinner  wave  considerations:  its  position  in  this  case  depends  also 
on  wavelength. 
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Assuming  that  during  the  removal/distance  beyond  the 


i 

* 

* 


horizon/level  the  attenuation  factor  rapidly  decreases,  then  it  is 
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possible  to  conditionally  consider  (as  this  is  made  in  Bhapter  15) 
that  the  horizon  distance  estimates  distance  radiowave  propagation. 

In  this  wa  is  obtained  simple  formula  for  the  superstandard  range  of 
radio  waves  during  the  superrefraction;  into  this  formula  enter  the 
height  of  receiving  and  transmitting  antennas,  the  wavelength  and  the 
parameters,  which  characterize  the  course  of  a  change  in  the  given 
refractive  index  with  the  height/altitude  (M  -  profile/airfoil). 
Especially  simple  form  takes  the  formula  of  distance  for  the  given 
refractive  index,  which  depends  on  the  height/altitude  according  to 
the  hyperbolic  lav  (paragraph  5  6hapter  15). 

The  analysis  of  hyperaistant  propagation,  produced  in  Chapter  15 
on  the  basis  of  the  concept  of  the  horizon/level ,  requires,  however, 
some  refinements.  It  is  first  of  all  desirable  to  explain,  what 
values  accepts  attenuation  factor  at  the  horizon/level  and  as  depends 
attenuation  factor  near  the  horizon/level  on  distance,  the  wavelength 
and  parameters  of  the  layer  of  inversion  (height/altitude  of  this 
layer,  its  average/mean  gradient,  etc.).  For  this  it  is  obvious  it  is 
necessary  to  calculate  attenuation  factor  in  some  special  cases, 
since  in  general  form  this  problem  to  solution  does  not  yield.  If  we 
in  this  case  will  explain  how  rapidly  attenuation  factor  decreases  in 
the  shadow  region  (beyond  the  horizon/level)  and  as  rapidly  it 
grows/rises  to  the  values  of  the  order  of  one  during  the 
removal/distance  from  the  horizon/level  into  the  illuminated  region, 
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then  ve  thereby  is  checked,  in  vhat  measure  the  horizon/level 
determines  in  the  practi  cal  cases  the  superstandard  range  of  radio 
waves . 

In  view  of  the  large  labor  expense  for  the  calculations  of 
attenuation  factor  during  the  superrefraction  it  is  possible  to 
conduct  calculations  only  for  a  small  number  of  typical  cases.  No 
comprehensive  calculations  as  during  the  normal  radiowave 
propagation,  cannot  be  here  fulfilled.  Therefore  we  were  restricted 
to  the  calculation  of  attenuation  factor  as  the  functions  of 
dimensionless  coordinate  $  in  four  cases,  which  allow  for  the. 
fixed/recorded  M-curved  and  at  the  fixed/recorded  heights/aitituaes 
of  the  corresponding  points  to  construct  the  dependence  of 
attenuation  factor  on  the  horizontal  distance  between  these  points 
for  four  wavelengths,  which  relate  as  by  1:3:9:27  (cf.  paragraph  7). 

In  this  way  it  proves  to  be  possible  to  refine  the  sense  of  the 
concepts  of  the  horizon  distance  and  superstandard  range  and  to 
answer  the  series/row  of  the  questions  presented  above,  in  particula 
by  a  question  about  the  dependence  of  the  phenomenon  of  hyperaistant 
propagation  on  the  wavelength. 

Let  us  recall  that  the  analysis  of  anomalous  propagation,  given 
in  Chapter  15,  let  us  use  if  and  only  if  one  of  the  corresponding 
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points  is  located  above  the  lowest  layer  of  the  atmosphere  of 
inversion,  another  point  can  be  both  within  this  layer  and  above  it. 
Therefore  during  calculations  of  attenuation  factor  we  were 
restricted  to  the  case  when  one  point  is  located  highly  above  the 
layer  of  inversion  and  another  -  within  the  layer  at  the  height  equal 
to  one  fifth  spot  height  of  inversion. 

Page  338. 

2.  On  the  concept  of  the  horizon/levei  in  the  presence  of  surface 
tropospheric  waveguide. 

Let  us  examine  in  more  detail  the  concept  of  the  horizon/level 
in  the  presence  of  surface  waveguide  (layer  of  inversion). 

Let  us  recall,  first  of  all,  the  beam/radiation  treatment  of 
normal  and  anomalous  propagation. 

For  homogeneous  atmosphere  the  given  refractive  index  is  a 
linear  function  of  height/altitude,  on  plane  s,h(s  -  distance  on  the 
earth/ground,  h  -  height/altitude),  the  rays/beams,  which  emerge  from 
source  Q,  take  the  form  of  the  curves,  turned  by  convexity  to  axis 
s(Fig.  1A).  Horizon/level  00’  is  determined  by  ray /beam  Q00' ,  which 
concerns  the  earth's  surface  at  point  0.  To  the  right  of  line  of 
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horizon  00'  is  located  shadow  zone,  where  field  penetrates  only  as  a 
result  of  diffraction,  to  the  left  -  illuminated  region.  For 
observation  points  which  are  located  in  the  illuminated  region  (it  is 
more  left  horizon/level  00'),  is  approximately  applicable  the 
reflecting  formula,  according  to  which  the  field  is  obtained  as  a 
result  of  the  interference  of  the  straight/direct  ray/beam  QP  with 
earth-reflected  wave  QP’P. 

Ravs/beams  ZhT  from  source  Q,  arranged/located  within  surface 
radio  duct  of  height/altitude  h.  (Fig.  lb),  have  within  the  waveguide 
a  convexity  upward  (from  axis  s),  and  higher  than  the  waveguide  - 
convexity  dow  (as  in  Fig.  1A) .  Therefore  ray/beam  Q1  is  passed  into 
the  space  above  the  waveguide,  and  ray/beam  Q2  proves  to  be  "closed" 
within  the  waveguide.  These  two  kinds  of  rays/beams  are  divided  with 
the  maximum  ray/beam  QO,  which  asymptotically  approaches  when 
s  —  oo  height/altitude  h  =  h,. 
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Fig.  1.  Illustration  to  concept  of  horizon/level :  a)  curing  normal 
refraction;  b)  during  superrefraction  -  on  geometric  optic/optics;  c) 
during  superrefraction  on  physical  optics. 
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Beside  the  straight/direct  rays/beams  into  the  space  above  the  layer 
of  inversion  fall  earth-reflected  wave,  for  example  Q1  *1',  separated 
from  the  closed  rays/beams  by  another  maximum  ray/beam  QO”  0’ , 
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asymptotically  approaching  height/altitude  hi  after  single  teirain 
echo.  Closed  prove  to  be  all  rays/beams,  which  emerge  from  the  source 
within  angle  of  0Q0" ,  formed  by  both  maximum  rays/beams. 

The  laws  of  geometric  optic/optics  in  this  example  lead  to  the 
conclusion  about  the  absenc  of  the  horizon/level  both  within  the 
waveguide  and  above.  Actually/really ,  through  observation  points, 
arranged/located  of  the  waveguide  more  to  the  right  of  rays/beams  1 
and  1',  pass  the  straight/direct  rays/beams,  which  emerge  from  Q 
within  the  angle  1Q0,  and  the  reflected  beams,  which  emerge  within 
angle  of  1  "QO" .  They  thread  entire  space  above  the  waveguide  more  to 
the  right  of  rays/beams  1  and  1',  and  therefore  geometrical  shadow 
and,  consequently,  also  the  horizon/level,  they  art-absent. 

It  is  easy  to  see,  however,  that  to  the  maximum  rays/beams  QO 
and  QO"  O'  and  to  the  rays/beams,  close  to  the  maximum  ones,  the  laws 
of  geometric  optic/optics  are  not  applied.  From  preceding/previous  it 
is  clear  that  precisely  these  almost  maximum  rays/beams  would 
transfer  (according  to  the  laws  of  geometric  optic/optics) 
electromagnetic  energy  to  the  large  distances  above  the  waveguide. 
Hence  it  follows  that  for  solving  the  ques  ion  about  the 
horizon/level  and  superstandairj  range  during  the  superrefraction  it  is 
necessary  to  draw  wave  considerations. 
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This  was  made  in  Chapter  15,  where  it  is  shown  that  in  the  space 
above  the  waveguide  is  certain  boundary  O’O'  (Fig.  1c),  more  to  the 
righc  by  which  earth-reflected  wave  penetrate  cannot.  This  boundary 
O'O’  is  the  hcrizon/level  when  the  layer  of  inversion  is  present, 
since  more  to  the  right  this  boundary,  i.e.,  into  the  shadow  region, 
the  field  (as  in  Fig.  1#^)  it  can  penetrate  only  by  diffraction. 

Together  with  boundary  of  O' O'  is  another  boundary  00,  more  to 
the  right  by  which  cannot  penetrate  the  straight/direct  rays/beams, 
which  did  not  experience  terrain  echo.  Boundary  O'O'  is  located  more 
to  the  right  boundary  00,  since  with  the  reflection  from  earth 
ray/beam  proves  to  be  more  to  the  right  parallel  to  it 
straight/direct  ray/beam  (cf.  rays/beams  Q1  and  Qln  1'  in  Fig.  lb). 

In  band  OO-O'-O'  straight/direct  rays/beams  do  not  pass;  therefore 
complete  field  to  this  band  to  beam/radiation  treatment  is  not 
subordinated.  More  left  boundary  of  00  compl'te  electromagnetic  field 
is  obtained  by  the  imposition  of  the  straight/direct  and  reflected 
beam. 

In  view  of  this  value  of  boundary  of  00  -  limit  of  the 
applicability  of  reflecting  formula  -  for  it  it  is  expedient  to 
introduce  the  special  name:  we  it  will  call  the  horizon/level  of 
direct  waves.  In  contrast  to  it  boundary  O'O'  we  will  name/cali  the 
horizon/level  of  the  waves  reflected.  While  durin  the  normal 
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propagation  these  "norizons/levels”  they  coincide,  in  the  case  of 
their  anomalous  propagation  it  is  necessary  to  distinguish. 
Horizons/levels  O'O’  and  00  in  Fig.  1c  replace  in  the  wave  picture 
the  maximum  rays/beams  QO"  O'  and  Q0  (Fig.  lb),  obtained  from  the 
geometric  optic/optics. 

These  general/common/total  considerations  will  be  refined  in 
paragraph  4. 
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3.  Basic  formulas. 


Attenuation  factor  V  in  the  unhomogeneous  atmosphere  index  of 
refraction  of  which  depends  only  on  height/altitude,  can  be 
represented  in  the  form  of  the  contour  integral 


l’  (x,  I ( ,  y)  =r  e  *’  4  Y^r  \  e‘"  F  (i ,  u'„  y)  d; .  (3. 0 1 > 


When  the  lowest  layer  of  the  atmosphere  of  inversion  is  present,  if 
one  of  the  corresponding  points  is  located  above  the  layer,  and 
another  -  within  it,  for  integrand  F  it  is  possible  to  take  the 
following  approximation  [see  the  formula  (5.08)  of  (chapter  15]: 

e<  IS  c.-'i— smSui! 


F  u '  -  in  = 


.  n.<(2) 


i  r  iv‘\  —  ;  i  .*•  if.  ■ 


i:<iv)  t- 


S  -  .  nv  7 


-I 
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Here  y'  and  y  -  dimensionless  heights/altitudes  of  source  and 
observation  point  (y*>y,  moreover  u  > y„  and  y  <yt.  where  U-:  — 
dimensionless  spot  height  of  inversion),  _  x  -  the  dimensionless 
horizontal  distance  between  source  and  observation  point,  and  p(y)  - 
the  function,  connected  with  the  given  refractive  index  M(h)  with  the 
formula 

n  lyi  =  \Uh)  -  2rr.-  ( n  —  ]  -  m=  (-y) 3  ,  (3.03) 

moreover  n  is  a  refractive  index  of  air,  a  -  radius  of  the  Earth. 

We  assume  that  function  M(h)  takes  the  same  form  as  in  Fig.  lb 
and  lc.  Therefore  with  this  t  the  equation 

p  (y)  —  i  =  0  (3.04) 

has  two  roots  of  y,  and  ys.  When  />(#,)  </<p(0>  these  roots  are  real 
and  positive,  when  *  </>('/.)  they  compositely  conjugated/combined;  when 
t  =  p  (yi)  they  pour,  and  then  yt  =  yt  =  y:-  Besides  these  two  roots  there 
can  be,  generally  speaking,  and  other  roots  (negative  or  composite), 
but  they  a  value  do  not  have. 
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Values  S(y),  S(y')  and  S„  are  given  by  the  formulas 


DOC  =  62001204 


PAGE 


V 

!  ■»  r~" 


S(y)=  \]'p(y)-tdy,  S(y')  =  \Vp(y)-tdy, 
'■  o 

,  h _ _  *, 

S"  =  — j  I  'p(y)-idy^--L  \  1  'p(y)-,dy, 

<•  - 


(3.05) 


moreover  for  »<p(y«)  radicals  J  p  [y)  —  /  with  real  positive  y  must  be 
taken  in  the  arithmetical  sense.  For  calculation  S„  when  t<p{y;i  it 
is  necessary  to  analytically  continue  radical  ]  p  (y)  —  i  into  the  region 
composit  y.  We  will  consider  that  p(y)  is  analytic  function  [cf. 
below  formula  (3.18)],  that  allows./assumes  this  continuation. 

Value  v  is  determined  by  the  formula 

Jr.  _ 

v  =  |  |  p  (y  i  — ;  du.  (3.06) 

Pi 

For  the  real  values  of  t  value  v  is  also  real,  moreover  sign  v 
is  selected  from  the  following  considerations.  When  y  ^  y-.  function 
p(y)  can  be  replaced  by  the  first  members  of  Taylor  series 

P  <i0  =  P  Or.)  — Oo)  (y  —  p"  (i !'<)  >  °> 

after  which  integral  (3.06)  can  be  calculated,  and  we  obtain  for 
t  ^  P  {yi)  the  approximation  formula 

(3-07) 

12 />'(}/,) 

in  accordancevith  which  we  we  count  O  0  when  t<p\y. »  and 

v  < 0  wnen  !  >  p  Uo-n  when  p  (t/.-i  </  <p  iuj  formula  (3.06)  is  deciphered 
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as  follows: 

v  — - —  I  I  I  —  p  (u)  dtj,  (3.08) 

V\ 


where  1  t  —  p  (y)  >  0,  and  y^yj. 


Function  t  is  determined  by  the  formula 


•/.('•)  = 


l/2ne 


—  7T  v-ri  (v— v  ]  g  v) 


(3.09) 


moreover  when  v>0(/  <p(j/,-))  for  la  v  is  taken  principal  value.  In 
this  case 


'/.(')— ’1  npd‘’v  — oc.  (3.10. 

Key:  (1).  with. 


During  the  calculation  of  attenuation  factor  for  high  values 
Hi  it  is  necessary  to  consider  that  when  y  —  o c  function  p(y) 

must  satisfy  the  relationship/ratio 

lim  [p  (//)  —  «]  =  o. 


(3.11) 
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Therefore,  representing  function  S(y')  in  the  form 


r  * 

Sun  =  f  1  y  —  t dy  -f-  f  [)  p U/)  —  t  —  1  y  —  t]  dy , 

6  o 

we  we  see  that  with  y1-*®  first  term  unlimitedly  grows/rises 
(infinite  part  it  is  equal  to  —//’’*—/  iV.j  and  the  second 
approaches  final  limit,  if  difference  p(y)-y  vanishes  sufficiently 
rapidly  [for  example  in  the  manner  that  for  function  p(y),  determined 
by  formulas  (3.18)  and  (3.19)]. 


Let  us  introduce  value  £0  as  the  limit 


=  iim  :  S  (;/’)  —  2S„  —  -7-  ;/v"  -  ;  |  if'  . 

f  ‘  *  - 


Utilizing  the  approximate  equality 


S  uf)  —  2S„  =  tf'  -  —  J  I  y  ~l„. 


valid  at  the  high  values  of  y',  and  replacing  in  the  denominator  of 
formula  (3.02)  for  value  J "p  uf)  —  i  on  t’7/'.  ve  we  obtain  attenuation 
factor  in  the  form 
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earth/ground  s=a-6  with  the  relationships/ratios 

\_ 

£  =  mO  =  m  -1- ,  m  =  *  (3.17) 

Let  us  note  that  the  point  of  contact  of  tangency  T  corresponds  to 

course  of  ray  in  homogeneous  atmosphere. 

The  infinite  duct/contour  C  in  the  plane  by  complex  variable  t, 
on  which  are  taken  the  integrals  for  V  and  V,,  to  a  considerable 
degree  is  arbitrary  and  must  be  selected  so  that  the  calculation  of 
integral  on  it  it  would  be  possible  to  produce  with  the  smallest 

labor/work,  in  particular,  so  that  the  main  section  of  integration 

would  be  as  far  as  possible  small.  In  this  case  duct/contour  C  must 
cover  entire  of  the  pole  of  integrand  in  the  positive  direction  so 
that  they  would  be  located  above  the  duct/contour  C.  As  the 
duct/contour  of  integration  most  comfortably  it  proved  to  be  to 
select  the  duct/contour,  depicted  in  Fig.  3;  the  salient  point  of 
this  duct/contour  was  arranged/located  either  when  /  =  pi',),  or  several 
more  left  (cf.  the  end/lead  of  paragraph  6). 

Into  integrand  T  (/,  u)  they  enter,  as  can  be  seen  from  formulas 
(3.05)  and  (3.06),  the  integrals  of  form  |  1  r  y/i  —  *  dy  with  different 
t  and  different  integration  limits,  including  of  composite  ones. 
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Fig.  2.  Fig.  3. 


Fig.  2.  Geometric  sense  of  value  $=m$. 

Fig.  3.  Duct/contour  C  in  plane  of  complex  variable 

Page  344 . 

For  facilitating  the  calculations  of  these  integrals  for  the  giver, 
refractive  index  M(h)  was  undertaken  hyperbolic  law  (6.23)  of  chapter 
15,  whence  function  p(  y)  is  obtained  according  to  formula  (3.03)  in 
the  form 

P<M)  =  P<Ui)  ~  7T  »  ,,lh 

f  yi 

moreover  from  relationship/ratio  (3.11)  it  follows  that 

p  (i/,i  =  y:  —  2h.-.  (3. 19'i 

By  force  of  (3.19)  in  formula  (3.18)  are  two  parameters  y,  and  y,. 
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moreover  y,  is  dimensionless  spot  height  of  inversion  expedient  to 
introduce  also  the  special  designation 


then 


Y  =  y,  -  y 


(3.20) 


p'  (&->  =  y-  •  (3.21) 

Let  us  note  that  in  the  case  of  hyperbolic  lav  equation  (3.04)  is 

quadratic  equation  with  two  roots  of  y1  and  ya  flowing  together  when 

l  -  P  <«:)• 


For  the  hyperbolic  lav  integrals  necessary  to  us  are  expressed 
as  the  elliptical  integrals  of  the  first  and  second  kind.  However,  in 
the  cases  examined  by  us  it  proved  to  be  more  convenient  to  calculate 
these  integrals  with  the  aid  of  the  expansions  according  to  the 
degrees  of  parameter  a3,  where 

c.=  _r^.  lo.22, 

In  these  expansions,  which  contain  also  logarithmic 
components/ terms/addends ,  sufficient  to  take  several  first  terms, 
since  the  main  section  of  integration  on  duct/contour  C  corresponds 
to  the  very  low  values  of  parameter  a3.  Further  terms  of  expansions 
are  essential  (at  the  high  values  of  parameter  Y  undertaken  by  us, 
see  the  beginning  of  paragraph  5)  only  in  such  sections  of  the 
ducts/ contours  where  entire/all  integrand  is  already  small. 


In  conclusion  let  us  pause  at  the  analytical  continuation  of 
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functions  F(t,  y',  y)  and  T  a.  y,  to  entire  plane  by  complex  variable 
t.  The  fact  is  that  values  S(y),  S(y'),  S„  and  y  (v)<  entering  these 
functions,  are  initially  determined  only  on  the  real  axis  when 
t<piy)  (v  >  0).  where  for  radicals  \'p  ijy')  —  t  and  'p  {yi  —  t  are  taken 
arithmetical  values.  However,  the  knowledge  of  integrand  when 
/  <  P  U/i>  proves  to  be  sufficient  only  during  the  calculation  according 
to  the  reflecting  formula  (paragraph  4).  For  calculating  the  contour 
integrals  it  is  necessary  to  know  integrand  with  composite  t,  which 
is  reached  with  the  aid  of  the  analytical  continuation. 

In  this  case  it  is  necessary  to  have  in  mind  that  at  point 
t-piyd  the  Precise  functions  F(t,  y',  y)  and  T it.  y\  peculiarities  do 
not  have.  Asymptotic  expression  (3.15)  for  function  ¥(/,«/)  has, 
however,  the  singular  points  (branch  point)  with  t=p(y)  (for 
expressions  }  p  ujj  —  /  and  S(y))  and  with  t=p(0)  (for  expressions 
S(y)  and  $„).  These  singular  points  are  obtained  as  a  result  of  the 
fact  that  we  use  asymptotic  expressions.  In  reality  of  branch  point 
they  are  absent,  since  a  precise  integrand  must  be  meromorphic. 
Therefore  we  is  bypassed  the  "apparent  singular  points"  from  below, 
counting,  for  example,  with  t>p(y),  that  arc  [p(y)-t]  =v  and 
\  pin)  —  i  —  i  I  •  —  p  ij/).  moreover  I  t  —  p  (w)  >0.  This  circuit/bypass  is 
actually  conditional,  since  formula  (3.02)  with  f^p(y)  is  not 
applied  in  view  of  Stokes's  so-called  phenomenon.  To  disregard  this 
phenomenon  is  possible  only  in  the  case  when  section  t>p(y)  gives  a 
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small  contribution  to  the  value  of  contour  integral,  which  in  the 
cases  examined  by  us  occurs.  The  carried  out  by  us  control 
calculations  with  the  aid  of  the  functions  of  the  parabolic  cylinder 
(see  Chapter  15  paragraph  3),  which  give  a  more  precise  asymptotic 
representation  of  subintegral  of  function  V  U,  u),  they  confirmed  not 
only  the  qualitative,  but  also  quantitative  correctness  of  the 
results,  obtained  with  the  aid  of  formula  (3.15). 

Function  H'it.y)  has  also  pole,  corresponding  to  the  roots  of 
equation  (6.01).  In  the  case  when  pole  closely  they  approach  the 
duct/contour  of  integration,  it  is  necessary  to  bypass  them  from 
below. 

4.  Reflecting  formula. 

In  the  illuminated  region  natural  to  calculate  attenuation 
factor  according  to  the  method  of  steady  state,  since  this  method 
gives  transition  to  the  laws  of  geometric  optic/optics,  applied 
sufficiently  far  from  the  horizon/level.  To  integral  (3.14)  the 
method  of  steady  state  can  be  used  as  follows.  Let  us  represent 
integrand  V  on  the  real  axis  in  the  form 

ur_J_  t'ot“ -***“' _  u.O! , 

*  0  4 -  .  , 

"■  $  Ptl'l —  t  ;MV) 


where 
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2  (.•')  =  ;0  —  S  (i'i  —  2S0  —  arc  •/  (vl. 

4» (’.)  =  lo  —  S {y)  -  2 St,  —  arc  x (v)  =  2  (?)  —  2S (u\. 

—  arc  y.  (v)  =  v  lgv  —  v  —  arc  f  (4 —  *v  ) 


A  =  — i— e“a'’’2iS*- 

Z  <v> 


Page  346. 


For  the  entire  integrand  in  integral  (3.14)  it  is  possible  to 
write  with  real  t  the  expression 


„‘:»ur _ L  *  * _ 

e  —  — , _ ~ 

“  i  p{v)— i  x 0  —  •'! 


where 


w(/)  =  ti  —  Q(;i.  j 
q  (;)  =  5/  —  4mm.  i 


Since  in  this  case  v  is  also  real,  first 


and  if  '*>&•  that 


X  (')  =  11  -  e~~'" 


Latter/last  formula  shows  that  when  v>0  a  <  p  <#,-).»  absolute 
value  A  less  than  one  (in  particular,  A  =  — 4^-  when  v  =  Os  and 

rapidly  vanishes  with  increase  v.  Therefore,  if  ve  search  for  the 
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point  of  steady  state  when  KPd/J,  we  can  disregard  the  phase  of 
denominator  1  —  a.  Then  the  points  of  steady  state  t2  and  t2  first  and 
second  term  in  the  right  side  of  formula  (4.05)  are  obtained  from 
equations 

w'  (f,)  =  0,  q/  (/,)  =  0  (4.09) 

or 


-C' (*»).  5  “-<!>' <*i).  (4-10) 

moreover  with  the  given  ones  $  and  y  of  value  t2  and  t2  are 
different. 


Calculations  show  that  the  functions  -SI' it)  and  —o' (?)  have  a 
maximum.  Therefore  we  find  two  values  of  t»  and  two  values  t2  (at 
least,  if  $  is  not  too  great).  Should  be  taken  only  values  of  t2  and 
tj,  which  satisfy  inequa  ity  i  <  P  (£;>  (a*<0),  since  only  for  these 
values  it  is  possible  during  the  determination  of  the  points  of 
steady  state  to  disregard  the  phase  of  denominator  1  —  A. 


After  finding  points  t2  and  t2  we  we  can  calculate  integral 
(3.14),  applying  the  method  of  steady  state  to  each 
component/term/addend  of  formula  (4.05).  In  this  way  we  come  to 
reflecting  formula  for  the  attenuation  factor  V2: 

AHA  «><?(<»)  A  (1A 

V,  (5.  y)=  T~  7  77==^  -  T=  ■■  -?=====  —  .  (4.1<) 

I  p{y)  —  l ,  1  —  21"  (h)  >  p(u)~  it  1  —  ^*1 


where 
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converts/transfers  into  the  usual  reflecting  formula,  which 
escape/ensues  from  the  laws  of  geometric  optic/optics  in 
heterogeneous  atmosphere.  The  latter,  similarly  is  applicable  to  the 
rays/beams,  is  sufficient  to  distant  ones  from  the  maximum  rays/beams 
QO  and  QO"  O'  in  Fig.  lb,  it  is  more  precise,  to  those  rays/beams 
which  have  v<?i)  and  sufficiently  large  positive  numbers.  For 

quite  maximum  rays/beams,  as  it  is  easy  to  consider,  has  v  ~  0,  and 
geometric  optic/optics  to  them  is  not  applied. 

Returning  to  general  reflecting  formula  (4.11),  let  us  introduce 
for  the  maximum  values  -fl'(t)  and  —  4*  <?l  the  designations 

<?)WC.  =■•  [  -4>' 0)i»-Kc- 

In  view  of  formulas  (4.02)  is  3lways  fulfilled  the  inequality 

?,<;«.  (4.  IT) 

Hence  we  see  that  to  find  the  points  of  steady  state  t2  and  t2  for 
both  components/terms/addends  in  formula  (4.05)  is  possible  only  with 
$<$2.  With  $>$2  the  equation  u’(t)=0  does  not  have  real  solution  and 
direct  wave  is  not  expressed  by  first  term  of  formula  (4.11). 
Therefore  value  $«=$2  determines  the  horizon/level  of  direct  waves 
(cf.  paragraph  2).  Analogously  value  $  =  $2  determines  the 
horizon/level  of  the  waves,  reflected  from  the  earth/ground. 

Physical  sense  $ 2  consists  in  the  fact  that  into  the  region  $>$2 
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the  electromagnetic  waves  are  drawn  through  only  via  diffraction; 
therefore  there  is  boundary  of  the  region  of  shadow. 

Page  348. 

Physical  sense  lies  in  the  fact  that  with  $<$!  is  applicable 
reflecting  formula  (4.11),  in  consequence  of  which  $=$i  there  is  a 
boundary  of  the  illuminated  region.  Region  $*<$<$,  is  the 
intermediate  region,  included  between  both  horizons/levels. 

Since  the  maximum  values  of  functions  -fi'(t)  and  —<!>'(/)  attain 
near  point  i  ~  p  (*/,-),  that  value 

S»«l-Q* (Oil-,,,,..  L- 1 -O' <01.-,.*.  (418> 

they  will  be  very  close  to  the  values,  determined  by  formulas  (4.16); 
we  will  show  this  based  on  examples  in  paragraph  5.  Therefore  the 
position  of  the  horizons/levels  it  is  possible  approximately  to 
determine  from  formulas  of  type  (4.18)which  is  much  simpler  than  to 
construct  plotted  functions  -fi'(t)  and  —<t>' (/;.  necessary  with  the  use 
of  formulas  (4.16)  For  hyperbolic  law  (3.18)  of  formula  (4.18)  they 
lead  to  the  expressions 

s,  =  G<,  —  G  («);  I.  —  G„  —  G  (//i,  ^4. 1 9i 


where 
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G„  = 


■  U;  - XX lg.LLni^^ -  XX  rc,  -  X lg (K», ; ;  (4.20) 

I  >  —  l  u, 


G  Uj )  =  —  i  yi  —  y  -  I  y}  — 


..  -LLL-0:.-  'X 


!  ig  -4-^-”==  -  ig  X^— XX 

i  6  VY  -  \r u  -  u  K  \  Y~Vy; 


(4.21) 


and 


Ci  =  C  —  7  lg  2  —  4  =  1 ,429 


(4.22) 


(C  -  Euler's  constant). 


Second  formula  (4.19)  upon  transfer  to  the  usual  (dimensional) 
coordinates  gives  formula  for  the  horizon  distance  of  the  waves 
[formula  (6.32)  reflected  of  Chapter  15].  However,  first  formula 
(4.19)  defines,  as  we  already  spoke,  the  horizon  distance  of  direct 
waves . 


Let  us  note  in  conclusion  that  reflecting  formula  (4.11)  is 
applicable  for  calculating  the  attenuation  factor  \\  almost  to  the 
very  horizon/level  of  direct  waves  Gj. 

Page  349. 

5.  Numerical  results  in  the  dimensionless  coordinates. 


iwHiwiwiifjrjl 
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During  the  calculation  of  attenuation  factor  for  the 
hyperbolic  law  of  inversion  we  took  the  following  numerical  values  of 
the  parameters,  entering  function  p(y)  [formula  (3.18)-(3.20) ] : 


Table  1. 


Functions  p(y)  at  the  selected  values  of  the  parameters  are 
depicted  in  Fig.  4.  The  selection  made  by  us  allows  with  the  defined 
M-prof ile/airfoil  (see  paragraph  7)  to  calculate  propagation  of  four 
wavelengths  which  relate  as  by  1:3:9:27.  In  this  case  the  first  row 
Table  1  corresponds  to  shortest,  and  fourth  -  to  longest  wave. 


In  all  cases  we  took  y=-j-.  i.e.  was  assumed  height/altitude  of 
one  of  the  corresponding  points  the  equal  to  one  fifth 
height/altitude  of  the  layer  of  inversion.  Another  point  we  took  at 
the  high  altitude  above  the  layer  of  inversion  -  such  large  that  it 
is  possible  to  use  the  attenuation  factor  V^S,  y) ,  connected  with  V 
according  to  formula  (3.13). 


The  four  curves  calculated  by  us  for  the  attenuation  factor  V, 
depending  on  variable/alternating  $  are  given  in  Fig.  5  1 . 
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FOOTNOTE  1 .  On  ail  drawings,  carried  out  on  the  logarithmic  scale, 
along  the  logarithmic  axis  are  plotted  decimal  (but  not  natural) 
logarithms.  ENDFOOTNOTE. 

Indices  1,  2,  3  and  4  in  the  curves  they  show,  to  what  row  table  to  I 
and  what  curve  in  Fig.  4  corresponds  this  curve  for  the  attenuation 
factor.  For  each  curve  the  point  Gx  notes  the  position  of  the 
horizon/level  of  the  direct  waves,  and  point  G2  -  position  of  the 
horizon/level  of  the  waves,  reflected  from  the  earth/grouna.  Points 
G0  near  the  origins  of  coordinates,  equipped  with  the  same  indices  1, 
2,  3  and  4,  determine  the  horizon/level  (boundary  of  straight/direct 
visibility)  with  homogeneous  atmosphere  the  corresponding  values 
they  are  obtained  from  simple  formula  »«  =  1  ~y. 


DOC  =  82001205 


PAGE  Z 


Fig.  4.  Plotted  functions  r  »«••>/-  for  the  values  of  the  parameters 

Table  1. 

Page  350. 

As  is  evident,  into  all  four  cases  examined  occurs  the 
hyperdistant  propagation,  most  strongly  expressed,  as  this  it  was  to 
be  expected,  in  curve  1.  Upon  transfer  to  curved  2,  3  and  4 
phenomenon  of  hyperdistant  propagation  monotonically  weakens; 
however,  on  curved  4  with  $~~5  attenuation  factor {V (proves  to  be 
order  0.1,  while  with  the  same  $  and  y,  but  in  homogeneous 
atmosphere,  |V*  \  takes  values  for  four  orders  below  (Jv,j — 0.000013). 

Table  2  gives  the  values  of  functionjVx|  at  the  horizons/levels 
Gx  and  G, . 

From  it  it  is  evident  that  the  values  of  attenuation  factor  at 
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both  horizons/levels  G,  and  G2  change  within  sufficiently  vide  limits 
-  3-3.5  times.  For  the  comparison  table  2  gives  valuesjV2|at  the 
horizon/level  Gc  during  the  normal  propagation  and  at  the  same  values 
of  y.  The  comparison  of  columns  shows  that  during  the  normal 
propagation  of  the  value  of  attenuation  factor  at  the  norizon/level 
they  have  because  of  the  dependence  on  y  much  the  same  spread,  as 
during  anomalous  propagation  at  the  horizons/levels  G2  and  Ga. 
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25  30  35  hC  Hi  50  r 


Fig.  5.  Dependence  of  attenuation  factor  V,  on  5.  The  numbers  of 
curves  1,  2,  3,  4  correspond  to  the  numbers  of  revs  Table  1  and  to 
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begins  to  monotonically  decrease  already  in  the  illuminated  region, 
it  is  more  left  both  horizons/levels.  This  leads,  in  particular,  to 
the  fact  that  the  factor  of  attenuation  at  the  horizons/levels  G*  and 
Gj  proves  to  be  on  Table  2  2-4  times  less  than  during  the  normal 
propagation  at  the  horizon/level  Ge.  This  behavior  of  attenuation 
factor  is  explained,  apparently,  by  the  fact  that  not  only  beyond  the 
horizons/levels  G*  and  Gj  but  also  are  more  left  them  they  have  a 
value  the  diffractive  (are  more  precise,  wave)  phenomena,  considered 
by  reflecting  formula  (4.11)  and  which  are  not  placed  in  the  laws  of 
geometric  optic/optics. 

For  explaining  the  applicability  of  simple  formulas 
(4.18)-(4.22)  for  calculating  the  horizon  distance  G,  and  G*  let  us 
compare  in  the  cases  examined  by  us  the  results  given  by  them  with 
the  results  according  to  formulas  (4.16). 


Table  3  shows  that  both  formulas  give  very  close  numbers. 
Therefore  for  the  practical  calculations  of  horizon  distance  it  is 
possible  to  use  the  simple  formulas  of  Chapter  15. 
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Table  3. 


St 

k! 

%» 

*> 

j 

49.:! 

40  i  j 

52,26 

52.2;. 

O 

,  28.56 

2o,5 o 

30.74 

3C.74 

3 

16.06 

15.99 

17.52 

17.50 

4 

6.67 

i 

6,45 

9.52 

9.50 

6.  Attenuation  factor  in  a  deep  shadow.  Series/row  of  deductions. 

Attenuation  factor  in  a  deep  shadow  is  conveniently  investigated 
with  the  aid  of  the  series/rov  of  deductions  which  is  obtained  from 
integral  (3.14)  in  the  usual  way  (cf.  chapter  14,  paragraph  O.  In 
order  to  obtain  the  series/row  of  deductions,  it  is  first  of  all 
necessary  to  determine  a  precise  position  of  poles  of  function 
^ J).  i.e.  the  roots  of  the  equation 

1  —  A  =  0.  (6.01.1 

These  roots  are  located  near  duct/contour  C  (Fig.  3)  or  within  it. 

Page  352. 

If  we  designate 

the  value  At  for  the  roots  found  by  us  they  are  given  in  Table  4 
whose  first  column  shows  the  number  of  row  in  Table  1,  and  the  second 
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-  number  of  root  for  the  present  instance. 

The  position  of  the  real  parts  three  first  roots  relatively 
curved  p(y)  for  the  first  case  is  represented  in  Fig.  6.  We  see  that 
only  the  first  root  corresponds  to  the  "seized"  wave  in  the  usual 
understanding,  other  two  roots  give  the  waves  which  from  the  point  of 
view  of  geometric  optic/optics  easily  would  exceed  the  limits  of  the 
layer  of  inversion.  However,  these  "being  drawn  through"  waves  have 
small  fading  and  actively  participate  in  the  process  of  hyperdistant 
propagation.  Let  us  recall  that  during  the  normal  propagation 
t1=1.17+i2.D2,  so  that  in  this  case  the  third  wave  attenuates  10 
times  more  slowly  than  least  damped  wave  under  normal  conditions  for 
propagation.  For  the  remaining  cases  all  roots  correspond  to  the 
"being  drawn  through"  waves. 

We  convert  equation  (6.01>  to  the  simple  approximate  form,  which 
allows/assumes  comparison  with  other  theories  of  hyperdistant 
propagation.  Let  us  begin  from  the  "seized"  waves  which  nave  almost 
real  t,  which  lie  between  ?($/,•)  and  p(0)  (as  the  first  root  in  Table 
4),  and  therefore  negative  values  v.  For  v<0  we  will  assume 


v  =  (— v)  eia.  Jg  v  =  lg'(— v)  —  in. 


(6.03) 
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Table  4. 


&  m 

Sz 

" 

•  i 

0.10653  —  i  0.00019 

i 

o 

-0.1733-  i  0.3293 

1  — 0.06364  t  i  0,05523 

3 

T  i 

—0.1036  -  i  0.2236 

1  3 

,  —0.1633  -  /  0.2107 

2  i 

—0.1863—  i  0.6934 

.  A 

—0.2495  —  i  0.3913 

*» 

1  * 

—0.0652  —  i  0.4661 

*  i  i 

s 

i  — 0.06336  —  i  £*.06516 

! 

o  ; 
*  j 

-0.1271 -‘i  i.1316 

Fig.  6.  Roots  w..  corresponding  to  seized  and  unseized  waves  for  data 

P  (w  —  p  yy.) 


Page  353. 


Then  together  with  (3.10)  we  will  have 


XM-1  with  v - 05 


(6.04) 


and  from  formula  (3.05)  we  obtain 
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where 

Vi 

•Si  =  J  1  P(ii)  —  i  dy.  (6.06) 

bv4-  Vi  indicates  the  smallest  positive  root  of  equation  (3.04). 

.'R'.i.  y  into  account  these  formulas  equation  (6.01)  takes  the  form 

ie-‘Si  =  y  (v).  (6.07) 

If  v  is  great  and  is  negative  ‘(strongly  seized  waves),  then  in 
view  of  relationship/ratio  (6.04)  we  obtain  the  simpler  equation 

S,  =  ( m - m=  I.  2,  .  ..  (6.0S> 

corresponding  to  the  known  characteristic  equation  of  the  seized 
waves. 

Let  us  visualize  now  that  v  positively  or  is  composite  with  the 
positive  real  part,  i.e.,  Re/-',  p  (.//,>  or  Re  At<0.  In  this  case  value  Si 
to  determine  with  the  aid  of  formula  (6.06)  is  impossible  at  least 
because  it  is  unknown,  which  of  the  composite  roots  of  Yi  and  Yi  one 
should  take.  However,  turning  formulas  (6.05),  we  can  always 
determine  S  with  the  aid  of  the  relationship/ratio 

■Si  =  S«  —  4-  v.  (6.09) 

and  then  from  equation  (6.01)  we  again  obtain  equation  (6.07).  With 
that  being  proper  [i.e.  corresponding  to  formula  (6.03)]  selection 
v  oo  we  always  have  Z  (with  exception  of  case  arc  v  - — jrj- 
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Since,  furthermore,  */  »0)  -»|  1,  that  for  the  "being  drawn  through" 
waves  it  is  possible  to  count  in  the  first,  roughest, 
approximation/approach  i  tv)  =  1  •  and  we  obtain  equation  (6.08). 


Let  us  note  that  simplified  equation  (6.08)  usefully  also  for 


the  normal  propagation  when  in  relationship  (6.06)  we  obtain 
p  (y)  =  ii  olhJ  pi  -  i.  Thus,  from  (6.08)  we  obtain 


i 

*  m 


l  \  13/ 

e 


n 

3 


(6.10) 


which  approximately  corresponds  to  the  roots  of  characteristic 


equation  for  homogeneous  atmosphere. 


Page  3b*. 

For  checking  equation  (6.08)  we  calculated  (table  5)  for  the 
roots  found  by  us  value  Si  according  to  formula  (6.09),  and  as  a 
result  were  obtained  these  numbers: 


table  5. 
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Jw 

m 

j  | 

s, 

i 

j  v 

1 

] 

2.350 

2.326  —  i  0.00! 

i  —0.768  -  i 

0.0014 

0 

5.498 

5.537  -  i  0.047 

,  0,459  —  i 

0,396 

3 

8.639 

8.646  -  i  0.009 

■  1,178—  i 

1.519 

H.781 

1 1,764  — i  0,005 

i  1,800  — i 

2,821 

O 

1 

2.356  . 

2.444  -  i  0,062 

i  0.317  —  t 

0,376 

9 

5,496 

5.50!  —  i  0.0  il 

i  0,867  _  i 

1.646 

3 

i 

,  2.356  J 

2.315  -r  i  0.0!  1 

0,360  -  i 

0.776 

9 

5.498  1 

5.516  —  i  0,014 

1  0.656  —  i 

2,402 

4 

i 

2.356  i 

2.436  -  i  0,079 

1  0,207  -  i 

1.120 

2 

5,496 

5.499  -  i  0,007 

!  0,319  -t 

2,718 

Thus,  calculating  SA  for  the  obtained  root,  we  with  the  aid  of 
approximate  relationship/ratio  (6.08)  can  ascribe  to  it  number  m. 


Fig.  7  depicts  the  attenuation  factor  in  a  deep  shadow, 
calculated  according  to  the  series/row  of  deductions  for  the  first 
case.  Fig.  7  shows  that  the  first  term  of  the  series/row  of 
deductions,  which  corresponds  to  pole  tj,  determines  attenuation 
factor  only  with  $>150,  i.e.,  for  the  wave  X=1  cm  with  s>1000  km. 
Since  the  first  term  has  negligible  fading,  then  at  such  large 
distances  the  absolute  value  of  attenuation  factor  will  be  almost 
constant  -  asymptote  in  Fig.  7  .is  almost  horizontal.  Let  us  note  that 
in  a  deep  shadow  in  Fig.  7  attenuation  factor  approaches  an 
asymptote,  completing  dying  oscillations.  These  oscillations  are 
caused  by  the  interference  first  and  of  the  second  "simple  wave". 
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6 
■o.c 

to 
?,5 
-2,0 

Pig.  7.  Dependence  on  $  that  computed  from  the  series/row  of  the 
deductions  of  attenuation  factor  Vj  in  a  deep  shadow. 

Key:  (1).  deduction. 

Page  355. 

Thus,  the  first  simple  wave,  which  has  the  smallest  fading,  is 
very  weakly  excited  by  the  wave,  which  falls  on  top  to  tropospheric 
waveguide,  thanks  to  which  tins  simple  wave  can  be  decisive  only  at 
the  very  large  distances.  Near  the  horizons/levels  Gj  and  G2  primary 
meaning  has  the  second  and  partly  third  term  of  the  series/row  of 
deductions.  This  phenomenon  must  have  general  character,  since  if 
simple  wave  is  "seized"  (see  above)  and  barely  it  is  drawn  through 
from  the  layer  of  inversion  (how  and  it  is  explained  its  negligible 
fading),  then  for  reasons  of  reciprocity  it  is  barely  excited  by  the 
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emitters,  arranged/located  higher  than  the  layer  of  inversion.  The 
waves,  which  possess  high  fading,  to  the  larger  degree  penetrate  the 
space  above  the  layer  of  inversion;  therefore  they  are  excited  more 
strongly  and  play  the  dominant  role  near  the  horizons/levels. 

Because  of  the  noted  fact  the  horizons/levels  Gj  and  Gs 
actually/really  define  (although  in  the  sufficiently  approximate 
sense)  the  superstandard  range  of  radio  waves  even  during  the 
strongly  expressed  superref raction  (as  can  be  seen  from  Fig.  7). 

Usually  in  the  examination  of  hyperdistant  propagation  they  are 
based  on  the  series/row  of  deductions.  In  this  case  they  assume  that 
only  seized  waves  (ReA/m>uj  can  have  a  small  fading.  In  fact  and  the 
waves  "being  drawn  through"  (ReA/m<0)  in  a  number  of  cases  also 
attenuate  weakly.  Therefore  waves  several  times  longer  than  the 
"critical"  wavelength  X0f  determined  according  to  Bremmer  [25],  it  is 
still  capable  of  the  hyperdistant  propagation  in  the  tropospheric 
waveguide. 

Let  us  note  in  conclusion  that,  as  showed  calculations,  several 
first  terms  of  the  series/row  of  deductions  already  make  it  possible 
to  calculate  attenuation  factor  by  Poti  up  to  the  joining  with  the 
reflecting  formula  and  they  free  thus  from  the  calculations  on  the 
quadratures  (cf.  paragraph  3). 
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7.  Numerical  results  for  the  specific  case. 

For  facilitating  the  physical  analysis  of  the  numerical  results, 
obtained  by  us  in  paragraph  5,  we  will  examine  the  here  corresponding 
to  them  specific  case. 

As  an  example  we  will  take  an  M-prof ile/airfoil ,  depicted  in 
Fig.  8,  and  will  construct  attenuation  factor  Vx  for  the  following 
waves:  1)  3.33  cm,  2)  10  cm,  3)  30  cm,  4)  90  cm.  Graph/curve  V,  is 
given  in  Fig.  9.  The  numbers  of  curves  Fig.  9  indicate  the 
wavelengths  enumerated  here.  Along  the  axis  of  abscissas  we 
plot/deposit  on  the  lower  scale  distance  of  s  in  the  kilometers,  and 
on  the  upper  -  angles  v  in  the  degrees  (cf.  Figure  2).  Along  the 
axis  of  ordinates  are  deposited/postponed  the  common  logarithms,  and 
on  the  right  scale  are  noted  other  values  |,VV- 

Let  us  note  that  in  our  calculations  is  not  taken  into 
consideration  the  dispersion. 
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We  assume  that  an  M-curve  takes  one  and  the  same  form  for  all  four 
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waves  for  which  Fig.  9  gives  attenuation  factor  vx. 


M-curve  in  Fig.  8  is  constructed  according  to  the  hyperbolic  law 

,  0-01) 


in  which 


;\1  {h;)=zizJ!iL.  (7.02) 


Hyperbolic  law  encompasses  two  parameters:  h=  and  1,  having  the 
dimensionality  of  height/altitude  and  connected  with  dimensionless 
constants  and  /,.  in  formula  (3.18)  by  the  relationships/ratios 

..  -  kh-  * 


177 


(f)T. 


r?.03i 


moreover  h%  is  spot  height  of  inversion  or,  which  is  the  same,  the 
height/altitude  of  radio  duct.  The  height/altitude 


H  ~  h 


,  fy  kH  \ 


(7.04. 


as  it  is  easy  to  show,  is  determined  radius  of  curvature  of  M-curve 
at  the  inversion  point. 


In  Fig.  9  distances  are  expressed  in  the  kilometers.  On  the  axis 
of  abscissas  is  noted  also  the  horizon/level  G„  during  the 
propagation  in  homogeneous  atmosphere. 
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Fig.  8.  Dependence  on  height/altitude  h  of  the  given  refractive  index 
(M-prof ile/airfoil)  for  the  values  of  the  parameters: 

h  -  4 c.5  v  /  H  k  5. 

A)  (‘fc)  =  !5i.=-  _  •.)  (Oi  -  A!  U  )  ^  fc.3H 


Fig.  9.  Dependence  of  attenuation  factor  Vx  on 
for  wavelengths:  curve  1  -  3.333  cm;  curve  2  - 
cm;  curve  4-90  cm. 


distance  of  s  (<*-  km) 
10  cm;  curve  3-30 
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This  horizon/level  is  determined  by  spot  height  of  observation  h  and 
does  not  depend  on  wavelength:  let  us  recall  that  is  everywhere 
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undertaken  For  each  curved  the  point  Gx  determines  the 

position  of  the  horizon/level  of  direct  waves,  and  point  Ga  - 
position  of  the  horizon/level  of  the  waves,  reflected  from  the 
earth/ground  (paragraphs  2  and  4).  The  horizons/levels  Gx  and  Ga  are 
changed  with  a  change  in  the  wavelength  and  therefore  for  each  curve 
their. 

In  all  cases  it  is  possible  to  state/establish  the  phenomenon  of 
hyperdistant  radiowave  propagation,  which  weakens  with  an  increase  in 
the  wavelength.  Taking  into  account  a  strong  change  in  the  wavelength 
upon  transfer  from  one  curve  to  another  (wavelength  they  relate  as  by 
1:3:9:27),  one  should  recognize  that  near  the  horizons/levels  the 
dependence  of  attenuation  factor  on  the  wavelength  is  comparatively 
weak. 


Into  formula  for  the  horizon  distance  (cf.  chapter  15,  paragraph 
6)  the  wavelength  enters  only  under  log  sign.  Therefore  horizon 
distances  form  the  arithmetical  progression,  if  wavelengths  as  in 
Fig.  9,  form  geometric  progression.  In  this  case,  however,  the  value? 
of  attenuation  factor  on  both  horizons/levels  Gx  G3  depends  from  the 
wavelength  to  the  same  degree,  as  during  the  normal  radiowave 
propagation  (cf.  Table  2). 

In  view  of  these  facts  to  identify  the  distance  of  propagation 
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of  radio  waves  with  the  horizon  distance  of  the  straight/direct  or 
reflected  waves  is  necessary  with  certain  precaution.  It  is  possible 
to  define  superstandard  range  otherwise,  for  example  as  such  distance 
at  which  the  attenuation  factor  has  absolute  value  of  0,  1,  moreover 
on  greater  the  distances  of  the  value  of  attenuation  factor  it  is 
still  less.  During  the  iatter/last  definition  the  "superstandard 
range"  is  included  between  the  horizon  distances  Ga  and  Ga  for  curve 
1  in  Fig.  9,  and  for  another  curved  this  distance  is  less  than  the 
distance  Gx;  in  the  rough  approximation  and  these  four  distances  form 
the  arithmetical  progression,  as  can  be  seen  from  figure.  Let  us  note 
that  for  evaluating  the  superstandard  range  according  to  value  of  0.1 
it  usually  proves  to  be  sufficient  to  produce  calculations  with  the 
aid  of  the  reflecting  formula  of  paragraph  4,  only  sometimes  applying 
the  extrapolation  of  curves  obtained  in  this  way. 

The  direct  purpose  of  the  reasonings  of  this  chapter  consisted, 
as  we  already  spoke  in  paragraph  1,  in  checking  of  formulas  for  the 
superstandard  range  of  radio  waves,  derived  in  the  preceding/previous 
chapter.  Above  we  showed  that,  introducing  the  horizons/levels  of  the 
straight/direct  and  reflected  waves,  it  is  possible  to  obtain  the 
simple  and  demonstrative  picture  of  hyperdistant  radiowave 
propagation  when  the  layer  of  inverse  is  present,.  However,  the 
distance  of  propagation  only  in  the  sufficiently  rough  sense  can  be 
identified  with  distance  of  one  of  the  horizons/levels. 
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The  fact  is  that  the  decrease  of  attenuation  factor  (after  the 
termination  of  oscillations  in  the  illuminated  region)  begins  earlier 
than  we  is  reached  first  horizon.  Therefore,  as  shown  in  paragraph  5, 
attenuation  factor  Vj  accepts  at  the  horizon/level  Cx  and  Ga  of 
value,  2-4  times  it  is  less  than  at  the  usual  horizon/level  G0  during 
the  propagation  in  homogeneous  atmosphere.  Furthermore,  near  the 
horizons/levels  G.  and  G-  attenuation  factor  decreases  it  goes 
without  saying  it  is  much  slower  than  during  the  normal  propagation. 


All  these  reasons  lead  to  the  fact  that  the  horizons/levels  Gj 
and  Gj  during  anomalous  propagation  characterize  the  superstandard 
range  of  radio  waves  more  roughly  than  the  horizon/level  G„  during 
the  normal  propagation.  However,  the  possibility  of  applying  the 
horizons/levels  Ga  ana  G2  for  the  approximate  estimates  of 
superstandard  range  does  not  produce  the  doubts,  as  is  evident  at 
least  from  the  comparison  of  attenuation  factor  near  the 
horizons/levels  and  of  a  deep  shadow  in  Fig.  7. 


f>rc-f «/e 

One  should  emphasize  that  selected  by  us  has 

sufficiently  weak  inversion:  difference  'HO 1  —  AJ  (ft.-)  does  not  exceed 


1 


DOC  =  82001205 


PAGE 


several  tenths.  This  inversion  in  certain  cases  can  in  practice 
remain  unestablished.  However,  our  calculations  show  that  even  this 
M-prof ile/airfoil  strongly  changes  the  character  of  radiowave 
propagation,  leading  to  the  hyperdistant  propagation. 
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Chapter  17. 

Radiowave  propagation  on  the  surface  tropospheric  waveguide  l. 

FOOTNOTE  1 .  Foch,  Weinstein  Belkina ,  1958.  ENDFOOTNOTE. 

A 

Chapter  is  dedicated  to  the  theory  of  radiowave  propagation 
between  the  corresponding  points,  which  are  located  in  the  lowest 
layer  of  the  atmosphere  of  inversion  (in  the  tropospheric  waveguide). 
To  this  case  are  applied  the  general  formulas,  derived  in  Chapter  15, 
moreover  they  are  utilized  both  the  functions  of  parabolic  cylinder 
and  Airy's  function.  With  the  aid  of  this  procedure  is  conducted  fcr 
a  series/row  of  specific  cases  the  calculation  of  the  radiowave 
propagation  of  different  length.  From  the  obtained  results  it  follows 
that  nyperdistant  (waveguide)  propagation  weakens  with  ar.  increase  in 
the  wavelength  only  very  slowly.  Thus,  for  instance,  wavelength  can 
exceed  the  so-called  "critical"  wavelength  by  an  order,  and 
hyperdistant  propagation  will  occur.  At  the  end  of  the  chapter  is 
more  precisely  formulated  the  criterion  of  hyperdistant  propagation. 


Introduction. 
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Present  Chapter  is  dedicated  to  the  theory  of  radiowave 
propagation  in  the  surface  tropospheric  waveguide  (layer  of 
inversion)  on  the  assumption  that  both  corresponding  points  are 
located  within  the  waveguide.  This  propagation  of  radiowaves  can  be 
named/called  intra-layer  or  intra-waveguide  ir.  contrast  to  the  case 
(examined  in  Chapter  16),  when  one  of  the  corresponding  points  is 
located  highly  above  the  layer  of  inversion. 


To  the  investigation  of  radiowave  propagation  in  the 
tropospheric  waveguide  is  dedicated  the  serios/rov  of  the  theoretical 
works  (for  example,  see  [23-25]);  however,  on  this  question  there  are 
still  many  vaguenesses.  In  this  chapter  we  investigate  intra-layer 
propagation,  developing  the  methods,  presented  in  Chapters  15  and  16. 
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1.  Basic  formulas, 


In  this  paragraph  we  will  compare  the  basic  formulas,  obtained 
in  Chapters  14,  15  and  16. 


The  attenuation  factor  of  the  V  with  the  arbitrary  dependence 
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refractive  index  of  the  atmosphere  on  the  height/altitude  is 
determined  by  means  of  the  contour  integral 

- -  7T  . 

\’u.  y,  tj\  q)  =  1' I e;x:FU,  y.  y',  q)di .  (1.0!) 

t  •>  £ 

where  duct/contour  C  covers  in  the  positive  direction  entire  of  the 
pole  of  integrand  (Chapter  14). 

We  are  limited  to  value  of  q=®,  which  corresponds  to  arbitrary 
polarization  on  the  decimeter  and  shorter  waves  and  to  horizontal 
polarization  on  the  waves  of  of  longer.  In  this  case  function  F(t,  y, 
y 1 ,  ®)=F(t,  y,  y')  takes  the  form  (for  y'>y): 

FU,  y.  y')= - t)  [kU’  n - fi  (y-  0_  .(1.021 

Here  y  and  y*  essence  the  dimensionless  heights/altitudes  of  the 
corresponding  points: 


m 


w: 

n. 


( 1.03) 


value  x  in  formula  (1.01)  is  the  nondimens ional  distance  between 
these  points,  counted  along  the  earth's  surface: 


x  = 


ks 


(1.04) 


but  parameter  m  in  the  latter/last  formula  has  a  value 

/  to.  \  1/3 

m={— )  •  (1.05) 

where  a  is  a  radius  of  the  Earth. 


DOC  =  82001205 


PAGE  ]<?£ 


Functions  fx  and  f2  -  the  independent  solutions  of  the 
differential  equation 

-jp--|PW  —  *)/  =  0.  (1.06) 

moreover  their  Wronskian  is  assumed  to  be  equal  to  -2i.  Function  p(y) 
is  connected  with  the  given  refractive  index  M(h)  with  the 
relationship/ratio 

P  (V)  =  -w-  M  W  =  2m“*  (n-l-v)  (l'07) 

(n  -  the  refractive  index  of  air)  and  is  its  dimensionless  analog. 
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In  Chapter  15  equation  81,06)  integrated  asymptotically  on  the 
assumption  that  p(y)  has  one  minimum,  i.e.,  that  the  given  refractive 
index  of  the  atmosphere  has  one  inversion  point.  This  character 
carries,  in  particular,  the  refractive  index,  depending  on 
height/altitude  according  to  the  hyperbolic  law 

•V!  (/!)  =  A! (h.)  -  -Y~^r  10-  ( 

(ftf—  the  height/altitude  of  inversion;  1  -  parameter),  to  which 
corresponds  the  function 


p  (i/)  =  p  (yP 


( 1 .09 1 
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This  function  is  accepted  in  all  calculations  both  in  the 
present  Chapter  and  in  Chapter  16. 


Attenuation  factor  V,  thus,  in  addition  to  dimensionless 
coordinates  x,  y,  y',  depend  on  parameters  y,  and  ."/  (or  Y= 

«)•  characterizing  function  p(y). 


As  shown  in  Chapter  15,  the  high-altitude  factors  f,  and  f ,  are 
expressed  as  the  functions  of  parabolic  cylinder  D>-  as  follows: 


m  (w.  ? i  =  r,  (vi  1  2  -jrgi  O 
/;(//.  /)  =  C..  (V  )  | 7  2  g:  Qi.  j 


( l.ll'l 


where 


£i  t 


(  — .  i1  \  ! 
D  =  d  ,  v.e  <:/. 


g;(;.  =  D 


(.•  ; 


{Mb 


Expression  f or  Dr.  (z)  in  the  form  of  series/row  is  given  below  [formula 
(2.05)]. 


Value  v  is  determined  by  the  relationship/ratio 

(1.12i 


V* 


v  =  — —  I  1  t  —  p(y)dy . 


y\ 


(y,  and  y,  are  the  roots  of  equation  p(y)-t=0).  As  a  result  of 
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p"  (yi)  > 0  it  is  possible  to  accept,  that  v>0  when  t<p(y,);  with  t> 
P  did  is  obvious  v<^®- 
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Variables  $  and  y  are  connected  with  "the  substitution 

%  _ 

S(y)-S0=  fj/4-P-rv^.  (1.13) 

0 

where 

t 

S(y)  =  I  I  p(y)  —  tdy ;  (1.14^ 

t  7  f  #■* 

S"  =  ill  P(y)  —  tdy~4-  I  I  'p(y)-tdy,  "  (1.15) 

t-  ~  ii 

and  finally 

i4v 

Cj(v)=e  <  6e  * 

•n  .  n  :  _i 

C;  (v)  =  e  4  fc'e  1 

For  the  large  positive  ones  $,  which  corresponds  to 
heights/altitudes  y,  arranged/located  higher  than  the  inversion  point 
y{  and' sufficiently  far  from  it,  high-altitude  factors  allow/assume 
the  asymptotic  representation 


•■rv  lev— s>t 


_4  vi£  %  -  Vi 


(1.16) 
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My.  i)  =  e- 

>  p  (y)  —  t 

A  <y.  /)  =  r~ —  -■  e-'  . 

t  p(y)  —  t 


(1.17) 


Lower  than  inversion  point  and  far  from  it  ($  is  great  and  is 
negative)  the  asymptotic  representation  f*(  y,  t)  and  f4(y,  t) 
the  form 


A  (y-  0  =  y.j  (v)  <-4,-%.  e-'  - 

,  piyi  —  i 

U  (y-  A  =  Xs  (V)  7- — —  e-'  <*-*.> 
»' p(y)  -  < 


e  '  « 

.  e-n'  _ _ p-is 

1  P  (£.’)  —  / 

.  a 
e‘  « 

-j_e-jtv_£ - eiS> 

I'p  (if)  —  i 


(1.18) 


where  it  is  placed 


X,  (v)  =  —  e--f + mv  -  v  .c  v, 


'(+-'•) 


fc(v)=-  . ■  e 


(iH 


31V 

— - *  (V  —  V  Ip  V) 


(1.19) 
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takes 


Substituting  expressions  (1.18)  into  formula  (1.02),  we  obtain 
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function  F(t,  y,  y’)  in  the  case  of  intra-layer  propagatioi \iu  y' '<$yt): 


F(t,  y.  u')  =  4-  — 


i 


V  p  (y)  — 1\'  p  (y')  —  / 

|e!“  U'l  _Ae-iS  ((.-)]  Ie-/S(*)_e/S(,,)1 

x  r^A 


(1.20) 


where 


A  =  — —  e— nv-:'s* . 
/!  (v) 


(1-21) 


Attenuation  factor  V  was  calculated  from  the  series/row  of  the 
deductions 


l’(x.  y.  y')  =  '2\  rue  4  ^ 


(1.22) 


Through  R„  is  designated  the  deduction  of  function  F(t,  y,  y’)  in  m 
pole  >  im.  If  we  use  approximation  formula  (1.20),  then  zm  is  the  m 
root  of  the  equation 

1  —  A  =  u,  ( 1 .2o) 

and  deduction  Rm  takes  the  form 


Rm=  2i  — 


si!i  Sjj'i  sin  5  (;••') 


(1.24) 


v'pim  —  tr-.l'p  (y'l  -  im  ^  d:  '■ 1  A  . 


The  majority  of  numerical  results  in  this  work  is  obtained  with 
the  aid  of  formulas  (1.22)— (1.24) ,  moreover  we  were  limited  by  three 
or  four  members  of  the  series/row  of  deductions. 


2.  Calculation  of  attenuation  factor  with  the  aid  of  the  functions  of 
parabolic  cylinder. 
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Asymptotic  representations  (1.18)  of  the  high-altitude  factors 
f*  and  f,  on  basis  of  which  are  obtained  formulas  (1.22)-(1.24) ,  are 
valid,  if  $  is  great  and  is  negative  or  it  is  located  in  certain 
sector  around  the  negative  real  semi-axis. 
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In  the  latter  case  absolute  value  $  must  be  large.  However,  as  can  be 
seen  from~”£abie  1,  in  a  number  of  cases  interesting  us  appears  the 
serious  doubt  of  the  applicability  of  asymptotic  formulas,  since 
value  $  proves  to  be  on  the  modulus /module  of  the  order  of  one  or 
even  less. 


In  order  to  explain  this  question,  we  calculated  function  V  also 
with  the  aid  of  the  high-altitude  factors,  expressed  directly  through 
the  functions  of  parabolic  cylinder  [formula  (1.10),  (1.11)],  but  not 
through  their  asymptotic  representations.  In  this  case  of  pole  of 
integrand  F(t,  y,  y’)  they  are  the  roots  of  equation  [compare  (1.02) 


and  (1.10)] 


Si  (so)  —  to* 


where  for  brevity  through  5,  is  marked  the  value  *,  which  corresponds 

y*o. 
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The  remaining  factors,  entering  f^O,  t),  zero  do  not  become. 
Deduction  0  is  calculated  from  the  formula  • 


Rm  o 


*  *r.)/i  (y'.  MO- 


(2.02) 


[see  also  (1.10)].  Value  $  is  found  from  the  transcendental  equation 


sh  2«  —  2u  =  IS  (y)  —  Sol.  (2.03) 

which  is  obtained,  if  ve  calculate  integral  (1.13)  with  the  aid  of 
the  substitution 
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- 

O') 

v  ■ 

l 

-  i 

» 

t*  \ 

A tn  no 

m>.ic  {1. 23,  | 

I 

AJ„.  no  4<op- 
;  kv.ic.J1. 23} 

■i  i  f  “ 

t- 

-:.CK-ie.045 

/<  .2.-; 

*  '  •} 

■ —* —  x  *• 

— (  .  *25* — . 

— (•  x‘i.466 

— 0.107— .•* 

- 

t.'iS - ...‘I? 

— t-.12fr-.M2. 

— ttvGST  * 
! 

— C.149  ~i ‘>.26? 

—C'.lSc*—  {(‘.2*v 

Key:  (1).  Number  of  pole.  (2).  according  to  formula. 
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The  necessary  root  of  equation  is  selected  from  the  condition  that  a* 
the  real  value  of  t<p[y,)  value  5  must  be  negative  real  number,  and 
from  the  considerations  of  continuity.  For  calculating  functions 
(1.11)  we  used  the  series/row 

rr. 

r,  *■  _  V  •  ni  _  r»  \ 

Dr  - e  ‘  T  - -7 - : — -  z"'-  (2.05) 

r  i  i—  m  n: 

It  is  first  of  all  interesting  to  compare  the  roots  of  equation 
(1.23)  and  more  precise  equation  (2.01).  We  bring  in"jEable  1  of  the 


value 


=  —  pu„i 


[see  (4.02)],  obtained  according  to  formulas  (1.23)  and  (2.01).  Let 
us  note  that  precisely  value  A;,,  is  utilized  actually  during  the 
calculations  on  the  series/rcv  of  deductions  and  is  determined. 
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therefore,  the  accuracy  of  results.  Therefore  we  bring  in  all  tables 
instead  of  very  poles  lpt  of  value  A/m.  ^bout  the  selection  of 
values  Y  see  paragraph  4).  The  coincidence  of  values  A/m.  of  those 
calculated  according  to  formulas  (1.23)  and  (2.01),  especially  their 
alleged  parts,  can  be  considered  satisfactory. 

Fig.  ltft  and  lb  gives  curves  1  for  the  attenuation  factor, 
calculated  according  to  asymptotic  formulas  (1.23)  and  (1.24)  (thin 
lines),  also,  with  the  aid  of  the  functions  of  parabolic  cylinder  of 
formulas  (2.01)  and  (2.C2)  (thick  lines). 

FOOTNOTE  x.  On  all  drawings,  carried  out  on  the  logarithmic  scale, 
along  the  logarithmic  axis  are  plotted  decimal  (but  not  natural) 
logarithms.  ENDFOOTNOTE. 
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Fig.  1.  Attenuation  factor,  calculated  according  to  the  formulas  with 
the  functions  of  parabolic  cylinder  (thick  line  and  thick  dotted 
line)  and  of  the  asymptotic  formulas  (thin  line  and  thin  dotted 
line):  a)  the  I  series;  b)  the  II  series,  ^fable  3. 

Page  366. 

From  the  comparison  of  these  curves  it  is  evident  that  both  methods 
give  even  closer  results,  than  it  was  possible  to  expect  on  the  basis 
Table  1,  especially  for  Y=23.11.  Therefore  the  following  presentation 
is  based  on  calculations  according  to  the  asymptotic  formulas  as 
siiqiler,  with  exception  of  one  case,  examined  in  paragraph  3. 

Analogous  calculations  we  conducted  also  for  the  case, 
investigated  in  Chapter  16,  moreover  it  turned  out  that  the 
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asymptotic  formulas  were  applicable  there  approximately/exemplari ly 
with  the  same  accuracy,  as  during  the  intra-layer  propagation. 

3.  Calculation  of  high-altitude  factors  with  the  aid  of  the  Airy's 
functions. 

In  paragraph  8  of  Chapter  14  the  derived  approximation  formulas 
for  the  high-altitude  factors  of  the  weakly  damped  "seized"  waves, 
for  which  the  pole  corresponding  to  them  is  arranged/located  between 
values  r  (_>,,>  and  p(0);  high-altitude  factors  are  expressed  through 
the  Airy's  functions  u(x)  and  v(x).  In  our  case  to  the  "seized"  waves 
correspond  the  first  of  pole  for  the  curves  with  Y>=208.01  and 
Y=109.20  (see  Table  4  and  Fig.  5  and  6),  and  deductions  in  them  can 
be  calculated  according  to  the  formulas  indicated. 

The  use  of  these  formulas  is  caused  by  the  fact  that  for  some 
values  of  y,  for  example  y  -  y-  2,  value  p(y)-tx  for  the  mentioned 
cases  proves  to  be  small  (see  Fig.  5  and  6)that  is  made  itself, 
impossible  the  use  of  asymptotic  formulas  of  paragraph  1.  However, 
the  use/application  of  series/row  (2.05)  in  this  case  is  difficult, 
since  in  the  cases  of  value  $  in  question  they  are  great  and 
series/ row  converges  very  slowly. 


The  approximation  formulas  for  the  high-altitude  factors  of 
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weakly  damped  waves  in  case  {/<«/,■  interesting  us  take  the  form 


where 


My.  0  =  2c1(v)e-»vi| hfe,  /),  \ 
f-AU-  t)  =  2c.,(v)e-nv\}>;{£/1  t),  ( 


(3-01) 


v.(»-  o=|'  [»<u +xe“'“  M .  | 

*  <*•  o  -  Vr^m  [» <fe>  r  «*”«  ft.)]  ■ 


(3.02) 


but  ^  is  determined  from  the  relationships/ratios 


n-M'O  dy=A-£!2 


(y  >  i/i). 

f  I  P  {y)~<-  dy  =  4-  (— hf"-  (y  <  «,). 


(3.03) 


moreover  through  ya  it  is  designated  less  of  the  roots  of  the 
equation 


p(y)  —  i  =  o. 


(3.04) 
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The  pole  of  integrand  in  this  case  are  the  roots  of  the  equation 

v(y  --f  e^Kfo)  =  0,  (3.05) 

where  through  is  marked  the  value  £,  which  corresponds  y=0. 


Finally,  deduction  Ri  is  equal  to 
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°jg  ~'nV,(!f.  <i)  (!/'■  <i>^j(Q-  <i» 


(3.06) 


High-altitude  factors  can  be  represented  vith  the  aid  of  the 
Airy's  functions  when  value  eiav  is  low.  This  condition  is  satisfied, 
if  Rev  is  negative  and  is  not  small  (-Rev  >  l).  The  latter  occurs 
for  poles  t!  with  Y=208.01  and  Y=109.20  C^able  2),  therefore,  this 
method  for  calculating  the  deductions  in  these  poles  is  applicable. 


It  is  possible  to  show  that  if  is  great  and  is  negative,  then 
formula  (3.06)  for  remainder  converts/transfers  in  (1.24).  Thus, 
formula  (1.24)  is  knowingly  applicable,  if  is  a  large  negative 
number  (or  it  lies/rests  at  certain  sector  around  the  negative  real 
semi-axis  and  great  in  the  absolute  value).  But  if  £*  is  small  or 
positive,  then  asymptotic  formula  (1.24)  becomes  unsuitable. 


This  is  confirmed  by  Table  2.  Both  with  Y=208.01  and  with 
Y=109.20  to  values  y  -  y.ib  correspond  v  and  £w  having  sufficiently 
large  negative  real  parts,  and  the  deductions,  obtained  according  to 
formulas  (1.24)  and  (3.06),  prove  to  be  close.  To  values  y  — 1/,/2  with 
the  same  v  correspond  small  h  (but  with  Y=109.20  real  part  £j  is 
already  positive),  and  deductions,  calculated  according  to  formulas 
(1.24)  and  (3.06),  sharply  they  differ.  In  this  case  the  deductions 
must  be  calculated  according  to  formula  (3.06),  and  formula  (1.24)  is 
unsuitable. 
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Thus,  if  value  p(y)-t  is  low  (£,  is  small),  to  calculate 
high-altitude  factors  according  to  formulas  (1.18)  is  impossible,  and 
it  is  necessary  to  use  formulas  (3-01)— (3.03) .  Let  us  note  that  the 
roots  of  equations  (3.05)  and  (1.23)  virtually  coincide. 

The  use/application  of  Airy’s  functions  makes  it  possible  also 
to  investigate  the  height/altitude  of  the  layer,  in  which  is 
propagated  the  seized  wave.  As  can  be  seen  from  formulas  (3.03),  real 
part  is  negative  only  with  y<y,;,  where  yx  is  a  smaller  root  of 
equation  (3.04).  With  y>y*  value  Re  h  becomes  positive  and  with 
increase  in  y  increases  (value  Im  always  remains  very  small,  since 
for  the  seized  waves  the  alleged  part  tx  is  negligible).  But  with  the 
large  positive  arguments  function  v(|i)  exponentially  decreases,  and 
function  u(^x)  although  increases,  cannot  become  more  than 
e-Jtv  (v  <  0)  and  therefore  it  is  extinguished  by  factor  e:,u  (compare 
equations  (3.03)  and  (1.12)  and  the  asymptotic  representations  of 
Airy’s  functions  in  addition  2).  Consequently,  function  and 

means,  and  deduction  Hi  after  transition  through  y=yj  rapidly 
decreases. 
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The  dimensionless  height/altitude  it  is  possible,  on  the 

strength  of  what  has  been  said  above,  to  name/call  the  effective 
height/altitude  of  this  seized  wave,  which  corresponds  to  pole  t», 
and  in  the  first  approximation,  to  consider  that  the  field  of  this 
wave  occupies  the  layer 

0  <y<yt.  (3.07) 

but  out  of  this  layer  it  is  negligibly  small.  In  reality 
height/altitude  ysyx  is  not  the  completely  clear  boundary  of  layer, 
but  near  this  height/altitude  occurs  the  smooth,  but  rapid  weakening 
of  the  high-altitude  factors  (and  thus,  the  fields  of  wave). 

The  noncaptured  waves  ( y a  is  composite)  do  not  possess  the 
property  indicated,  and  their  field  is  distributed  in  entire  layer  of 
inversion  and  even  it  is  higher. 

4.  Intra-layer  propagation  of  waves  (numerical  results  in  the 
dimensionless  coordinates). 

Attenuation  factor  V  depends  on  the  dimensionless  coordinates  x, 
y,  y'  and,  furthermore,  from  function  p(y),  which,  in  turn,  is 
determined  by  parameters  //.  and  //,.  It  is  convenient  to  conduct  the 
still  dependent  parameters,  characterizing  function  p(y),  namely: 
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y  -  u,  ~  Hi- 

P  (!<.)  =  2v;  —  {),.  (-1-02) 

P  (0|  —  p  (ij;)  =  -pf  -  (4.03) 

Calculations  were  conducted  for  two  series  of  the  values  of  the 
parameters  of  function  p(y) ,  given  in*~£able  3. 
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Keys  (1).  according  to  the  formula. 
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Table  3- 


Hoitr 

Kpimoi. 

b: 

I'd-,) 

P  (0)  —  P  (.1/,-) 

i  \ 

20S.il! 

10.40 

1  cepiis  /j-i 

197.61  :  215.41 

0.542 

2  ! 

KWi.ofi 

5.09 

!  93.01' 

:  105.00 

0.260 

3 

45.0? 

2.40 

|  45.6? 

'  50.46 

0.125 

4  ‘ 

23.1! 

1.16 

|  21.95 

;  24.27 

0.06u 

i  \ 

109.20 

10.40 

11  CCpiifi  (3) 

98.80  1 19.66 

1.095 

0  ? 

52.50 

5.00 

\  47.50 

i  57.5o 

0.526 

3  | 

25.24 

2.40 

!  22.84 

|  27.64* 

0.252 

Key:  (1).  Number  of  curve.  (2).  series. 
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Curves  p  <;.•)  —  p  {ut\  are  plotted/applied  in  Fig.  2a  (I  series)  and  2b 
(II  series),  moreover  curves  are  marked  by  the  numbers  under  which 
the  parameters  corresponding  to  them  are  given  in  the  tables. 

The  parameters  of  function  p(y)  have  selected  we  so  that  with 
this  M-prof ile/airfoil  it  is  possible  to  calculate  attenuation  factor 
for  four  wavelengths,  which  relate  as  1:3:9:27  (I  series)  or  three 
wavelengths,  which  relate  as  by  1:3:9  (II  series). 

Let  us  note  that  the  I  series  of  the  parameters  was  accepted 
also  in  chapter  16  where  was  examined  the  case  when  one  of  the 
corresponding  points  is  located  highly  above  the  layer  of  inversion, 
and  other  -  within  the  waveguide.  For  the  latter  we  took 
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y  =  yt;bt  where  y>  is  spot  height  of  inversion. 

In  present  chapter  the  attenuation  factor  V  is  calculated  for 
the  intra-waveguide  propagation,  moreover  the  heights/altitudes  of 
transmitter  and  receiver  are  everywhere  selected  equal  to  (y^y'). 
Calculations  are  produced  for  two  cases: 


Si 

5  ’ 

(4.04) 

Vi 
o  * 

(4.05) 

Attenuation  factor \v| is  depicted  as  the  function  of 
nondimens ional  distance  of  x,  moreover  along  the  axis  of  ordinates 
are  deposited/postponed  the  values  of  common  logarithm] Fig.  3 
gives  the  curves  of  the  I  series,  in  Fig.  4  -  the  curves  of  the  II 
series.  Unbroken  curves  relate  to  case  (4.04),  broken  lines  -  to  case 
(4.05).  The  curves  of  the  I  series  are  calculated  according  to  four, 
while  the  curves  of  the  II  series  -  on  three  members  of  the 
series/row  of  deductions. 
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Fig.  2.  Plotted  functions  p(y)  —  p(yfc  a)  the  I  series,  b)  the  II 

series,  |:able  3. 

Page  371. 

Table  4  gives  values  &in  [formula  f2.06)],  entering  all 
calculations  and  interesting,  furthermore,  also  fact  that  they 
characterize  the  location  of  poles  tm  relatively  curved  p(y).  With 
Re  At„ > 0  we  deal  concerning  "seized",  very  weakly  damped  wave. 

This  pole  is  in  the  first  curved  each  series  (Y=208.01  and 
Y=109.20,  see  Table  4,  and  also  Fig.  5  and  6),  that  also  explains  the 
striking  special  feature/peculiarity  of  upper  curves  in  Fig.  3  and  4. 
Instead  of  decreasing,  as  this  seems  natural  for  the  attenuation 
factor,  these  curves  increase.  In  fact,  as  can  be  seen  fronf^able  4, 
the  alleged  parts  of  the  first  poles  in  this  case  are  very  small,  so 
that  in  the  region  of  the  applicability  of  one  member  of  the 
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series/row  of  the  deductions  where  there  remains  only  "seized"  wave, 
exponential  factor  virtually  fading  does  not  give.  However,  an 
increase  in  function  V  occurs  because  of  factor  jl  [see  (1.22)]. 
The  field  of  this  wave  carries  cylindrical  character,  since  is 
proportional  not  to  1/x,  but  I'y'x.  .  More  demonstrative  there  will 
physically  be  here  not  function  V,  but  function  T,  by  the  specific 
relationship 

V=  2]  Hxe  <4-06» 

whence 

T=  V.  Rmeixs'*  (-5.07) 


[compare  1.22)]. 
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Fig.  3.  Fig.  4. 


Fig.  3.  Dependence  of  attenuation  factor  V  on  nondimens ional  distance 
of  x  (I  series). 

Fig.  4.  Dependence  of  attenuation  factor  V  on  nondimens ional  distance 
of  x  (II  series). 

Page  372. 

It  is  the  attenuation  factor  of  cylindrical  wave,  which  is  spread  in 
the  surface  tropospheric  waveguide.  Here  seemingly  there  is 
reflecting  layer  (its  height/altitude  is  approximately /exemplar ily 
equal  to  y  ^beyond  limits  of  which  the  seized  wave  barely  is  drawn 
through  (compare  the  end/lead  of  paragraph  2).  Therefore  amplitude 
Rm  of  this  wave  with  y^yA  can  with  increase  of  y  decrease,  as  it 
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takes  place  in  Fig.  3  with  Y-208.01  and  4  with  Y=109.20;  in  all  other 
cases  the  attenuation  factor  monotonically  increases  with  increase  in 

y- 

As  can  be  seen  from  able  4  (see  also  Fig.  5  and  6),  in  cases  of 
Y=208.01  and  Y=109.20  the  second  of  pole  also  have  small  alleged 
parts.  Thus,  the  second  simple  wave  for  these  cases  attenuates 
although  considerably  more  rapid  than  the  first,  it  is  nevertheless 
very  slow.  Therefore  the  oscillations  of  attenuation  factor,  which 
are  obtained  as  a  result  of  the  interference  of  the  first  and  second 
simple  waves,  for  long  do  not  disappear. 
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Fig.  7  and  8  well  show  damping  oscillations  of  a  function  T:  around 
the  almost  horizontal  asymptote,  which  is  the  first  term  of  the 
series/row  of  deductions  (4.07). 
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Let  us  note  that  in  the  case  when  one  of  the  corresponding 
points  is  located  highly  above  the  layer  of  inversion  (chapter  16), 
the  first  simple  wave  which  slowly  attenuates  and  provides  therefore 
hyperdistant  propagation,  it  is  excited  very  weakly.  Because  of  this 
it  determines  field  only  far  in  the  shadow  when  several  following 
simple  waves,  which  are  decisive  for  the  field  near  the 
horizon/level,  manage  to  damp.  However,  in  the  case  of  intra-layer 
propagation  the  first  simple  wave  (Yc208.01  and  Y=109.20)  not  only 
greatly  weakly  attenuates,  but  also  is  excited  with  the  large 
amplitude. 


Ff 


if 
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Fig.  5.  Fig.  6. 


Fig.  5,  Location  of  roots  .  corresponding  to  seized  and  uncaptured 
waves,  relative  to  graph/curve  puo-puj, i  (I  series). 

Fig.  6.  Location  of  roots  corresponding  to  seized  and  uncaptured 
waves,  relative  to  graph/curve  r  «i/>  -  r  (II  series). 


Fig.  7.  Fig.  8. 


Fig.  7.  Dependence  of  attenutation  factor  of  cylindrial  wave  'p  on 
nondimensional  distance  of  x  with  Y=208,01. 

Fig.  8.  Dependence  of  attenutation  factor  of  cylindrical  wave  t  cn 
nondimensional  distance  of  x  with  Y=109-20. 
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At  the  remaining  values  of  parameter  Y  there  are  no  "seized" 
waves  already.  However,  here  in  the  composition  of  fields  are 
sufficiently  weakly  damping  simple  waves,  and  field  falls 
considerably  slower  than  in  the  absence  of  inversion.  For  the 
comparison  we  give  in  Fig.  4  dot-dash  line  of  the  fields,  obtained 
under  the  assumption  of  the  absence  of  refractions  for  y=y'=l  and 
y=y ,!=0.48f  "£hey  go  Y=52.5  and  Y=25.24,  which  give  field  at  the  same 
heights/altitudes,  as  dct-dash  curves,  but  in  the  presence  of 
inversion. 

Simple  dotted  line  in  Fig.  3  and  4  noted  those  sections  of 
curves  in  which  used  by  us  the  members  of  the  series/row  of 
deductions  give  already  sufficiently  rough  results,  although  even  and 
accurate  it  is  qualitative.  These  sections  very  closely  approach  the 
vertical  lines,  which  note  for  each  curve  the  geometric  boundary  of 
light/world  and  shadow  without  the  refraction. 

5.  Numerical  results  for  the  specific  case. 

On  the  basis  of  the  dimensionless  curves,  examined  in  paragraph 
4,  it  is  possible  to  obtain  numerical  results  for  the  specific  cases 
of  the  radiowave  propagation  of  different  length.  For  an  example  we 
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selected  M-prof ile/airf oil ,  obtained  from  the  first  dimensionless 
functions  p(y)  I  and  of  the  II  series  (jfable  3,  the  I  series,  row  1 
and~]iable  3,  the  II  series,  row  1)  of  formula  (1.07)  (Fig.  9),  Fig. 

10  and  11  give  the  appropriate  attenuation  factors  V  for  the 
wavelengths  3.33  cm  (curves  I),  10  cm  (curves  2),  30  cm  (curves  3) 
and  90  cm  (curve  of  4  Fig.  10).  Height/altitude  h,  identical  for  the 
points  of  radiation/emission  and  reception/procedure,  it  is  indicated 
in  Fig.  10  and  11  in  the  meters.  Distance  s  is  given  in  the 
kilometers. 

M-prof ile/airf oil  I  has  the  same  height/altitude  of  inversion 
hit  as  M-prof ile/airfoil  II,  but  doubly  more  strong  inversion 
M  (0)  —  M  (hi).  Consequently,  we  can  compare  attenuation  factors  for 
these  profiles/airfoils,  since  Fig.  10  and  11  correspond  to  the 
identical  heights/altitudes  of  the  corresponding  points/items  and  to 
identical  wavelengths.  Thereby  we  will  obtain  representation  about 
the  effect  of  inversion  (at  the  same  values  /,(  and  X)  on  the 
radiowave  propagation. 
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Fig.  9.  Dependence  of  the  given  refractive  index  M  on  height/altitude 
h  (M-prof ile/airf oil  I  and  II). 

Page  375. 

For  curved  1  and  2  superstandard  range  proves  to  be  greater  in 
the  layer  with  the  more  strong  inversion,  as  this  it  was  to  be 
expected.  However,  for  curved  3  occurs  a  somewhat  unexpected  effect: 
at  short  distances  the  curves  barely  differ,  but  from  a  certain  point 
onwards  field  it  proves  to  be  more  strongly  in  the  case  of  weaker 
inversion.  The  physical  sense  of  this  result  we  will  examine  at  the 
end  of  paragraph  6. 
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Fig.  10.  Fig.  11. 

Fig.  10.  Dependence  of  attenuation  factor  V  on  distance  of  s  (in  km) 
for  wavelengths:  1  -  3.33  cm;  2  -  10.0  cm;  3  -  30.0  cm;  4-90  cm. 
Unbroken  curve  -  for  a  height  h=9.31  m.  Dot-dash  line  -  for  a  height 
h=23.27  m.  (M-prof ile/airfoil  I). 

Fig.  11.  Dependence  of  attenuation  factor  V  on  distance  of  s  for 
lengths:  1  -  3.33  cm;  2  -  10.0  cm;  3  -  30  cm  for  heights/altitudes 
h*9.31  m  and  h=23.27  m  (M-prof ile/airfoil  II). 

6.  Attenuation  of  waves  in  the  tropospheric  waveguide. 
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As  we  saw  above,  attenuation  factor  V  is  represented  in  the  form 
of  the  series/row  of  deductions  (1.22).  Complex  numbers  determine 
the  dependence  of  the  separate  terms  of  the  series/row  of  deductions 
(or  "simple"  waves)  on  nondimens ional  distance  of  x.  They  are  the 
roots  of  equation  (1.23)  or  more  precise  equation  (2.01).  We 
investigate  in  more  detail  the  dependence  of  numbers  tm  on  different 
factors. 

As  it  was  shown  in  the  paragraph  of  6  chapters  16,  equation 
(1.23)  approximately  is  reduced  to  the  simpler  equation 

S,=  (m- 4-).t.  (6-0!) 

Here  m  is  number  of  this  root  t  ~  im,  and  Sx  -  integral 

Si  =  f  1  p(y)  —  i  dy,  (6.02) 

o 

where  through  yx  is  designated  the  smaller  root-  of  the  equation 

p  (y)  —  t  =  o.  (6.0i) 

* 

Page  376. 

Formulas  (6.01)  and  (6.02),  as  it  is  easy  to  show,  remain  applicable 
when  the  roots  of  equation  (6,03)  are  composite,  if  we  by  yx 
understand  root  with  the  positive  alleged  part.  With  their  aid  it  is 
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possible  to  approximately  find  roots  corresponding  both  weakly  to 
the  damped  and  strongly  damped  simple  waves  in  series/row  (1.22).  We 
will  use  these  formulas  in  order  to  explain,  on  what  parameters 
depends  the  attenuation  of  simple  waves.  This  will  make  it  possible 
to  obtain  representation  about  what  parameters  of  inversion 
M-prof ile/airfoil  determine  first  of  all  the  superstandard  range  of 
radio  waves  and  under  what  conditions  occurs  hyperdistant 
propagation. 


At  first  glance  seems  that  the  basic  parameters  are 
height/altitude  inversions  and  complete  increase  M  <()>  —  M  Ui  \  in 
the  layer  of  inversion.  Actually/really,  they  determine  the  so-called 
"critical  wavelength"  >.m  for  the  m  simple  wave  in  the  surface 


tropospheric  waveguide: 


-7—  h.  !  2  j.U (O.i  —  M  («,)! 


(see  Bremmer's  book  [25]). 


In  reality  the  assumption  about  the  primary  meaning  of  these 
parameters  is  incorrect.  In  order  to  study  this  question,  let  us 
introduce  instead  of  y  and  p(y)  the  new  given  variable/alternating: 

u  h  .  ,  p  ti/i  —  p  <«,)  ,  AJ  (ft)  —  A!  tk.f 

2  .  0(21  =  4  -- - -  =  4  -tt-t- — Trrrr  *  (6-Uo) 

y.  h .  •  plO)  —  r «?/)  A!  (0;  —  A1  (/</>  ' 

so  that  function  q(z)  always  satisfies  the  relationships/ratios 

q  (0)  =  4,  q  (1)  =  0.  (6-06) 
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Instead  of  by  the  variable/alternating  t  let  us  introduce  value 

'  =  4FW^h-  (6.07) 

Then  equation  <6.01)  will  take  the  form 


? _ (—4.). 

j) 


(6.08) 


where 


k-h; 


G  =  ^\p(0)-P  &,)]  10-*  [Ai  (0)  -  M  (ft,)).  (6.09) 

and  Zj  is  a  root  of  the  equation 

q  (?)  —  -t  -  0,  (6.10) 

that  corresponding  yx. 


Page  377. 

For  this  M-prof ile/airfoil  with  the  inversion  critical 
wavelength  }.n  is  determined  from  the  relationship/ ratio 

!  ( ri  — ~  ]  i 

I  |  '(/  u)  dz  ='- - ^ —  .  (6.  i  1 ) 

•  1  o... 

where  Gm  is  value  of  G  when  /.  =  /.„..  If  we  assume 

1 

f)  =  |  l  q  u)dz,  (6.12) 

f 

for  the  majority  of  M-curves  value  0  will  be  close  to  unity  as  a 
result  of  formulas  (6.06).  Thus,  for  M-prof iles/airf oils  examined 
above  we  have  0=0.991  (I  series)  and  0=0.983  (II  series).  From 
formulas  (6.11)  we  obtain 
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{n:~±), 

i  Gn  ~ - j7 — — ,  (6.13i 

whence  with  the  aid  of  relationship/ratio  (6.09)  we  find 

M  2  [.U  (0)  —  .11  (6.14) 

~  T 

This  expression  differs  from  (6.04)  only  in  terms  of  the  presence  of 
factor  p. 

Name  "critical  wavelength"  is  introduced  in  the  literature  for 
following  reasons.  When  G'-G...  i.e.  when  ;Vj.  equation  (6.08)  has 
real  root  *  ^  lying  within  the  limits  0<t<4.  To  this  root  in 
series/row  (1.22)  corresponds  the  undamped  simple  wave.  When  G  G„ 
i.e.  when  equation  (6.08)  has  the  composite  root  which 

determines  the  damped  simple  wave. 

We  will  subsequently  use  term  "critical  wavelength", 
understanding  it  as  the  abbreviated/reduced  designation  of  expression 
(6.14).  However,  at  the  critical  wavelengths  it  does  not  occur  any 
qualitative  jump  (see  below) ,  so  that  this  designation  is 
conditional. 

Equation  (6.08)  can  be  written  in  the  following  form: 
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Hence  it  is  apparent  that  for  this  M-prof ile/airfoil  [and,  therefore, 
for  the  fixed/recorded  function  q(z)]  all  roots  r  must  depend  only  on 
relation  y.mr,.  Instead  of  this  relation  to  more  conveniently  use  its 
common  logarithm 

/=log'f  =  log ]/£.  (6.16) 

Thus,  if  different  M-curves  have  coinciding  functions  q(z),  then  the 
corresponding  roots  t  must  be  placed  to  one  and  the  same  curve 

t  =/(/).  (6.17) 

The  drawn  conclusions  utilize  equation  (6.01)  for  roots  ;m. 
Actually  during  the  calculations  we  this  equation  did  not  use,  since 
it  is  too  rough;  however,  the  values  t  obtained  of  the  more  precise 
equations  are  approximately  subordinated  to  these  laws.  Fig.  12  gives 
functions  q(z)  for  two  series  of  M-curves,  examined  by  us  above. 

These  functions  with  0<z<2  virtually  coincide,  although  M-curves 
themselves  noticeably  differ  (compare  Fig.  9),  Therefore  values  r  for 
all  calculated  earlier  roots  tm  (see  Table  4)  lie  down  near  certain 
common  curve  (6.17).  This  is  evident  from  Fig.  13,  where  the  light 
small  circles  plotted/applied  values  Im  r,  which  correspond  to  roots 
tm  in  Jable  4  (series  I),  and  by  black  small  circles  -  values  Im  r, 
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which  correspond  to  roots  tn.  to  those  refined  according  to  the 
formulas  of  paragraph  2  (see  Table  1);  triangles  noted  points  on 

“fables  4  (series  II). 


Fig.  12.  M-prof ile/airfoil  in  new  given  variable/alternating 
[function  q(z)]  for  I  and  II  series. 

Fig.  13.  Alleged  parts  of  roots  r  of  equation  (6.15)  depending  on  by 
variable/alternating  1  (6.16). 

Page  379. 


Numerals  about  the  small  circles  and  the  triangles  indicate  number  by 
the  curve  p(y)  (first  numeral)  and  number  of  root  m  (second  numeral). 


i 
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Thus,  for  instance,  3.2  indicate  the  second  root  for  third  by  the 
curve  p(y).  Let  us  recall  that  the  alleged  parts  of  numbers  and 
*<r.  determine  the  attenuation  of  simple  waves  in  series/row  (1.22). 

We  see  also  that  with  1=0,  i.e.,  when  it  does  not  occur 

any  qualitative  change  in  the  attenuation  of  the  m  simple  wave.  This 
fading  according  to  a  strict  theory  is  both  when  and  when 

/.>/.„  and  sufficiently  slowly  it  increases  with  an  increase  X. 
Because  of  this  hyperdistant  propagation  can  to  be  observed  (although 
into  a  somewhat  weakened  form)  at  the  wavelength  X.  by  an  order 
exceeding  the  greatest  "critical"  value  X,.  Actually/really,  already 
curve  2  in  Fig.  10  corresponds  to  wavelength  X,  larger  than  X . , 
although  hyperdistant  propagation  has  the  place  even  for  the  longer 
waves  ( curves  3  and  even  4 ) . 

If  function  q(z)  for  the  undertaken  M-prof ile/airfoil  strongly 
differs  from  the  functions,  depicted  in  Fig.  12,  then  shape  of  the 
curve  (6.17)  will  qualitatively  be  the  same,  but  numerical 
relationships/ ratios  will  be  completely  different.  Thus, 
assuming/sett ing 


we  with  n=l/2  and  n=l/5  we  come  to  the  cases,  examined  in  works  [23] 
and  [24].  These  functions  q(z)  are  depicted  in  Fic.  14.  They  strongly 
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differ  from  the  functions,  plotted/applied  in  Fig.  12,  and  the  values 
r  for  them  with  the  same  1  are  obtained  approximately/exemplarily  by 
an  order  less  than  in  Fig.  13.  Latter/last  confirmation  follows  from 
the  figures  of  Hartree's  article  with  co-authors  [24],  if  we  consider 
che  relationships/ratios  between  the  designations  of  Hartree  and  our 
designations. 


Thus,  Fig.  13  is  not  universal,  and  to  calculate  with  its  aid 
the  attenuation  of  simple  waves  with  any  form  of  M-prof ile/airfoil  is 
impossible;  hypothesis  about  the  fact  that  the  hyperdistant 
propagation  is  determined  only  by  the  height/altitude  of  inversion 
hK  and  by  increase  M  (0)  —  Mi(ht),  proves  to  be  invalid. 
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Fig.  14.  Functions  q(z)  from  (6.18)  with  n=l/2  and  n=l/5. 

(M-prof iles/airfoils,  accepted  in  [23]  and  [24]). 

Page  380. 

Hence  it  follows  that  there  is  at  least  one  additional  parameter, 
which  characterizes  M-curve  and  which  has  fundamental  importance  for 
the  hyperdistant  propagation.  As  this  parameter  is  natural  to  take 
curvature  of  M-curve  at  the  inversion  point,  i.e.,  value  M'  </;,).  since 
this  parameter  substantially  affects  the  infiltration  of 
electromagnetic  field  from  the  layer  of  inversion. 

If  we  proceed  from  the  assumption  that  value  m *  (h,)  is  also  the 
basic  parameter,  then  instead  of  values  q(z)  and  r  should  be 
introduced  the  new  given  values 
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where  for  the  hyperbolic  law  of  the  inversion 

<?"(]>  =  6y-,  ( 6.20 1 

but  for  power  law  (6.18) 

<?' <1)  -  4?z.  (G.2i) 

In  Fig.  15  the  curve  A  corresponds  to  the  center  lines,  near 
which  are  grouped  the  values  t,  given  in  Fig.  13.  The  curve  B  is 
constructed  according  to  the  values,  undertaken  from  Hartree’s 
article  [24]  (for  the  function  (6.18)  with  n=i/5). 

Comparing  between  themselves  the  curves  A  and  B  in  Fig.  15,  we 
see  that  two  completely  different  types  of  M-prof iles/airfoils 
(compare  Fig.  12  and  14)  give,  if  we  use  given  values  (6.19),  not  too 
differing  values  of  the  attenuation  factors  of  simple  waves  in  the 
surface  tropospheric  waveguide.  Therefore  given  values  (6.19)  are 
more  convenient  for  the  comparison  of  the  properties  of  different 
layers  of  inversion,  than  value  (6.05)  and  (6.07).  This  result 
confirms  the  hypothesis  made  above  about  that  which  M'i/'m  is  also 
the  basic  parameter,  which  are  determining  the  infiltration  of 
electromagnetic  energy  from  the  layer  of  inversion.  Let  us  note  that 
real  parts  T  for  two  types  of  inversion  examined  are  distinguished 
much  more  strongly. 

Since  according  to  series/row  (1.22),  the  reduction  of  the 
amplitude  of  the  m  simple  wave  is  determined  by  the  factor 
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e 

where  s  is  a  horizontal 
the  more  demonstrative 
gives  the  coefficient 


=  e-'<  (6.22) 

distance  be  tv-*  :'.e  corresponding  points, 
representation  ioout  the  attenuation  of  wave 


y-n 


—  10  (h,)e. 


where  value  6  is  equal  to 


(6.23) 
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Pig.  15.  Comparison  of  the  alleged  parts  T  (6.19)  for 
M-prof iles/airfoils  of  different  types  (A  -  type  of  the 
profile/airfoil  Fig.  12;  B  -  type  of  the  profile/airfoil  Fig.  14). 
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This  value  is  plotted/applied  in  Fig.  16  as  function  from  1,  moreover 
the  curves  A  and  B  are  constructed  according  to  to  the  curves  A  and  B 

Fig.  15.  Formula  (6.23)  can  be  rewritten  as  follows; 

v...  '  ik  - — -  '  ]()""  t->  n-  •  *,< 

•*  '  f.  I  :*i  '.I  :-i  M  I  .• 

Since  in  formulas  (6.23)  and  (6.25)  on  wavelength  depends  only 
0,  then  Fig.  16  directly  indicates  the  character  of  the  dependence  of 
fading  on  the  wavelength.  First  with  an  increase  in  the  wavelength 
fading  increases,  but  then  it  reaches  maximum  and  begins  to  decrease. 
Tne  latter  is  easy  to  understand  physically,  since  on  the  propagation 
of  sufficiently  long  waves  the  layer  of  inversion  has  a  small  effect, 
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and  such  radio  waves  are  spread  further  beyond  the  hori2on/level , 
than  their  length  it  is  more,  as  it  takes  place,  also,  during  the 
normal  refraction. 

The  obtained  results  show  that  upon  consideration  of  three 
parameters  of  M-curv'e  i/i.. A!  i0) i/i,»  and  Art/i,))  the  attenuation 
factors  of  simple  waves  for  the  strongly  differing  types  of 
M-prof iles/airfoils  although  are  obtained  one  order,  they  differ 
quantitatively  (in  certain  cases  even  1.5-2  times,  see  Fig.  16).  This 
means  that  to  precise  account  for  the  effect  of  M-prof ile/airf oil  on 
the  propagation  it  is  necessary  to  take  into  consideration  even  some 
parameters.  This  question  needs  supplementary  investigation. 

Propagation  in  the  surface  waveguide  can  be  rated/estimated, 
comparing  attenuation  factor  x,  with  the  attenuation  of  the  first 
simple  wave  in  the  theory  of  radiowave  propagation  during  the  normal 
refraction.  This  comparison  makes  it  possible  to  explain,  how  this  a 
propagation  can  be  considered  "hyperdistant".  The  criterion  of 
hyperdistant  propagation  is  somewhat  indefinite;  however,  this 
uncertainty/indeterminancy  is  caused  by  the  very  nature  of  the 
phenomenon,  which  does  not  undergo  abrupt  changes  with  a  continuous 
change  in  the  wavelength  and  layer  of  inversion. 
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Fig.  16.  Dependence  on  1  attenuation  factor  0  of  simple  wave  (6.24) 
during  the  propagation  in  the  tropospheric  waveguide 
(M-prof ile/airfoil  of  types  A  and  B). 

Page  382. 

In  calculations  examined  above  in  the  second  series  of  M-curves 
is  undertaken  doubly  larger  value  m  tO)  —  A!  <  /:.»  at  the  same 
height/altitude  of  the  layer  of  inversion,  in  consequence  of  which 
increased  critical  wavelengths  (6.14).  Generally  speaking,  short 
waves  better  are  propagated  in  the  waveguides  with  the  larger 
critical  wavelength  (II  series),  than  in  the  waveguides  with  the 
smaller  critical  wavelength  (I  series).  However,  in  the  second  series 
increases  curvature  of  M-curve  at  the  inversion  point,  thanks  to 
which  is  facilitated  the  infiltration  of  electromagnetic  energy  of 
longer  waves  from  the  layer  of  inversion.  As  a  result  for  the  shorter 
waves  distant  propagation  proves  to  be  best  in  the  first  series, 
whereas  for  the  longer  waves  position  will  be  reverse/inverse  (curve 
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3  in  Fig.  11  corresponds  to  a  somewhat  worse  propagation  up  to  the 
greater  distances  than  curve  3  in  Fig.  10). 

During  the  calculation  of  field  in  the  layer  of  inversion  it  is 
most  convenient  to  use  the  series/row  of  deductions  (1.22). 
Physically  this  means  that  more  distant  and  hyperdistant  radiowave 
propagation  between  two  corresponding  points,  which  are  located 
within  the  layer  of  inversion,  can  be  considered  as  the  transmission 
of  cylindrical  waves  along  the  peculiar  transmitting  "line"  -  to  the 
layer  of  inversion.  The  attenuation  of  these  waves  is  caused  by 
radiation  losses,  their  amplitudes  depend  on  field  distribution 

according  to  the  height/altitude. 

The  results  obtained  above  illuminate  the  series/row  of  the 
laws,  which  relate  to  the  transmission  of  radio  waves  on  the  surface 
tropospheric  waveguide. 


